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Introduction 


Introduction 

We have in Ventus: Complex Functions Theory a-1 characterized the analytic functions by their 
complex differentiability and by Cauchy-Riemann’s equation. We obtained a lot of important results 
by arguing on line integrals in C. In this way we proved the Cauchy’s Integral Theorem and Cauchy’s 
Integral Formula. 

In this book we shall follow an alternative approach by proving that locally every analytic function 
is described by its Taylor series. Historically this was the original definition of an analytic function, 
introduced by Lagrange as early as in 1797. The advantage of this approach is that it is easy to 
calculate on series. The disadvantage is that this approach is not global. 

By combining the two aspects of analytic functions it is possible in the following to use Cauchy- 
Riemann’s equations , when they are most convenient, and series when these give a better description, 
so we can benefit from that we have two equivalent, though different theories of the analytic functions. 


Complex Functions Theory is here described in an a series and a c series. The c series gives a lot of 
supplementary and more elaborated examples to the theory given in the a series, although there are 
also some simpler examples in the a series. When reading a book in the a series the reader is therefore 
recommended also to read the corresponding book in the c series. The present a series is divided into 
four successive books, which will briefly be described below. 

a-1 The book Elementary Analytic Functions is defining the battlefield. It introduces the analytic 
functions using the Cauchy-Riemann equations. Furthermore, the powerful results of the Cauchy 
Integral Theorem and the Cauchy Integral Formula are proved, and the most elementary analytic 
functions are defined and discussed as our building stones. The important applications of Cauchy’s 
two results mentioned above are postponed to a-2. 

a-2 The book Power Series is dealing with the correspondence between an analytic function and 
its complex power series. We make a digression into the theory of Harmonic Functions , before 
we continue with the Laurent series and the Residue Calculus. A handful of simple rules for 
computing the residues is given before we turn to the powerful applications of the residue calculus 
in computing certain types of trigonometric integrals , improper integrals and the sum of some not 
so simple series. We include a residuum formula for the computation of the Mellin transform of 
some simple functions, and finally we show that the sum of some series can also be found easily 
by using Complex Functions Theory. 

a-3 The book Stability, Riemann surfaces, and Conformal maps is planned to be written soon. It 
will start with the connection between analytic functions and Geometry. We prove some classical 
criteria for stability in Cybernetics. Then we discuss the inverse of an analytic function and the 
consequence of extending this to the so-called multi-valued functions. Finally, we give a short 
review of the conformal maps and their importance for solving a Dirichlet problem. 

a-4 The book Laplace Transform will be the next one in this series. It will focus on this transform and 
the related 3 -transform, which in some sense may be considered as a discrete Laplace transform. 
Both transforms are of paramount importance in some engineering sciences. This book will be 
supported by examples in Ventus: Complex Functions Theory c-11. 
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a-5 and a-6 Future plans. The plan is then to continue with a book on Polynomials. Contrary to the 
common thought, the theory of polynomials is far from trivial. It is important, because polynomials 
are always used as the first approximations. Also, the topic Linear Difference Equations is of 
interest and far from trivial. However, the latter two books are postponed for a while. 

The author is well aware of that the topics above only cover the most elementary parts of Complex 
Functions Theory. The aim with this series has been hopefully to give the reader some knowledge of 
the mathematical technique used in the most common technical applications. 


Leif Mejlbro 
17th August 2010 
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1 Power Series 

1.1 Accumulations points and limes superior 

We shall later on need limes superior so we shall recall the definition from Real Calculus. 

Let (c„) be any real sequence. An accumulation point c € M of ( c n ) is a real number, such that for 
every e > 0 there exists an element c n from the sequence, such that \c n — c\ < e, or formally, 

Ve > 0 3n £ N : \c n — c\ < e. 

We extend for convenience this definition to also include the following cases, where we consider +oo 
(or —oo) as a (generalized) accumulation point of the sequence (c n ), if for every constant C > 0 there 
is an n £ N, such that c n > C (or c„ < —C), i.e. formally, 

VC'>03n£N:c„>C' for +oo, 

VC>03n£N:c Tl < —C for — oo. 
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Example 1.1.1 The sequence 

1, 1 ^ 1. 1, 2. 2, . . 3, 3, . . .... rij rt, . , .... 

2 2 3 3 n n 


has clearly the accumulation points — oo, 0, +oo. It is not hard to prove that when we include Too 
as possible accumulation points, then every real sequence has at least one accumulation point. In the 
present example we have got three accumulation points. 0 


If (c n ) —"► c for n —> Too, then the limit c is the only accumulation point. When c ± 00 , we say that 
( c n ) converges towards c. When c = Too or —00, we say that (c n ) diverges towards c. Notice that a 
divergent sequence does not necessarily diverge towards Too or —00. Two simple counterexamples are 
c n = (—1)" (a bounded, though not convergent sequence with the two accumulation points ±1) and 
c n = (—1 ) n n (an unbounded sequence, where Too and —00 are the two (generalized) accumulation 
points). 

We mention without proof the converse result. 


Theorem 1.1.1 Let (c n ) be a real sequence, where ( c n ) has only one accumulation point c. 

1) If c £ R, then (c n ) converges towards c for n —> Too, i.e. lim „^ +00 c n = c. 

2) If c = Too (or = — 00 ), then ( c n ) diverges towards c, i.e. lim n ^ +00 c n = c. 


Example 1.1.2 Usually a real sequence has many accumulation points. A very extreme example is 
the following. It is well-known that all rational numbers in the interval [0,1], say, are countable, so 
they can in principle be written as a sequence (<?„), q n £ Q D [0,1]. The countable set Q U [0,1] is 
dense everywhere in [ 0 , 1 ], hence every point in [ 0 , 1 ] is an accumulation point! 0 

Based on the discussion above we finally introduce 

Definition 1.1.1 Let ( c n ) be a real sequence. Then we define its limes superior, limsup „^ +00 c n , as 
the largest accumulation point c of ( c n ). 

If c el is finite, then for every e > 0 there are only finitely many n € N, for which c n > c T e, and 

infinitely many n € N, for which c — e < c n < c + e. 

If c— Too, then for every C > 0 there are infinitely many n, for which c n > C. 

If instead c = —00, then for every C < 0 only finitely many c n > C. 

Similarly, we can define limes inferior lim inf „^ +00 c n , as the smallest accumulation point c of the 

real sequence (c ra ), so 

lim inf c n = — lim sup { — c n } . 

n >+00 n —,-j-oo 

However, we shall not need lim inf „_,_(_ 00 c n in the following. 

It should be emphasized that the introduction of limes superior relies heavily on the usual ordering 
of R. For complex sequences, limes superior does not make sense at all. We shall only need lim sup 
to define the radius of convergence of the complex series in the following, and this only requires the 
lim sup of a real sequence. 
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1.2 Power series 

We shall typically deal with power series of the type 

-f-OO 

(1) ^ a n ( z ~ ZaT , Z £ C, 

n =0 

where we in general define ao (z — zq)° in Complex Functions Theory as oq. The coefficients a n are 
complex numbers, and the expansion point Zo £ C is fixed for all terms of (1). 

From the symbol (1) we define the corresponding sequential sequence of functions s n = s n (z), given 

by 

n 

(2) s n = s n (z) = ^2 a j ( z ~ z 0 Y , Z £C, 

3=0 

i.e. the n-th element s n (z) is the sum of the first n T lterms of (1). 

Definition 1.2.1 Consider the series (1) with its corresponding sequential sequence (2), and let Cl yf 0 
be an open set. We say that the series (1) converges towards the limit function f(z) for z £Cl, if 

n 

lim s n (z ) = lim > a, (z — zq) j = f(z) for all z £ Cl. 

n —»+oo n —^+oo ' ^ 

3=0 

The convergence of the series (1) is therefore derived from the corresponding sequential sequence 
(2). It must here be emphasized that the sequential sequence s n {z) = ^" =1 aj (z — zq) j must not be 
confused with the sequence (a n (z — zo) n ) neNo , which is obtained from ( 1 ) by just deleting the sum 
sign. Such a misunderstanding may cause some disastrous conclusions. 

We mention the well-known result that if a real series of continuous functions X]rS) /«( x ) bas a 
convergent majoring series Yln=o Cn < +°°j be. all c n > 0 are constants, and \f n (x) < c n for all 
relevant x, then fn( x ) is absolutely and uniformly convergent , and its sum function is continuous. 


We immediately extend this result to complex series of continuous functions, because we have 

m n ( Z )\ | +oo 

} < \fn(z)\ < C n and ^ c„ C Too, 

|S/„(2:)| J n—0 

and we can use the argument above on the real series 3? f n {z ) and J2n=o ^ f n {z ). 

We shall now more generally turn to the complex power series. Given (1), i.e. Y^n =o an ( z ~ z o) n , an( i 
consider the real sequence of the absolute value of the coefficients (|a„|). We introduce the number A 
by 

(3) (0 <) A := lim sup \/|a ra | (< Too). 

n—>■+oo 

Then we have the following theorem 
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Theorem 1.2.1 The power series Y^ n =o °n ( z ~~ z o) n absolutely convergent for every z G C, for 
which X\z — zq\ < 1, and divergent for every z G C, for which A \ z — ^ol > 1- 

Proof. 

1) If A \z — zo\ < 1, then 

limsup y \a n (z - z 0 ) n \ < 1. 

n— »+oo v 

It follows from the definition of limes superior that we can find a constant k G [0,1 [ and an N £ N, 
such that 

yj | a n ( z — -So) ,l | < k, thus \a n (z — Zo) n \ < k n for all n> N. 

Since k G [0,1 [, the sum Y2n=i v h n is convergent, hence a « ( z ~ z n) n is absolutely convergent 

for every such z € C, satisfying A \z — zq\ < 1. 

2) If instead A \z — 2 q| > then 


limsup y \a n (z - z 0 ) n \ > 1 , 

n —>-+oo v 


so | a n (z — Zo) n \ > 1 for infinitely many n G N, and the necessary condition , \a n (z — ^o)"l ~^ 0 , 
n —> +oo, for the convergence of ( 1 ) is not fulfilled. □ 


It follows from Theorem 1.2.1 that if 0 < A < +oo, then the power series Y^n=o an ( z ~ z o)” is 


absolutely convergent in the open disc T? ( ^ ’ anc ^ ^ divergent in the (open) complementary 


set C \ B 


z 0 , 


A 


of the closed disc B 


z 0 ■ 


A 


. It will be shown below in Example 1.2.1 that by the 


primitive test of Theorem 1.2.1 alone nothing can be said about the convergence/divergence of the 

power series on the circle | z — Zq \ = — , which separates the open domain of convergence from the 

A 

open domain of divergence. 


For completeness, if A = 0, then A \ z — Zq\ = 0 < 1 for all z G C, so the power series is convergent in 
all of C, and if A = +oo, then A \z — zq\ < 1 is only satisfies at the point z 0 , which is not an open set. 


The investigation above leads us to define the radius of convergence of the power series as the number 

(4) g := i = ---- —t==, 0 < g < Too. 

A hmsup„_ >+0O Vl«nl 

If g > 0, we call the open disc B ( zq , g) = {z G C | \z — zo| < g\ the disc of convergence. For conve¬ 
nience we say that B (zq, Too) = C is “a disc of radius Too”. 
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Example 1.2.1 On the circle of convergence \z — zo\ = g we do not get further information from 
Theorem 1.2.1. We mention with only sketches of proofs the following four (not exhausting) possi¬ 
bilities of convergence/divergence, where we for comparison in all four cases have chosen zo = 0 and 
0 = 1 - 


1) The series — 2 absolutely convergent for \z\ = 1. 

2) The series J2n=i divergent for |^;| = 1. 

3) The series — ~ n is divergent for z = 1, and it is conditionally convergent (i.e. the convergence 

depends on the order of the terms) for \z\ = 1 and z/ 1. 

4) For every aelwe let [a] £ Z denote the integer part of a, i.e. the largest integer n £ Z, for which 
n < a. The power series 


+00 . 

' n 


is conditionally convergent everywhere on the circle of convergence \z\ = 1. 

The former two examples are easily proved. In the latter two one has to apply Dirichlet’s criterion , 
known from real calculus. This is straightforward in 3), but difficult in 4). 0 
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If eventually all a n ^ 0 (e.g. for n > N), then it is sometimes easier to apply the following result 
instead of Theorem 1.2.1. 


Theorem 1.2.2 Given a power series J2n=o ttn ( 2 ~ z o) n > where a n =/= 0 for all n > N. If the quotient 
sequence 



is convergent, then it has the radius of convergence 


( 6 ) g = lim 

n —>-+oo 


&TI 

^n +1 


Remark 1.2.1 Notice that (5) may be defined without being convergent, and yet the series may have 
g > 0 (which then must be found e.g. by using Theorem 1.2.1 instead). One such example is 


+oo 


£{ 2 +(-!)"}" 


n—0 

where 

<2271-1 


<2277 

for n —> Too. 0 


for g = 


1 


= —[-* 0 and 

3 2n 


< 22 n 


<2277+1 


= 3 2n -» Too 


Proof. We consider the real series J2n=o l a «l ' \ z ~ z o\ n - Let z ^ zq, and write 


< 2 ti +1 


b n = \a n \ • I z — z 0 \ and X' = —- = lim 

g 71—7+00 

where we shall prove that g' = g , or, equivalently, A' = A. We get 

\z — zq \ —s► A' \ z — zq\ for n — > Too. 


^77 + 1 

1 < 277+11 • 

1 in+1 

\z-z 0 \ 

&n+l 


|<2tt. • 

hi 

0 

— 3 

1 



If A' \z — zo\ < 1, then choose k, such that A' \z — zq\ < k < 1. Due to the convergence there is an 
N £ N, such that 


Ti+l 


< k for all n > n > N , 


from which we conclude that 


0 < 6jv+ p < k ■ &jv+p-i < • • • < k p ■ for all p £ N, 

and since 0 < k < 1, the series 

+oo +oo 

bn = X! l a ”l ' I 0 “ z °i” 

n—0 n=0 
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is convergent in this case. 

If instead A' \z — zq\ > 1, then choose k, such that 
X' — ^o| > k > 1. 

There is an N G N, such that 

” +1 > k for all n > N, 

On 

hence 

bN+p > k ■ bw+p-i > • • • > k p ■ bN —> oo for p — > +oo, 

and the series ^ n c l ear 0 divergent in this case. 

The uniquely determined number A' satisfies the same condition as A in Theorem 1.2.1, hence X' = A, 
and thus g' = g. □ 

Example 1.2.2 Important! The simplest example of a power series, which is not a polynomial, is 
the geometric series 


+oo 

71=0 

In this case, all o„ = 1, so A = lirnsup r( ^ +OCJ ^/|a n | = 1, and g = — = 1 and zq = 0. Thus, the 

A 

geometric series is absolutely convergent, if \z\ < 1 and divergent for |^| > 1 , because then \z\ n > 1 
for all n £ N and the necessary condition of convergence is not fulfilled in this case. 

The geometric series is important, because it in some sense is the prototype of all power series of 
finite radius of convergence. We shall therefore find its sum function in the open disc \z\ < 1. 

More precisely, we claim that the sum function is 

f(z) = for |*| < 1 . 

1 — z 

In fact, by the usual algorithm of division we obtain 
1 

f(z) = - -= 1 + z + z 2 -\ -h z n +- - for \z\ < 1. 

I — z 1 — z 

The corresponding sequential sequence is given by 
$n(z) = lT T T • • ■ T z n -> 
and we see that 

( 7 ) I/O) - s n (z)\ = |/ 0 ) - {l + z + z 2 + ••• + z n } 


^n+l 


1-Z 


,|n+l 


< 


l-oi 
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ln+1 


Now, - —- —> 0 for n —> +oo, when \z\ < 1 is kept fixed. It therefore follows from (7) that 

1 - z 


( 8 ) 


1-z 


+oo 

= £•■ 
n =0 


pointwise for \z\ < 1. 


Let I\ C 5(0,1) be any compact set of the open unit disc. There is an r £ [0,1 [, such that also 
K C B[0,r]. Then we conclude from (7) for all z £ K that 


I f(z) -S n (z)| = | f(z) - {l + Z + Z 2 -\ - h Z n } \ < 


n+1 


1 — r 


0 for n —► +oo, 


so the convergence is uniform over every compact subset K of 5(0,1). Then f'{z) can be found by 
Theorem 3.4.2 in Ventus: Complex Functions Theory a-1 by termwise differentiation, i.e. 


1 

1-z 

2 


+oo 


= £*«> d -^2 = £ nz " 1 = £( n+1 > ra > 

n—0 ^ n= 1 n =0 

+oo +oo 

= n(n — 1 )z n ~ 2 = + 2)(n + 1 )z n , 


+OO 


+oo 


( 1 -*) 


n =2 


n —0 


k\ 


etc. for \z\ < 1, so the coefficients of . ., ,. 

1 (1 - z) k+1 

if a series is given by polynomial coefficients 


are polynomials of degree k in n. This implies that 


Pk{n) = a,kn k + • • • + a\n + do, for all n £ No, 
then the sum function of Pkinfz™ in B( 0,1) is a linear combination of 


k\ 


1-z’ 


(1-z) 2 ’ (1-z) 3 ’ 


_ z\k -\-1 


(i- 2 ) 


0 


Theorem 1.2.3 Let f(z) = Y^n=o a nZ n an d d( z ) = Y^n^o^nZ™ be two power series with the same 
expansion point zq = 0, and assume that they are both absolutely convergent for |;?| < r. Then the 
power series of their sum is given by 

+oo 

(9) (/ + g){z) = f(z) + g(z) = ^ (a n + b n ) z n at least for \z\ < r. 

n =0 

In some cases, (9) may be convergent in an even larger disc. 


The easy proof is left to the reader. That the sum may be convergent in a larger disc can be seen 
from the following example, where we choose 


+oo 

f(z) = J2 z 


n —0 


+oo 

and g(z) = — ^ z n both convergent only for \z\ < 1. 

n —0 
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Clearly, 


+OO 

(/ + S)(s) = f{z) + g(z ) = ^0./=0 for all * £ C. 

n =0 

Note that we have only proved that / + g = 0 in the disc \z\ < 1, but the strong property of being 
analytic implies that 0 is the unique analytic continuation to the largest possible set C. This shows 
that if we only argue on series and the situation is not as clear cut as the above, then we could get 
into some situations, where Cauchy-Riemann’s equations would be better to apply. 

The following theorem is difficult to apply in practice, and the unexperienced reader should avoid to 
use it. We shall, however, later on need a part of the proof, and it is furthermore quite naturally 
to show a theorem on multiplication, once we have obtained Theorem 1.2.3. Therefore, the reader 
should check the proof and is at the same time warned against using Theorem 1.2.4 in practice. Such 
applications are only for very skilled persons. 
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Theorem 1.2.4 Let f(z) = )C(5) a ra~" an d d( z ) = S)S) b e as * n Theorem 1.2.3, i.e. zo = 0, 

and they are both absolutely convergent for z < r. Then their product f ■ g has also a power series, 
and this is given by Cauchy multiplication, 

+oo 

(10) (/ ■ 9){z) = f( z ) ■ g( z ) = CnZn ’ 

n —0 

where the coefficients c n are given by the discrete convolution of the sequences ( a n ) and (b n ), which 
is defined by 


(11) c n := ^2 a k b n - k , for n £ N 0 . 
k =o 


Proof. It is given that 

limsup >/|a ra | = — < - and limsup y/\b n \ = — < -, 

n —*+00 Ql T n —>+oo Q 2 T 

so Qi, Q 2 > r. Choose any 0 < s < r. There exists a constant C = C s only depending on s, such that 

, C , „ . C 

\a n \ < ~ and b n < —. 
s n s n 

We shall first estimate (11), 


n n ^y 

— — 5Z ~gk ' 


c 

gn—k 


C^y^ 1= (n + 1)C 2 
S n ' ^ s n 


VfcJ < - \/n+ 1 • ^/C 2 , 

where lim„^ +00 f/n + 1 • •\/C 2 = 1, so we conclude that 
limsup ^/fcj < -. 

n —>+oo S 

This holds for all s < r, so we also have 
limsup ^ 

n —>+oo T 

and the series of the right hand side of (10) is indeed absolutely convergent for \z\ < r, hence 


+oo n 

(12) £5> fc H&n-fcH*r 

n= 0 k —0 

is convergent for |^| < r. 
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We note formally , if we collect the terms according to their power that 

+oo +oo +oo +oo 

f(z) ■ 9(z) = E a i z3 ' E bmZm = E E a J b mZ J+m 

j =1 m =0 j —0 m =0 

+oo ( n +oo 

= E i E akbn ~ k } = E CnZ,n > 

n—0 i. A;—0 J n—0 

and ( 12 ) shows that this formal series is absolutely convergent for \z\ < r. 

We shall now prove that the c„ given by (11) in reality gives the right series expansion of the product, 
and not just formally. 

We put 

/W(~) = ao + a,i z + • • • + OnZ N and 9n{z) = b 0 + biz + ■ ■ ■ + b N z N . 

Let \z\ < r. Then clearly, 

f(z) = lim fN(z) and g{z) = lim g N (z). 

N^+oo N— >+oo 

We have 

+oo n 

(13) \(fg) N {z) - f N (z)g N {z)\ < E Ei-mu-n", 

n=N-\-l k —0 

because all terms of degree < N have disappeared on the left hand side, and no term from ( fg ) n{z) 
enters the right hand side. 

Due to (12), for fixed z, \z\ < r, and every e > 0 there is an Nq £ N, such that the right hand side of 
(13) is smaller than e for every N > Nq. This shows that 

f(z)g(z) = lim f N {z)g N {z) = lim (/ • g) N {z) = (/ • g){z). □ 

N— »+oo N— >+oo 


The following important theorem contains a lot of information, much more than one would guess at 
a first glance. 

Theorem 1.2.5 Let ( z — z o) n be a power series of radius of convergence g > 0. Then the 

power series is uniformly convergent on every compact subset K C B (zo, g). The sum function 

+oo 

(14) f(z) = E a n {z - z 0 ) n , for z £ B (z 0 , o) 

n —0 

is analytic in B(zo,g), and the derivative f'(z) is obtained by termwise differentiation 

+oo 

(15) f(z) = E na « ( z - 2o) n_1 , for z £ B (zo, g). 

n =1 

The sum function is differentiable of any order p £ N with e.g. its derivative of order p given by 

+oo 

(16) / (p E) = E n ( n - 1) ■■■(n- p + T)a n (z - z 0 ) n ~ p , for z £ B (z 0 , g). 

n=p 
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Proof. Assume that K C B (zo, g) is compact, i.e. K is closed and B (zo, g) is open. Hence, there is 
an r £ [0, g[, such that K C B [zo, £>]. Then 

limsup v/| a n \ = - < -, 

n —>-+oo Q T 

and it follows in exactly the same way as in the proof of Theorem 1.2.4 that there is a constant C = C s 
corresponding to s £ ] r, g [, such that 

C 

\a n \ < — for all n £ N. 
s n 

If z £ K, then we get the estimate 


+ OO 


^ a n (Z - Z 0 ) T 


n —0 


+00 


n—0 


+00 


<E-- n =^IT=^ 

- S U ^ \ S J 


n =0 


because 0 < r/s < 1, so the latter series is convergent, and its sum is independent of the choice of 
z £ K , proving the uniform convergence. 

It follows from Corollary 3.4.3 in Ventus: Complex Functions Theory that (14) represents an analytic 
function of derivative (15). Finally, (16) is obtained by p successive termwise differentiations. □ 

A very simple application of Theorem 1.2.5 with an unexpectedly large effect is to put z = zq into 

(16) , in which case we only get a contribution from the term n = p. Thus, 

(17 ) f ip) (z 0 ) = p\ a p , i.e. a p = ^ / (p) (z 0 ). 

p\ 

Then by insertion of (17) into (14) we get 


+00 . 

f{z) = f (n) ■ ( z ~ z °)” ’ 

n=0 

and we have proved 


for z £ B (z 0 , g ), 


Corollary 1.2.1 Let f(z) be the sum function of a power series X])5) a « ( z — z o)” °f radius of con¬ 
vergence g > 0. Then f(z) is given by its Taylor series, expanded from the centre Zq, in B(zo,g), 

i.e. 

+00 . 

(18) f(z) = V — f {n) (z 0 ) ■ (z - z 0 ) n , for z£ B (z 0 , g). 
n\ 

n —0 


It follows immediately from Corollary 1.2.1 that we have 

Theorem 1.2.6 The Identity Theorem. Assume that the two power series 

+ OO +OO 

^2a n (z - z 0 ) and ^ A n (z - z 0 ) n , 

71=0 71=0 

expanded from the same point zq have positive radii of convergence and share the same sum function 
f(z) in their common domain. Then the two series are identical, i.e. a n = A n for all n £ N. 
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Proof. This follows directly from (17), because 

a n =~7 f (n) (z 0 ) = A n . □ 
n! 


Then we turn to the indefinite integrals. 


Theorem 1.2.7 Let f(z) be the sum function of a power series Y^n=o ttn ( z ~ 2 o) n . expanded from zq 
and of radius of convergence g > 0. The indefinite integral F{z) of f(z), for which also F (zo) = 0 , 
is in the disc B(zq,q) given by the termwise integrated series 


+oo 


(19) F(z) = 


n —0 


n + 1 


, {z - Zo) 


n+1 


Proof. It follows from the definition (4) that 


lim sup 

n —>-+oo 



lim sup 

n —>-+oo 





n —>+oo Q 


so the series of F(z) and f(z) have the same radius of convergence g > 0. They are both expanded 
from the same point zq, and it follows from Theorem 1.2.5 that F'{z ) = f(z), and the claim is proved. 


□ 
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Referring to Theorem 3.3.5 of Ventus: Complex Functions Theory a-1 we conclude that every indefinite 
integral f(z) in B (zq, g) has the structure F(z) + c for some uniquely determined constant c6C. □ 


Example 1.2.3 Using Theorem 1.2.7 on the geometric series 

indefinite integral F(z) of- in this disc, for which F( 0) = 

1 — z 


Y^n =o ~ n f° r \ z \ < 1 if follows that the 
0, is given by 


+oo 


nz) = E 


1 

„ n + 1 

n =0 

On the other hand, 
G(z) := Log 


+oo 


n+1 


= E 


z . 


n 


1-z 


= —Log(l - z) 


is also analytic in this disc, and we have 


1 


G'(z) = - -= f(z) 

1 — z 


for \z\ < 1, 


Hence, G(z) is also an indefinite integral of f(z), so G(z) = F(z) + c for some c £ C. Finally we see 
that F'(O) = G(0) = 0, so c = 0, and we have proved another important result, 


(20) Log 


1 - 


+oo .. 

Log(l - z) = E - z n , for \z\ <1. 0 

n 


We shall emphasize in this Ventus: Complex Functions Theory series that one must always specify the 
domain of convergence of a series, because otherwise one could easily jump to very wrong conclusions. 
It is of course legal to try to find a formal solution of a problem, but once a formal series solution 
has been found, one should immediately find the domain of validity, outside which the result is not 
reliable. 


1.3 Expansion of an analytic function in a power series 

We proved in Section 1.2 that the sum function f(z) of a power series expansion from zo and of radius 
of convergence g > 0 is analytic in B ( zq , g) and that f(z) in B ( zq , g) is given by its Taylor series 
expanded from the center Zq of the disc. Furthermore, Theorem 3.4.2 of Ventus: Complex Functions 
Theory a-1 showed that every analytic function is infinitely often (complex) differentiable. 


Remark 1.3.1 The situation is different for real functions in C°°(]R), because far from all of them 
can be extended to an analytic function by “just writing z £ C instead of x £ K”, a wrong statement 
which is frequently met. It is possible and even not too difficult to construct a real C°° function which 
cannot at any point Xo £ R be extended to an analytic function in any complex neighbourhood of 
Xq £ R. We shall give an example in Remark 1.3.2 where this phenomenon occurs in one point, from 
which this general result can be derived by some advanced, though standard mathematical procedure. 
0 
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We shall now show a converse result, namely that every analytic function locally is a sum function of 
a power series. Once we have proved this result, we have shown that every analytic function can be 
treated either by Cauchy-Riemann’s equations , or by local power series. 


Theorem 1.3.1 Let f : > C be an analytic function in an open domain f 1, and let zo £ be any 

fixed point. The Taylor series of f(z) expanded from zq is convergent in (at least) the largest open 
disc B (zq, g) C Ll of centre zq. In B (zq, g) the sum function of the Taylor series is f(z), thus 

+oo .. 

(21) f( z ) = £ -T / (ri) (- 0 ) • (2 - Z 0 ) n for all B (z 0 , g). 
n\ 

n —0 


Proof. It follows from Cauchy’s inequalities , cf. Ventus: Complex Functions Theory a-1 , Theo¬ 
rem 3.4.5, that 


/ (n) Oo) 
where 


• 77 ,! 

< —-- 1 for every n € No and r £ ] 0 , g[, 


M r = max{|/0)| | \z - z 0 \ = r}. 

If M r = 0, then the Taylor series is the zero series, which of course is convergent. 
If M r > 0, then the radius of convergence of the Taylor series is at least 


£+00 ^ M r n\ 


= r ■ lim — 7 - = r. 

n^+00 yM r 


This holds for every r £ ]0, g[, so we conclude that the Taylor series of f(z) has at least g as radius of 
convergence. 

Choose any r £]0,p[ and any point z £ B(zo,r). Then \z — zo\ < r, so if £ lies on the circle 
|£ — zq\ = r, then we have the estimate 


+ OO , y 

E {z ~ Zq) 

(/- „ 'in+l 


n^O (C - *>)* 


+ OO 




n —0 


\z - Zo\ 


r — \Z— Zq I 


from which follows that the series -+T uniformly convergent in £ for |£ — zo\ = r. 

(C - Zo) n 

Using (95) of Theorem 3.4.2 of Ventus: Complex Functions Theory a-1 we get 

/(C) 


' w <*>=£/, 


|C-Zo|=r (C ~ z 0) 


n+1 


d£ for n G Nq. 
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Insert the Taylor series, then interchange summation and integration followed by a reduction and then 
finally apply Cauchy’s integral formula to get 


+ OO 


+oo 




n—0 


n —0 


/(C) 


2tt* 7| C _ Zo | =r (C - z 0 ) 


n+1 


dC • (i 


*>)" 


1 

2lti 


/(C) 

IC-zol=r „=0 


+oo , >1 

K z ~ z 0) 


^0 (C - z o) n+1 d( 



1 

C - z 0 


1 

z - z 0 

c^ 


dC 


1 

2iri 


/■ 

■'IC—*ol=r 


d C =/(*)• 


Finally, notice that to every z G B(zq,q) we can choose r < g, such that z G I3(zo,r), where the 
computation above is valid, and the theorem is proved. □ 


Some basic power series. 


It follows from Theorem 1.3.1 that all known real Taylor series are immediately extended to complex 
Taylor series, because the Taylor series only depends on its sequence of coefficients, { f^ ( 20 )}, 


derived by differentiation. We therefore get the complex Taylor functions of the following well-known 
functions. The reader is highly recommended to learn all these by heart, as the appear over and over 
again in the following, as well as in applications outside these books. 


1) 

exp z = J2n= 0 ^ zn , 

zee, 

2) 

™ J = E -0(1,)! * 2 ”’ 

zeC, 

3) 


2G<C, 

4) 

cosh z — V +0 ° 1 * 2n 

^ n =° (2n)! ’ 

2G<C, 

5) 

sinh z-V + °° 1 z 2 n+i 

sum 2 — z^n=o (2n _|_ 1 ) ! ’ 

2 g C, 

6) 

Log(i+2)=£:r 0 ^^ +i , 

\ z \ < 1 , 

7) 

(i + *)“-EIS (“)*". 

N < !> 

8) 

1 y^+°° 

1 _ z ~^n=0 Z - 

|2| < 1. 


Formula 7) is strictly speaking the definition of what later is called the principal value of the (usually) 
multiply defined function (1 + z) a . If a = n G No, then (1 + z) n is of course a polynomial instead, so 
it is uniquely defined for z G C, and not multiply defined in this exceptional case. 
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The Taylor series l)-5) are found by means of Theorem 1.3.1, because it is in all these cases easy to 
find /M(0). 

Formula 6) is obtained from (20) by writing — 2 : instead of 2 and then change sign. We notice that 
Log(l + z) itself is defined in the open domain = C\] — 00 , — 1], so the largest open disc contained 
in of centre 0 is B( 0,1). 



Figure 1: The domain of the Taylor series of Log(l + z) expanded from zq 


0 is the open unit disc. 


Formula 7) is here considered as a definition of the principal value of (1 + z) a , where we define the 
general binomial coefficients by 


a. 

n 


a(a — 1) • • • (a — n + 1) 


a £ C, nt No, 


with n factors in both the numerator and the denominator. Notice that if a = n £ No, then the series 
of (1 + z) n is a polynomial of degree n, and the domain is all of C. 


Remark 1.3.2 Again the situation is different in the real case, C” 00 ^). It is not hard to construct 
a real ip £ C 00 ^) and a corresponding point Xo £ R, such that the (real) Taylor series of <p(x) is 
convergent everywhere in R, and such that 

+OO 

<p(x) / ^ (^0) • (x — Xo) n for every i£l \ {a?o} • 

n=0 

One simple example of such a function is 
1 


<p(x) = 


expl "w 


for x £ M \ {0}, 
for x = 0. 


Clearly, <p(x) is C°° outside x = 0, and for x = 0 we use the definition of a converging sequence of dif¬ 
ference quotients and one of the rules of magnitudes of functions (exponentials dominate polynomials) 
to prove that 


- ^(°) = 1 . exD 
x — 0 x 


for 
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and then by induction also for higher derivatives to get y/ n )(0) = 0 for all n £ No- In all cases we get 
<^( n )(0) = 0. Thus, the Taylor series is the zero series, and it is obvious that <p(x) > 0 for x yf 0. 0 


Theorem 1.3.1 shows that if a function f(z) is analytic in the disc B (zo,r), then the function f(z) 

is in this set alone given by the sequence j ^7 f^ (z 0 ) 1 . This sequence is sometimes called the 

l n! JneNo 

germ of the analytic function f(z) in B ( Zo,r ) expanded from Zq- 


Remark 1.3.3 From an application point of view it is strange that an analytic map f(z) is uniquely 
determined in a whole disc B (zo,r), if we just know its germ j — /(") (zq) 1 at the centre Zo¬ 


rn 


raG N 0 


This is again a warning to the reader that the analytic functions may be easy to handle in practice, 
but they have there limitations, and they cannot provide us with a universal model of the real world. 
In particular, it is annoying that they can never directly describe causality. 0 
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Remark 1.3.4 The real function f(x) 
series, 


1 

1 + x 2 ’ 


x G M, is of class C°°(R). 


However, its Taylor 


fix) 


1 

1 + x 2 


+oo +oo 

V (-*T = E(-i)"x 2 " 


n—0 


n=0 


is only convergent for |x| 


< 1. This looks like an enigma, as long as we only consider if 1. 



Figure 2: The largest open disc of centre 0 not containing the two singularities ±i. 


In the complex plane we see that f{x) is a real rational function, so it is uniquely extended to the 
analytic function 

f{z) = , * 2 , for 2 G C\ {±i} = SX 
1 i - z 

The largest open disc contained in C C of centre 0 is B( 0,1), so in the complex plane we see why 
the radius of convergence is only 1, when Zq = 0. The complex singularities ±z have therefore a 
profound influence on the convergence of a real Taylor series. This strange phenomenon has puzzled 
many students, who had no knowledge of Complex Functions Theory. 0 


We finally prove 


Theorem 1.3.2 Weierstrafi’s Double Series Theorem. Let { 5 , n ( 2 )}„ eNo be a sequence of functions 
which are all analytic in the same disc B( 0, p). Assume that the series 

+oo 

fi z ) = ^9n{z) 

n =0 

is uniformly convergent in every smaller closed disc B[0,r], r < g. Then we obtain the power series 
of f(z) by first expanding all the g n {z) and then collect all terms of the same power of z. 
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Proof. It follows already from Corollary 3.4.3 of Ventus: Complex Functions Theory a-1 that f(z) 
is analytic in 5(0, g) and that 


+OO 


(22) / (P) (0) = E^ P) (0). 


n —0 


Hence, for z £ 5(0, g), 


+oo 1 +oo f +00 1 

p =0 ** p—0 \ n —0 ** 


□ 


Example 1.3.1 Clearly 


cos z 


+°° l On 

_ \ ' y i l Jln 
^ (2n)! 

n=0 y ’ 


for z £ C, 


is continuous, so there is a g > 0, such that 
11 — cos | < 1 for ^ £ 5(0, g). 

If we put 

9n(z) ■= (1 - cos z) n , 

then we can in principle find the power series of g n (z) by Cauchy multiplication. Choosing r < g we 
get in 5(0, r) that 


1 


1 


+oo 


+oo 


COS z 


1 — (1 — cos z) 


= - C0S -) n = ^29n(z), 


n —0 


n —0 


so using (22) we can in principle find the Taylor coefficients of-, i.e. its germ. We therefore obtain 

1 COS 0 

a power series expansion of- in the disc 5(0, r). 


cos z 


We shall not go into details with the sketch above, because the computations are fairly big and 
extremely tedious. We shall, however, point out the following unexpected result: Since - is 


cos z 


analytic in the set C \ j — +p7r p £ z|, it follows from Theorem 1.3.1 that the Taylor series is 
convergent in the disc 5 ^0 , — ^. This is far from trivial, because 


(r) 


cos ( — *) — 1 


> ^ >1 ’ 


for all 2 £ B(0,g). 


It follows by the continuity that the chosen g above, which was used to compute the Taylor series, 

7r 

must satisfy the inequality g < —, and yet the final domain of the Taylor series constructed on 5(0, g) 


is the larger set 5 ^0 , ^ j ■ 0 
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1.4 Linear differential equations 

One of the big issues of the theory of power series is the solution method of finding an analytic solution 
of a linear differential equation with polynomial coefficients. 

Existence and Uniqueness Theorems. 

We start with the following theorem of existence and uniqueness of the solution. 


Theorem 1.4.1 Let Ll be an open subset of C, and let ao(z), ..., a n (z) and g(z) be analytic functions 
in fi. Let zq £ Ll be a point for which ao ( 20 ) 7 ^ 0. 

For any given complex numbers Cq, c±, ..., c n -i, there exists one and only one function f(z), which 
is analytic in a neighbourhood u> of zq, such that 

d n f df 

(23) a 0 ( 2 ) — H-h a n _i{z) — + a n (z)f(z) = g(z) for all z 

and such that 

(24) / (zo) = Co, f'(z 0 ) = c t , ..., / (n ~ 1} ( 20 ) = c n -i. 

One may choose lo as the largest open disc B (zq, g) C f l, which does not contain any zero 0 / 00 ( 2 ). 


Sketch of proof. The equation is only considered in B (zo,g), where 00(2) 7^ 0 . To ease matters 
we norm the differential equation (23), which means that we divide it by ao(z), so that the coefficient 
of the highest order term is 00 ( 2 ) = 1. 


First assume that f(z) is indeed a solution in B ( 20 , g). Let 0 < r < g. By Cauchy’s inequalities (cf. 
Theorem 3.4.5 of Ventus: Complex Functions Theory a-1) there exists a constant M r , such that for 
all j = 1, ..., n, and all k £ N 0 , 


(25) c4 fc) (2 0 ) 


< 


M r k\ 


and 


9 (k) (zo) 


M r k\ 

rjr*lc 


Using that 00 ( 2 ) = 1 we get by a rearrangement of (23) that 


f { n \ z ) = 5(2) -^ 2 a n - j ( z ) f u \ z ), 
j =1 

thus by k differentiations, 

n —1 k / 7 \ 

(26) f {n+k) (z)=g^(z)-Y J Y,[ a ) f U+q) (z). 

3 =1 9=0 V q ' 

Using that 

f(zo)=co, f(z 0 ) = c 1 , ..., / (n ^ 1} (20) = C n _i, 

are given, we find /^ (zq) for all p € N 0 . 


It follows that if an analytic solution exists, then it must be unique, because its Taylor coefficients 
are uniquely determined. 
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It is seen by insertion of the formal Taylor series into (23) that if the series is convergent, then it must 
be a solution. 


The difficult part is to prove the existence. The idea is again to apply Cauchy's inequalities. We 
assume that 


(27) 


f {q) (zo) 


<C q 


j'Q 


for q = 0, 1, ..., n + p — 1, 


where the C q are given constants. Then proceed in the following way (the complicated proof is left to 
the interested reader). Put (25) and (27) into (26) to find constants C n+P as small as possible, such 
that 


(28) 


f {n+p) (zo) 


— C n ^-p 


(n + p)\ 

r n+p 
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The sequence ( C q ) is clearly increasing. What is more important, however, is that it is bounded 
from above. (The difficult proof is left to the interested reader.) Hence, there exists a function C(r), 
depending on r < g, such that 


f (n) (zo) 


< C(r) ■ 


n\ 


for all n. 


Then it follows from Cauchy’s inequalities that the Taylor series in convergent in B (zo, t). 


This procedure can be performed for every r < g, so we conclude that the Taylor series is convergent 
in at least the disc B (zq, g), and the existence follows. □ 


Remark 1.4.1 It should be mentioned that the Taylor series in some cases (though not in all) can 
be extended over (some of) the zeros of the coefficient ao(z), in which case the radius of convergence 
of the Taylor series becomes larger than g. 0 


Corollary 1.4.1 Every linear and homogenous differential equation (23) withg(z) = 0 and of analytic 
coefficients has in a neighbourhood u of every point zq € for which a o (zo) 7 ^ 0 , precisely n linearly 
independent solutions. 


PROOF. This follows immediately from the fact that every solution is uniquely determined by the n 
constants Co, c±, ..., c n _ 1 of (24). □ 

Practical procedures for solving a linear differential equations of analytical coefficients. 


In this subsection we shall in some examples demonstrate three standard procedures os solving a linear 
differential equation of analytical coefficients. These are 

1) Inspection 

2) Calculation of the germ { f^ (zo)}„ eN 

3) Method of power series. 

Of these, inspection is the most difficult one. However, when it succeeds, it is also the most elegant 
method. It requires some skill in manipulation. 

The method of calculation of the germ may not always be applicable, but when it succeeds, it is 
usually straightforward. 

The method of power series is the most commonly used method, because it is easy to understand, 
and for the novice it is felt as the procedure. Shortly described, one inserts a formal power series 
f(z ) = J2n=o^ n ( z ~ z o) U into (23), where the constants b n are the unknowns. Collecting the terms 
of the same power of the result we obtain a recursion formula (or a difference equation) in the b n , 
which then is solved. 

Finally - and this is very important - the constructed series is still formal, so we must always finish 
the task by computing the radius of convergence. In fact, the formal series may in some cases be 
divergent for every z with the exception of the expansion point zq itself, because g = 0 . 
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Referring to Theorem 1.4.1 we see that g = 0 is possible, when cto ( zq ) = 0, in which case we cannot 
conclude anything about the existence of a convergent power series solution. One may insert a formal 
series and obtain a result, but then we must check if it has a positive radius of convergence. 

Remark 1.4.2 We mention without proof that the zeros of ao(z) determine all the possible radii 
of convergence, so one can by knowing these make a qualified guess of g, which is either one of the 
numbers \zj\, where the z :j are all the zeros of ao(z), or +oo. Note also, that the number of zeroes 
could be infinite, e.g. for ao(z) = sin z, cos z, sinh;z or cosh z. 0 


Example 1.4.1 We shall demonstrate the three methods on the simple equation 
(29) f’(z)-f(z)= 0. 

First we see that according to Corollary 1.4.1 there is, apart from a constant factor, just one solution. 


Then use Remark 1.4.2. Since ao(z) = 1 does not have any zeros in C, we may expect that the radius 
of convergence is g = +oo, which by Remark 1.4.2 is the only possibility. 

First method, inspection. By checking our arsenal of known common analytic functions we imme¬ 
diately see that 



so f(z) = e z satisfies the differential equation (29). We therefore conclude by the beginning of this 
example that the complete solution is given by 

(30) f(z) =c-e z , ceC. 


Variant of the first method. We may instead multiply (29) by a so-called integrating factor. By 
this method we shall use the well-known rules of calculations, 

(31) f'{z ) • g(z) + f(z ) • g'(z) = ■ g) and f(z ) • g(z) - f(z) ■ g'(z) = {g(z)} 2 -^~ f-\ , 

az az \gJ 

in the apparently “unusual direction”. In other words, one shall search for structures in the equation 
of sums of products of the type 

f(z)g(z) + f(z)g'(z) or f(z)g{z) - f(z)g'(z), 

and then apply (31). 

In the present case, the integrating factor is e~ z 0 for all z € C, hence (29) is equivalent to 

rtf de~ z d 

(32) 0 = e-m - e-/W = ■ /(*) = K7Ml • 

Then we immediately get by indefinite integration of (32) that 

e~ z f{z ) = c, i.e. f(z) = c ■ e~ z for c € C, 
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and we have again found (30). 

This variant has the advantage that it is often possible by some small computations to solve the original 
equation, and there will usually be no problem of finding the domain of the solution, because series 
do not occur in this method. Its disadvantage is that it is not always possible to find an integrating 
factor by inspection, and furthermore, when this can be done, it requires some experience and skill. 

Second method. Determination of the germ at zq = 0. the equation is of first order, so it 
suffices to assume that /(0) = c £ C. It follows by induction from equation (29) that 

/(") ( z ) = / (n_1) (z) for all n £ N, 

and then by recursion that 

/(")(0) = /(0)=c for all n £ N. 


The formal power series is then 


-i i 

m = E -r / (n) (°) = c E -i = c • c £ c, 

n\ n\ 

n —0 n —0 

where we recognize the power series of the exponential. 

This recognition also implies that since e z is defined in all of C, the radius of convergence must be 
g = Too. 


This also follows from (6), because if c ^ 0, then 

c (n+l)l 


g = lim 

n —t-j-oo 


l n -\-1 


= lim 

n —>-+oo 


= lim (n T 1) = Too. 

n —»+oo 


Notice also that the official definition (4) gives 
1 lim„_, +00 VnJ. 


0 = 


lim 


._= lim Vn!. 

n ^ +00 t/|c| n ^+°° 


limsup n ^ +00 Vl a. 

In order to compute \fn\ we need Stirling’s formula 




meaning that 


(")' 


1 for n —> Too 


or better, the estimate 


^-0 


- exp 
e 


v 277- -|- 1 

Then in the present case, 
g = lim \/rd = lim \ 

n—> -)-oo n —^+oo L 


<n\ < 


V2™.0) 


- exp 


( — 
\ 12n 


for all n £ N, 


nn ■ — > = Too. 

e J 
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The advantage of this method is that if the initial conditions are given, then it sometimes is easy to 
find the germ j —- /( n ' ) (zq) 1 directly. The disadvantage is that the recursion formula in other cases 


may be extremely complicated and unsolvable in practice. 


Third method. The power series method. We assume that a solution of (29) is given by a power 
series 


+oo 

f( z ) = '52 a n.z n for |*| < g, 

n =0 

where the task is not finished, before we also have found g and checked that g > 0. First note that by 
the change of index n rx n + 1, followed by a corresponding change of the lower bound n = 1 to n = 0 
(always check, if the first terms in the two series are equal) 

+oo +oo 

f(z) = 'Y^,na n z n ~ 1 = ^(n+ l)a„ + i *". 

n =1 n=0 




In the past four years we have drilled 

A ^ 

81,000 km 

A 

That's more than twice around the world. 


Who are we? ^ 

O-i> 

We are the world's leading oilfield services company. Working 

O ■ H 

globally—often in remote and challenging locations—we invent, 
design, engineer, manufacture, apply, and maintain technology 
to help customers find and produce oil and gas safely. 

Who are we looking for? 

We offer countless opportunities in the following domains: 

■ Engineering, Research, and Operations 

■ Geoscience and Petrotechnical 

■ Commercial and Business 


If you are a self-motivated graduate looking for a dynamic career. 

apply to join our team. 

What will you be? 

careers.slb.com 

Sdilumberger 
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Then we get by insertion into (29) that 

+oo +oo 

0 = f(z) - f(z) = ^2 na n z n ^ 1 - ^ a n z n 

n=1 n=0 

-(-OO +00 + 00 

= ^2( n + l)a n+ lZ n - ^2 a nZ n = Y2 + l ) a n+l - An} Z n . 

n—0 n—0 n—0 

(Check that the lower bounds of the sums are identical, before we add them; if not, add or remove 
some terms to get the same lower bounds.) The unique power series expansion of 0 is then written in 
two ways, 

+OO +OO 

^2 {( n + ^“n+i -a n }z n = 0 = ^20- z n . 

71 = 0 71=0 

By the identity theorem the coefficients of the two series are equal, so by identification we get the 
following recursion formula 


(33) (n + l)a„+i — a n = 0 for all n € No (i.e. in the common range of summation). 

The easiest way to solve (33) is to multiply it by n\ ^ 0 and put b n = n\a„, because then we first get 
by a rearrangement, (n + l)!a n +i = n\a n , and then by a very simple recursion, 


b n + i = b n (= n\a n ) = i = • • • = b 0 = 0!a o = c, 

Q 

from which we get a n = — as previously. 

n\ 

Another method is the following. Assume that a Q = c. Put n = 0 into (33) to get 

c 

(0 + l)ao+i = ao = c, i.e. ai = -. 
for n = 1 we get 

c c 

(1 + l)ai+i = 2a 2 = ai = -, i.e. a 2 = ^y. 


For n — 2 we get 

c c 

(2 + l)a 2+ i = 3ci3 = 02 = 2 ! > i- e - a 3 = 

Based on these three results we assume that 

Q 

(34) a n = — for some n € No, 
n\ 

and then we shall prove that (34) also holds for the successor, i.e. when n is replaced by n+ 1, because 
then (34) by induction holds for all n £ No- 
1 and 2, so it holds indeed for some n £ No- 

(n -\- l)tt n _|_i = a n = ., i.e. Un+i 

n\ 

which is precisely (34) with n replaced by n 


Clearly, we have just proved that (34) is true for n = 0, 
However, it follows from the recursion formula (33) that 

c 

(n+ 1)!’ 

+ 1, and (34) follows by induction for all n £ N 0 - 
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The formal power series solution is then given by 


+oo .. 

f(z) = — z n = c • e 2 , for all z £ C, 

n=0 


where we have once again recognized the power series of the exponential. If we do not see this, we 
must instead apply one of the methods from the second method above to find g = +oo. 0 


The simple Example 1.4.1 above had only constant coefficients. We shall in general only consider 
linear differential equations of polynomial coefficients, in which case it is always possible to find a 
linear recursion formula by the method of power series, although this recursion formula still may be 
difficult to solve. We shall by the following two examples also demonstrate the impact of the zeros of 
ag (z) on the radius of convergence. 


Example 1 . 4.2 Solve the differential equation 

(35) (1 - z)f'(z) = f(z), 

using zo = 0 as point of expansion. 

First method. Inspection. We get by a rearrangement of (35), 

0 = (1 - z)f\z) - 1 ■ f(z ) = ^{(1 - z)f(z)}, 
hence by an indefinite integration, (1 — z)f(z) = c, for a constant c £ C, and thus 

(3 6) f(z) =—°— for :6C\{1}. 

1 — z 

Only the zero solution (for c = 0) can be extended to all of C. 

Second method, determination of the germ — j( n )(°). The expansion point is zq = 0, where 

n\ 

ao(0) = 1^0. The only zero of ao(z) = 1 — 2 is z = 1, so the power series solution is at least 
convergent in the open disc B( 0,1). When we differentiate (35) and then rearrange the result, we get 

( 1 -*)/"(*) = 2 /'(*). 

A comparison with (35) suggest that the general structure is possibly 

(37) (1 -z)f M {z) = nf ( ~ n - 1) {z) forneN. 

This is at least true for n = 1 and for n = 2. When (37) is differentiated, we get 
(1 - z )f^ +1 \z)-f^\z)=nf^(z), 
hence 

(1 ~ z)f( n+1 \z) = (n+ l)/(")(^), 

which is (37) with n replaced by n + 1, and (37) follows in general by induction. 
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Putting z = 0 into (37) it follows by recursion that 

/ (n) (0) = 0) = n(n - l)f^~ 2 \0) = ■■■ = n!/(0), 

and the Taylor series is given by 


+oo 


+oo 


m = E^r / (b) (o)*" = /(°) E ^ ^ = /(°) E 


ti=0 


n =0 


+oo 


n—0 


The radius of convergence is of course g = 1, and the sum function of the geometric series Y^ n =o ^ 


is 


1 — z 


, cf. also Section 1.3. The complete solution is 


m = 


1-z 

0 


for z G C \ {1} and c/0, 
for z € C and c = 0. 


Third method. The method of power series. Assume that the series 

+oo 

m = E a n z n of radius of convergence g > 0, 

n=0 

is a solution of (35). Then we get for \z\ < g , 

+oo +oo +oo 

0 = (1 — z)f'(z) - f(z ) = E - E na nZ n ~ E 

n= 1 7i— 1 n=0 

+oo +oo +oo 

= E^ n + 1 ) a n+lZ n - E(" + ^)0-nZ n = E(" + 1 ) {“"+1 - a n} 

71=0 71=0 71=0 

Hence, the zero function can be written in two ways as a convergent power series for \z\ < g for some 

g> o, 

+oo +oo 

E(™ +!) {°n+i -a n }z n =0=E 0 '/. 

71=0 71=0 

It follows from the identity theorem that corresponding coefficients are equal, hence we get the following 
recursion formula for n € No, 

(n + 1) {a n +i — a n } = 0, i.e. a„+i = a n , because n + 1 / 0. 

Thus by recursion, 

a n + 1 = a n = a„_i = • • • = a 0 for all n £ N 0 , 
and the series is given by 

+oo 

f(z) = aoE^- 

71=0 
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For ao = 0 we get the zero series, which is convergent for all z € C. If ao yf 0, then the radius of 
convergence can be found by either (25, 


Q 


1 


limsup n _ +00 



1 , 


or by (26), 


lim 

a n 

= lim 

a 0 

n — f— |— oo 

dn+l 

n — f+oo 

a 0 


lim 1 = 1. 

n —f+oo 


Alternatively , the result 

f(z) = —, z € C \ {1} for c £ C \ {0} constant, 

1 — z 

follows from Example 1.2.2. 0 


We still need to give an example, in which none of the standard procedures above is applicable with 
success. This is given by the following. 


hv IwMvij / i 

Maastricht University 


Join the best at 
the Maastricht University 
School of Business and 
Economics! 


Visit us and find out why we are the best! 

Master's Open Day: 22 February 2014 


Top master's programmes 


• 33 rd place Financial Times worldwide ranking: MSc 
International Business 

• I s * * place: MSc International Business 

• 1 st place: MSc Financial Economics 

• 2 nd place: MSc Management of Learning 

• 2 nd place: MSc Economics 

• 2 nd place: MSc Econometrics and Operations Research 

• 2 nd place: MSc Global Supply Chain Management and 
Change 

Sources: Keuzegids Master ranking 2013; Elsevier 'Beste Studies' ranking 2012; 
Financial Times Global Masters in Management ranking 2012 


Maastricht 
University is 
the best specialist 
university in the 
Netherlands 
(Elsevier) 


www.mastersopenday.nl 
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Example 1.4.3 Apply the three standard procedures to the inhomogeneous equation 
(38) z 2 /'(z) - f(z) = -z, 


with the expansion point zq = 0. 


The coefficients are polynomial. However, a o (zo) = ao(0) = 0, so nothing can be concluded from 
Theorem 1.4.1. 


First method. Inspection. When z ^ 0 we multiply (38) by the integrating factor , 2 exp 
Then we get, reading the equation above from the right to the left, 



^0. 


1 

— exp 

z 





1 


Clearly, the problem would be solved, if we could find an indefinite integral of — exp — in a 


neighbourhood of the expansion point Zq = 0, but this is not possible with the available methods 
known so far in this book. 


Second method. Determination of the germ 

/(0) = 0. Then by successive differentiations of (38), 



Putting z 


0 into (38) we get 


z 2 f"(z) + (2z - l)f(z) = -1, /'(0) = 1, 

z 2 f^(z) + (4z - 1 .)/"(*) + 2 f'{z) = 0, /"(0) = 2, 

z 2 / (4) (z) + (6z - l)f {3 \z) + 6/"(z) = 0, /W(0) = 12 


It is left to the reader to prove by induction that in general, 

z 2 / {n+1) (z) + (2 nz - 1 )/ (n) (z) + n(n - l)/ (n " 1} (z) = 0 for all n > 1. 


Hence, for z = 0, 

/ (n) (0) = n{n - l)/ (n-1) (0), for n > 2. 

We divide this equation by n\{n — 1)!, and then we get by a simple recursion, 

/(”) (0) _ / (ra_1) (0) _ _/®(0)_2_ 

n\(n — 1)! (n — l)!(n — 2)! 2!1! 2 

so the Taylor coefficients are 

1/W(0) = (n-1)! forneN. 

Then the formal Taylor series becomes 

+oo 1 +oo 

Y^^f^(0)z n = J2(n-iy.z n . 

n =0 n =1 
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Its radius of convergence is computed by (6), 

(n — 1)! 1 

g = lim -:— = lim — = 0, 

n —►-)-oo Til n —»+oo fl 


so the Taylor series is only convergent for z = 0, and this formal 


solution is useless in the applications. 


Third method. The power series method. Assume that (38) has the convergent power series 
solution 

+oo 

f(z) = '52 a nZ n for \z\ < g. 

n=0 

It follows by insertion into (38) that 

+oo +oo +oo +oo 

z 2 f{z) - f(z) = z 2 ^2 na n z n ~ 1 - ^ a n z n = ^ na n z n+1 - ^ a n z n 

n=l n =0 nOl n =0 

+oo —|—oo +oo 

= ^(n - l)a n _i z n - '^2 a nZ n = —a 0 + ^2 {( n _ l)an—l - a n } z n . 

n=2 n=0 n= 1 

This expression is equal to —z, if —ao = 0 and 0 • ao — a\ = —1, i.e. a\ = 1, and in general, 
a n = (n — l)a„_i, for n> 2. 

Then we get by recursion (the details are left to the reader), 
a n = (n.l)lai = (n - 1)!, 
so the formal series solution is 

+oo 

y^(ra — 1 )!z”, unfortunately with g = 0. 

71=1 

Thus, this series is divergent, whenever z ^ 0, and the method is not applicable. 0 


Example 1.4.3 shows that if ao (z o) = 0, then Theorem 1.4.1 does not apply, and the problem may not 
be solvable by any of the three suggested standard procedures. We shall later on in connection with 
Laurent series prove that we may in some cases be able to solve such a linear differential equation, 
even if ao (zo) = 0> while in other cases this is not possible, because the equation in reality should be 
solved on a so-called Riemann surface cf. Ventus: Complex Functions Theory a-3. 


1.5 Zeros of analytical functions 

We have already in Section 1.4 seen that the zeros of an analytic function may have some influence on 
the behaviour of the function. We shall in this section see by using some abstract topological results 
shown in Ventus: Complex Functions Theory a-1 that the zeros of an analytic function also have other 
unexpected consequences, which are stronger than possibly similar results for real C°°(K) functions. 
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Theorem 1.5.1 Assume that f : Cl —> C is analytic and not the zero function. For every given point 
z £ Cl there exists an n £ No, such that f^ n \z) 7 ^ 0. 


We here adopt the convenient notation, f^(z) := f(z ) (no derivatives of f{z)). 

Proof. Given an open domain Cl in C, and an analytic function / : Cl —> C. We define a subset 
EC by 

(39) E:={z£Cl | / (n) (-) = 0 for all n £ N 0 } . 

This means that E is the set of points z £ Cl, for which both f(z ) = 0 and all its derivatives /= 0. 

If z 0 € E , then clearly the germ of the Taylor series is just the zero sequence. According to Theo¬ 
rem 1.3.1, the Taylor series has its sum function f(z) in the largest open disc B (zg, g) C Cl, hence 
f(z) = 0 for every z £ B (zg, g), so by the definition (39) of the set E we have proved that B (zg, g) C E, 
proving that E is an open subset of Cl, cf. Definition 2.1.1 in Ventus: Complex Functions Theory a-1. 

On the other hand, the complementary set 

n\E={J {* e n |/<">(*)^o}= (J {/<">} “ (c\{o» 

n —0 7i—0 

is also open, because C \ {0} is open and every derivative /^ is continuous, cf. Definition 2.1.2 in 
Ventus: Complex Functions Theory a-1 , so each set {/(”)} (C \ {0}) is open, and every union of 
open sets is again open. 

Then the open domain Cl = E U {Cl \ E} is written as a disjoint union of two open sets. Since 
every domain by definition is connected, it follows from Corollary 2.1.1 in Ventus: Complex Functions 
Theory a-1 that either E = Cl or E = 0. 

If E = Cl, then / is identically zero, which was excluded in the assumptions. Hence E = 0, so for 
every given z £ Cl there exists an n £ No, such that f^{z) 7 ^ 0, and the theorem is proved. □ 

Assume that f(z) is analytic and not identically zero. Let zq £ Cl be a zero, i.e. / (zo) = 0. It follows 
from Theorem 1.5.1 that there is at least one n £ N, such that /^ (zg) 7 ^ 0. 


Definition 1.5.1 Let Zo £ Cl be a zero of the analytic function f : Cl —> C, where / 7 ^ 0. We say that 
the zero zg has the order, or multiplicity n £ N, if n is the smallest integer for which f ^ (zg) 7 ^ 0. 


In order to motivate this definition we consider the power series expansion of / with the zero z g of 
order n as expansion point. We get in a neighbourhood of zg that 

f(z) = - ] f W {z 0 )-{z-Zg) n + r ^-f^ 1 Hzo)-{z-Zg) n+1 + --- 

n\ (n+1)! 

= (z- Zg) n ■ | ^ f (n) {Zg) + / ( ” +1) (*o) • (z - Zg) + • • • | = (z - Zg) n ■ g{z), 

where g{z) is analytic in the same neighbourhood of zg as f(z), and where furthermore, g (zg) 7 ^ 0 . 
This means that z = Zg is precisely n times a zero of /, explaining the notation. Furthermore, notice 
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that since g{z) is continuous, we can find a (possibly smaller) open neighbourhood ui of Zo , such that 
g(z) 0 for all z £ ui. Then clearly also f(z) ^ 0 for every z £ ui\ {-Zo}, and we have proved 


Theorem 1.5.2 Let f : Cl —> C, / 0, be analytic. Every zero Zo of f is an isolated point, which 

means that there is an open neighbourhood of zq, such that Zq is the only zero in to. 


Theorem 1.5.2 immediately implies the following stronger version of Theorem 1.2.6, The Identity 
Theorem. 


Theorem 1.5.3 The Identity Theorem. Let f : LI —> C and g : LI —> C be analytic functions in the 
same domain Cl. If the set {z £ Cl \ f(z) = g(z)} has an accumulation point lying in the set Cl . then 
f and g are identical. 



> Apply now 


REDEFINE YOUR FUTURE 

AXA GLOBAL GRADUATE 
PROGRAM 2015 


redefining /standards 
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Remark 1.5.1 The underlined assumption that the accumulation point lies in the set fl is very 
important. We shall in Example 1.5.1 show that what without this assumption, we cannot conclude 
that the two functions are equal. 0 

PROOF. First recall that z 0 £ LI is an accumulation point of a set A C LI, if for every r > 0, 

{z £ A | 0 < \z — zo\ < r} / 0. 

Notice that the inequality 0 < \z — zq\ excludes zq from this set. 

Assume that zq £ LI is an accumulation point of the set 
E :={z£Ll | f(z) = g(z)}. 

Choosing r = — . n £ N, and then z n £ E, such that / (z n ) = g (z n ), we define a sequence (z n ), for 
n 

which z n —> zo for n —> +oo. Now, / — g is continuous, and / (z n ) — g (z n ) = 0 for all n £ N. Since 
Zq £ LI, it follows from the continuity that also / ( 20 ) — 3 { z o) = 0, so zq £ E. 

It follows from the definition of the sequence (z n ) that Zq is not an isolated zero of the analytic 
function f — g, so we conclude from Theorem 1.5.2 that / — g is the zero function, and thus g = g by 
a rearrangement. □ 


Example 1.5.1 Consider the zero function 0 and the function sin - for z £ Q = C\ {0}. Clearly, the 

z 

two functions are different from each other. On the other hand, both functions are zero on the set 


E := 



n £ Z \ {0} 


which has zq 
0 


0 ^ as an accumulation point, because-> 0 for n —> + 00 , and also for n —> — cxd. 

titt 


Another unexpected consequence is the following theorem. 


Theorem 1.5.4 Let ip : I —> C be a function defined either on a real interval I, or on a differential 
curve I in C. If H is an open domain in C, which contains I, then there is at most one analytic 
function f : fl —> C, such that 


f{x) = ip{x) 


for x £ I, an interval 


f(z) = <p(z) 


for 


z £ /, a piecewise C 1 curve in C. 


Proof. In both cases all points of I are trivially accumulation points of /, so the theorem follows 
immediately from Theorem 1.5.3. □ 

The reader should be surprised that the values of an analytic function on a one-dimensional curve 
uniquely determines / : Cl —> C in its domain. This means that if we change the analytic function 
/ on e.g. a real interval to another analytic function g , then we change without any time delay / 
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to g over all of 12. This implies that analytic functions are not suited directly to describe causality , 
because e.g. the impact of a sudden force applied to a system will evolve in time through the system 
and not immediately. Fortunately, it is possible to describe indirectly the causality (approximately) 
by analytic functions. However, this is not the right place to go further into this discussion. 

One should also note that Theorem 1.5.4 states that there is at most one analytic function / : 12 —> C, 
the restriction of which to an interval is a given function tp : I —> C. It is not hard to construct a 
function in C 00 (M) which cannot be extended to any analytic function, not even locally! 


1.6 Simple Fourier series 

It follows from the definition e* e = cos© + i sin© for 0 £ K that e* e has the period 2i r, so e m0 has 

2n 

for fixed neM the period — in 0. 

n 

The classical Fourier Series Theory in real calculus states that every piecewise C 1 ([0,27r]) function 
can be represented in the complex form 

+oo 

(40) *>(0) ~ Y, c " ein0 ’ 


where 


r 2n 


C-n • — 


2n 


<p(Q) e 


—in@ 


d© 


for n £ Z. 


Remark 1.6.1 The symbol ~ in (40) indicates that <p(Q) is equal to its Fourier series in the sense 
of L 2 , i.e. 


lim 

N —>-+oo 



N 

¥>(©) - Y c " e 

n=—N 


2 


d0 = 0, 


which can also be interpreted as convergence in energy. It can be proved that if / Q 27r |cp| 2 cl0 < +oo, 
or just / Q 27r |<p(0)| p d0 < -|-oo for some p > 1, then (40) holds with pointwise equality sign for almost 
every 0 £ [0,27r]. (For p = 1 this statement is wrong as proved by Kolmogorov in the early 1920s.) 
The proof of the statement above is extremely difficult, and the result is of limited value, because one 
usually cannot specify for which 0 (40) holds with equality sign, even if we know that it holds “for 
almost every 0. 0 


Consider an analytic function / : 12 —> C, where 12 is an open neighbourhood of 0. Then we can 
construct the convergent Taylor series 

+oo 

(41) f(z) = Y a n zU for M < 6- 

71=0 
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Using polar coordinates z = re 10 for fixed r £ ]0, £?[, we obtain the Fourier series of a function ip r , 
given by 

+ OO 

(42) ^(0) := / (re ie ) = ]T {a n r n } e inQ , 

n= 0 

so when we compare with (40) we see that we have the Fourier coefficients 

c n = a n r n for n £ N 0 and c n = 0 for n £ Z_. 

In (42) the sum is only over n £ No, while we in (40) use all n £ Z in the sum. This suggests that 
at least concerning Fourier series, it would be quite natural also to allow power series (41), where 
negative exponents occur. Such series are indeed very useful in the applications. However, they 
cannot be defined everywhere in a disc B(0,g), because they are at least divergent for z = 0 , where 
lim^o z ~ n = oo. 

We shall later in Chapter 3 study such power series of negative exponents. Such series are called 
Laurent series. The motivation for their introduction is here given by the fact that we are missing 
some very natural terms in the Fourier series (42), but it will turn up that there is far more in these 
Laurent series than one would expect at a first glance. Their main applications are in the so-called 
residue calculus, which is a powerful device to compute many definite integrals and infinite sums, 
including some which cannot be computed by methods from the real calculus. The Laurent series 
with no positive exponent are furthermore used in the theory of the 3 transform,, which is a discrete 
form of the Laplace transform. 


Example 1.6.1 Let us play a little with this connection between Complex Functions Theory and 
Fourier series. We know that 

e z = e x cos y + ie x sin y, z = x + iy £ C. 

Using polar coordinates, 

z = r e 1 i.e. x = r-cosO and y = r- sinO, 
we get by insertion, 

e 2 = e r(cos0+i sine) = co S e| c0s ( r gin Q) + { sin ( r gin©)}, 
and 


+00 +00 +00 +00 

= V — z n = V — r n e ine = V — cosn© + i V — sinn©. 
' n\ ' n) T)\ n) 


n\ 


n\ 


n —0 n—0 n =0 

By identifying the real and the imaginary parts, 

+OO 


n\ 


jr cos 0 


COS 


n —0 


(r sin©) = y — cos nO, 


+oo 


0 r cos © 


sin(r sin©) ) —- sin?r©. 

^ n\ 


Notice that the summation starts at n = 1 in the latter sum, because sin 0 • © = 0 for n = 0, no matter 

0. o 
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Example 1.6.1 is only a demonstration of what we already can obtain. The following results are more 
important, so for this reason they are not relegated to an example. 

This time we use that 


1 


+OO +OO 


+OO 


+ OO 


— = = ^Ve ine = ^r n cosnQ + i ^Vsinn© for r = \z\ < 1, 


1 — z 

n =0 n =0 

and that also 

1 1 — x 


n =0 


n=l 


(43) 


1 — z (1 — x) 2 + y 2 (1 — x) 2 + y 2 
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Using polar coordinates the dominator becomes 


(1 — x) 2 + y 2 = 1 — 2x + x 2 + y 2 = 1 + r 2 — 2 r cos 0 , 


so 


1 


1 — r cos 0 


+ i ■ 


sin© 


+oo 


+oo 


— = r n cos n® + i r n sin n®, 0 < r < 1. 


1 — z 1 + r 2 — 2 r cos 0 1 + r 2 — 2 r cos 0 

n =0 n =1 

When we identify the real and the imaginary parts we obtain the following convergent Fourier series, 


(44) 


1 — r cos 0 

_ y'N+OO 

A-/71=0 

1 + r 2 — 2 r cos 0 

r sin 0 

Y^+OO 
_ 2^71=1 

1 + r 2 — 2 r cos 0 


for 0 < r < 1 and 0 £ 


We shall later also need the Fourier expansion of 
l + z 2 


- 1 


for hi < 1. 


1-z 1-z 

It follows immediately from (44) that 


l + z 
l-z 


= 2 


— 1 + i 


2 r cos 0 


(45) 


1 — r cos 0 
1 + r 2 — 2 r cos 0 

1 — r 2 

-H - 'f' - 

1 + r 2 — 2 r cos 0 1 + r 2 — 2 r cos © 

+oo +oo 


1 + r 2 — 2 r cos 0 
2 r sin 0 


= 1 + 2 r n cos n 0 + 2 i r" sin n© for 0 < r < 1 and 0 £ 


We shall finally prove a special case of Parseval’s equation. This result will be applied later. 


Theorem 1.6.1 Parseval’s equation. Let 

+oo +oo 

f ( z ) = ^2 anzn and 3(~) = bz 71 


n —0 


71—0 


be two analytic functions defined in 1?(0, g). Then for every r £ [0, g[, 

27T _ +00 

(46) — / / (re l& ) g (re ie ) d® = ^a n b n r 2n . 

^ J 0 n = 0 


Proof. Since both series are absolutely convergent for r = \z\ < g, we may apply termwise multipli¬ 
cation, from which we get 


+oo +oo 


(47) / (re* e ) g (re i0 ) = ^ 


n+m i(n-m)Q 
, 1 C ? 


771=0 71=0 
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which for fixed r is absolutely and uniformly convergent, because )T) m |a n & TO | r n+m is a converging 
majoring series. Hence, we may integrate (47) termwise. Then notice that 

2tt 

e i(n-m)0 d Q _ q, 

so the only relevant terms of (47) by this integration, are the terms given by n = to, and we get 
trivially, 


1 

27T 



/ (re 10 ) g ( re ie ) d 0 


1 

27T 


+oo 

y, a n b n r 2n 


n —0 



+oo 

yOnbnr 2 ". 


n=0 


□ 


Corollary 1.6.1 Lef f(z) = Yjn=o a n^ n f° r kl < 6- Then 


c2tT 


+oo 


(48) ^ / |/(^ e )| = E 


aJ 2 r 2n . 


n—0 


Proof. Just put g = f and b n = a n into (46). □ 


1.7 The maximum principle 

Another strange property of a non-constant analytic function / : fl —> C is that the continuous 
function |/(z)| can never have a local maximum at an interior point Zq G fi. We have more precisely 


Theorem 1.7.1 The maximum principle Let / : 12 —> C be analytic in an open domain 12. If \f(z)\ 
has a local maximum at a point Zo G 12, then f(z) is constant in 12. 

PROOF. Assume that \f(z)\ has a local maximum at an interior point zq G f l. Then we have a 
convergent series expansion 

+oo 

(49) f(z) = y a n (z - Z 0 ) n for \z - z 0 \ < g. 

71—0 

Choose r G ]0, p[ and put z = Zq + re 10 . It follows from Parseval’s eguation (48) that 

+oo 

/ (z 0 + re l& )\ 2 d0 = y \a n \ 2 r 2n . 

71—0 

By assumption, |/(^o)l is a local maximum, so due to the continuity we can choose r G]0, p[, such 
that \f(z)\ < \f (zq)\ in the closed disc B [z 0 ,r], and we get the estimate 

+°° -| -j p2n 

/ \f (z 0 + re ie )\ 2 d 0 < — / \f(z 0 )\ 2 d 0 = |/( z 0 )| 2 = |a 0 | 2 , 

n —0 ^ JO 71 Jo 
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hence, 

+ OO 

y. \a n \ 2 r 2n < \a 0 \ 2 . 

n— 0 

Since we have chosen r > 0, this is only possible, when all a n = 0 for n £ N, and we get by insertion 
into (49) that f(z) = ao in Ll. □ 

There is of course a corresponding minimum principle. This is, however, slightly more complicated, 
because every zero of f(z) for natural reasons is a minimum point of the function \f(z)\. The simplest 
way to formulate it is 


Theorem 1.7.2 The minimum principle. Assume that the analytic function f : Ll — > C is not 
constant in the open domain f l. If \f(z)\ has a local minimum at an interior point zo € f l, then 
f (zq) = 0 , i.e. z o must necessarily be a zero of f{z). 

PROOF. Assume that / (zq) y^ 0. Then there is an open subdomain wCf2, such that f(z) y^ 0 for all 

z £ to, where also zq £ ui. The function g(z) = ——— must then be analytic in u>, so when \f(z)\ > 0 

/ 0 ) 

has a local minimum at zq £ u>, then \g(z)\ has a local maximum at the same point zq £ to. By the 
maximum principle this is not possible, so we conclude that / (zq) = 0 . □ 


Theorem 1.7.3 Let f : LI —> C be analytic in a bounded open domain, and assume that f : LI —> C 
is continuous on the closure f l of LI. Then |/(z)| has its maximum lying on the boundary. 


PROOF. It follows from the assumptions that LI is compact. The function |/| : Ll —> C is continuous, 
so it has a maximum on Ll. (Main theorem for continuous real functions.) By the maximum principle, 
the maximum cannot be attained at an interior point, unless the function is constant, so the maximum 
is attained at a boundary point, no matter if |/| is constant or not. □ 

The conclusion of Theorem 1.7.3 is not true for unbounded closed domains. A simple counterexample 
is 


f(z) = e z for 2 £ Ll := {z £ C | IRz > 0}. 

Clearly, |/(z)| = e z is unbounded in Ll, so no maximum exists. On the boundary, however, i.e. on the 
imaginary axis, we have \f(iy)\ = |e ly | = 1 , which is bounded. 

The following is a partial inverse. 


Theorem 1.7.4 Phragmen-Lindelof’s theorem. Assume that the function f(z) is analytic in the half 
plane IRz > 0, and assume that f(z) is bounded and continuous on the boundary (the imaginary axis), 
i.e. \f(iy)\ < M for some constant M > 0 and all real y £ R. 

Furthermore, assume that there exist real constants a < 1 and K > 0, such that 

(50) \f(z) | < K ■ exp (r a ), for all z = r e l& for which 3 ftz > 0. 

Then \f(z)\ < M everywhere in the half plane IRz > 0. 
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PROOF. Choose any b S]a,l[ and any e > 0. When we consider the function g(z) = f(z) exp (— ez b ), 
we get the estimate 


IsC 


01 = l/( 2 )|exp(-er b cos( 60 )) < |/( 2 )|exp ^-£r 6 cos 


(51) < K ■ exp |^r a — e cos 

Let zq be any point in the half plane and choose R, such that |zo| < R, and such that \g(z)\ < M on 
the semicircle \z\ = R, 3?;: > 0. 
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This is possible, due to the choices of the constants e and b , and due to the estimate (51). 

The closed half disc is compact, so it follows from Theorem 1.7.3 that |g ( 20 ) | < M , hence 

|/(zo)| = Iff ( 20 ) I • |exp (s Zq) I < M ■ exp (e • r 0 cos(60 o )). 

This holds for every e > 0, hence also by the limit process e —> 0+, by which we get |/ ( 20 ) | < M. 
Since zq was chosen arbitrarily in the half plane, the theorem is proved. □ 

Notice that (50) is not fulfilled for any a < 1 for the function f(z) = e z , so the theorem cannot be 
extended to a = 1. 

There exists, however, another version of Phragmen-Lindelof’s theorem, in which we allow a > 1, 
though the conclusion of course becomes weaker. 


Theorem 1.7.5 Weak Phragmen-Lindelof’s theorem. Let f(z) be analytic in a vertical strip 

X\ < Sftz < X 2 (which means more precisely that f(z) is analytic in some open domain containing this 

strip). We assume that \f(z)\ < 1 on the boundary of this strip. 

If there exist constants a > 0 and K > 0, such that 

(52) |/( 2 )| < A'• exp (|z| a ), when Xi < $iz < x 2 , 

then |/( 2 )| < 1 everywhere in the vertical strip. 


PROOF. Let A > a. Since x = 5?2 is bounded in the strip, there exists an y 0 , such that 
(53) \f(z)\ < exp (|y| A ) , when \y\ > y 0 and x S [x 1 ,x 2 ] ■ 

Choose p € N, such that to = 2 + 4p > A. 



Figure 4: Proof of variant of Phragmen-Lindelof’s theorem. 


If 2 = re ,e in the strip is large in the sense that |Cs 2 1 = \y\ > yo, then 
2 m = r m (cos to© + i sin to0), 
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where m0 lies close to 7 r. More precisely, there is a q £ Z, such that |m0 — (7r + 2qn)\ is sufficiently 
small. 


Choose any e > 0 and consider the function 


(54) g e {z) = 7 ( 2 ) exp (e 2 ™). 


If 2 lies in the strip x\ < JJ 2 < X 2 , then we get the estimate 
\g e {z)\ < exp ( \y\ x ) exp (er™ cos m0). 

This shows that the function g e (z) is bounded for large yo = yo(e), and |y e (^r)| < 1, when |?/| = yo in 
the strip. The latter follows from the choice of to, by which we get cos to.0 < — a < 0. 

Since \g e {z)\ < exp (ex™) on the edges of the rectangle, we conclude that \g e {z)\ < exp (ex™) for all 
2 inside the rectangle. Finally, since we already have proved that |g e ( 2 )| < 1, when 2 in the strip 
satisfies 1321 > yo, we conclude that 

\de{z)\ < exp (ex™) for all e > 0 and X\ < 3?2 < X 2 . 

Hence, it follows from (54) for every fixed 2 in the strip that 

|/( 2 )| < exp (ex™) • exp (e| 2 |™) for all e > 0. 

By taking the limit e —+ 0+ for every fixed 2 in the strip it follows that \f(z)\ < 1. □ 

The following result is often useful in the most unexpected situations. 

Theorem 1.7.6 Schwarz’s lemma. Let f : B(a,R) —> C be analytic, where f(a ) = 0 and |/( 2 )| < M 
for every z € B(a,R). Then 


(55) \f(z)\ < — \z — a\ for all z £ B{a,R). 

R 


If we have equality in (55) at some point 2 £ B(a, R), z yf a, then f(z) has the structure 


(56) f{z) = e iQ • — • (2 — a), for all z £ B{a , R). 

R 


Proof. It follows from the assumption /(a) = 0 that 



f(z) = ^2 an ( z ~ a )" = {z - a) ^2 a n+i( z ~ a ) n for z e B ( a > R )- 


n =1 n=0 


If we put g(z) = Yln=o a «+i ( z ~ «)"> then 



2 — a 



for 2 = a, 


and g : B(a, R) —> C is analytic. 
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Let r £ ]0, R[ and \z — a\ = r. Then we get the estimate 

, , M 
\g{z)\ < 

hence by the maximum principle, 

M 

\g(z)\< — for all z £ B(a,r). 

Since r < R can be chosen arbitrarily close to R , we conclude by taking the limit r —> R— that 
M 

\g{z)\ < — for all z £ B{a,R), 

K 

hence 

M 

\f(z)\<-.\ z -a\. 


Then assume that we have equality at a point zo £ B(a, R), zo ^ a. Then |g(z)| has its maximum at 
the point zq in the interior of B{a,R ), so g{z) is constant by the maximum principle, hence 

g{z) = e 10 • and thus f(z) = e ie -^-(z-a). □ 



Figure 5: Proof of Corollary 1.7.1. 


One nice application of Schwarz’s lemma is the following 


Corollary 1.7.1 Let f l be an open domain with a £ LI. If there exists an analytic map f , which maps 
Lt bijectively onto 13(0,1), such that f(a) = 0, then every bijective analytic map g : Ll —» 13(0,1), for 
which g(a) = 0, is given by 

(57) g{z) = e l ®f(z) for some 0 £ [0,27r[. 

Thus, these maps are uniquely determined apart from a rotation of the disc 13(0,1). 
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Proof. Clearly, every map of the form (57) satisfies the given conditions. 

Let g : —> B{ 0,1) be any bijective analytic function for which g{a ) = 0. 

The composite map go/ -1 : B{ 0,1) —» B{ 0,1) is bijective, and go/ _1 (0) = g(a) = 0. Using Schwarz’s 
lemma it follows from | g o f~ 1 (z) | < |;?|. 

Now, the inverse of g o / _1 , i.e. / o g -1 , fulfils precisely the same properties, so we also have 
\f ° | < \w\. 

Then, by putting w = g o 

M = |/ 0 5 _1 M| < \w\ = \g ° / _1 (^) | < \z\, 

so we must have equality. Then by Schwarz’s lemma , 

e l& ■ (/ ° 5 _1 ) (w) = w. 

Finally, putting w = g(Q , we get g(£) = and the corollary is proved. □ 
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2 Harmonic Functions 


2.1 Harmonic functions 


The harmonic functions are solutions of the Laplace equation , 


(58) 


Au := 


d 2 u 
dx 2 


d 2 u 

dy 2 


= 0 , 


which has many applications in the plane i the applied sciences. We shall see that the harmonic 
functions are also closely connected with the analytic functions, so many results on harmonic functions 
can easily be derived from the theory of analytic functions. 


Let / : fi —> C be an analytic function in the open domain f 1. Then we have proved in Ventus: 
Complex Functions Theory a-1 that / is of class C°°(Ll). If we split / into its real and imaginary 
parts, f = u + iv, then we get by Cauchy-Riemann’s equations 

. du dv dv du 

l{z) = di + 'Si = ay-'W 

. d?v_ _ d 2 v _ . d 2 u 
dx 2 % dx 2 dydx * dydx 

etc., from which we conclude that 3?/ = u{x,y) and Lsf = v(x,y) are both of class C°° in the real 
variables ( x,y ) £ f l. In particular, we may interchange the order of differentiation. Since f(z) is 
analytic we can apply Cauchy-Riemann’s equations, 

du dv 

dy dx ’ 

d 2 u d v d 2 v 

dy 2 dydx dxdy ’ 

so by adding these two expressions we get the Laplace equation (58). 

Similarly, we prove that Au = 0. 

We introduce formally, 


du dv 

= W and 
dx dy 

hence by differentiation, 

d 2 u d 2 v 
dx 2 dxdy 


and 


Definition 2.1.1 Let Ll C R 2 be an open domain. A real function u £ C 2 (Ll) is called harmonic in 
f }, if Au = 0 in fi. 

It follows from the above that 

Theorem 2.1.1 If f(z) = u + iv is analytic in an open domain LI C C, then its real and imaginary 
parts are both harmonic functions in Ll Cl 2 . Here the planar domain Ll is considered both as a subset 
of C and of Si 2 , depending on the context. 

The harmonic functions are of course important in every two dimensional potential theory. We list 
some elementary harmonic functions in Table 1, page 66. 
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Definition 2.1.2 Assume that u and v are harmonic in the open domain Cl. If the pair ( u , v) satisfies 
Cauchy-Riemann’s equations, then we say that v (the second coordinate) is the harmonic conjugated 
of u (the first coordinate), or that the pair (u,v) is an harmonic conjugated pair. 


We notice that if ( u, v) is a harmonic conjugated pair, then {—v, u ) is also a harmonic conjugated pair, 
while (v, u) is only harmonic conjugated, when both u and v are constant functions. Thus, harmonic 
conjugation is not a symmetric relation. 

ft follows immediately from Cauchy-Riemann’s equations that 


Theorem 2.1.2 If (u,v) is an harmonic conjugated pair in Cl, then f = u + iv is an analytic function 
in Cl. 


Assume that u(x,y) is harmonic in an open and simply connected domain Cl, i.e. 

— = forced 

dx 2 dy 2 tor fz, y) £ U. 

ft follows that the differential form 

du , du , 

- — dx + — dy 
ay ox 

is closed, so it is exact, Cl being simply connected. 

Consider a fixed zq £ Cl. ft follows from the above that the function 


(59) v(x,y) := f 

J z c i 


du du . 
-dv iX + dx dV >' 


z = x + iy £ Cl, 


is uniquely defined in Cl. 


du du 

Since u is of class C 2 (Cl ), it follows that — and — are of class C 1 (d), so we conclude that v £ C 2 (Cl). 

dx ay 

In particular, 

. . dv dv du du 

(60) dv = —dx+ — dy = — dx + — dx, 
dx dy dy dx 

hence Cauchy-Riemann’s equations follow by identification. Since u, v £ C 2 (Cl), we conclude that 
/(z) = u + iv is analytic, and we have proved 


Theorem 2.1.3 If u(x,y) is harmonic in a simply connected open domain Cl, then there exists an 
analytic function / : fi —> C, such that u(x,y) = Sft/(z) in Cl. In particular, u is of class C°°(Cl), and 
(u,v) is a harmonic conjugated pair, when the function v{x,y) is given by (59). 

Example 2.1.1 If Cl is an open domain, which is not simply connected, and u(x,y) is harmonic 
in Cl, then there does not necessarily exist an analytic function / defined in all of Cl, such that 
5R/(z) = u(x,y). A very important counterexample is 

u(x, y) = - In ( x 2 + y 2 ) = In |z| = lnr, 
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which is harmonic in the (not simply connected) domain C\ {0}. By considering the analytic function 
Log z in the sliced and connected domain C \ {R_ U {0}} we see that 

3? Log z = ^ In (x 2 + y 2 ) for z £ C \ {R_ U {0}} , 

where Log z cannot be extended analytically to all of C \ {0}. 0 

Assume that u, Vi and n 2 are harmonic functions in the same domain Cl, and that both (u,v i) and 
(u,v 2 ) are harmonic conjugated pairs. Then fi(z) := u + iv 1 and / 2 (z) := u + iv-2 are both analytic 
in Cl, so 

f(z) = fi{z) - h(z) =i{v 1 - v 2 } 
is analytic in Cl. Then by Cauchy-Riemann’s equations, 
d , d , 

Si {V1 - n] = Sy {V1 - V2} = O’ 

hence v\ = V 2 + C for some constant C £ M, and we have proved 


Corollary 2.1.1 If u(x,y) is harmonic in a simply connected open domain Cl, and zo £ Cl is a fixed 
point, then all harmonic conjugated functions of u are given by 


(61) v(x,y).= J^-^dx+^dy} + C, 


C £ K arbitrary. 


If v(x, y) is harmonic in a simply connected domain Cl, it follows similarly that all harmonic functions 
u, for which (u,v) is a harmonic conjugated pair, are given by 

r f dv dv j 

(62) u{x,y) := / < — dx — —— dy > + C, C £ R arbitrary, 

J z 0 l Qy J 

(notice the change of sign compared with (61)), and f(z) = u + iv is analytic in Cl. 


Example 2.1.2 We shall show that the function u(x,y) = x 3 — 3 xy 2 is harmonic, and find all its 
harmonic conjugated functions. Clearly, u £ C°° (R 2 ), and 

du 2 o 2 d 2 n r. , d u r. d 2u 

— =3x -3 y , = ox and — = -6 xy, —^ = -Ox, 

ox ox z ay oy z 

thus 

. d 2 u d 2 u „ 

Ait = 7 — + 7—7 = 6a; — 6x = 0, 
ox z oy z 

and u{x, y) is harmonic. 

du du 

Then insert the expressions of 7 — and 7 — above into (61). Choosing zq = 0 we get 

ox dy 


r* n x ,y> 

v(x,y) = / { — {fsxy) dx + (3a : 2 — 3y 2 ) dy} + C = / { (6xy dx + 3x 2 dy) — 3y 2 dy} + C 

Jz 0 - 1 ( 0 , 0 ) 

r{ x iv) 

= d (3 x 2 y — y 3 ) + C = 3x 2 y — y 3 + C, C £ R arbitrary. 

J< 0,0) 
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Here we have applied that the differential form 6xydx + (3a ; 2 — 3a/ 2 ) dy is exact, so it can be written 
in the form dv. We have in particular applied the rule of the differential of products, 

d(ipijj) = ip dip + ip dip 

in the reverse direction, i.e. from the right to the left. 

Alternatively we may integrate along a simple curve composed of axiparallel segments. 0 


Remark 2.1.1 The beginner often makes the error that f 6xy dx and J {3.t 2 — 3 y 2 } dy are computed 
separately, where the other variable erroneously is considered as a constant. By an addition we get 
the wrong result 6 a Ay — y 3 , which can easily be checked by Cauchy-Riemann’s equations, which are 
not fulfilled in this case. 0 


It is in some cases possible from a given harmonic function it in a simply connected domain fl directly 
to find the corresponding analytic function f{z), such that 

u(x,y) = 


First, it follows from Theorem 2.1.3 that 


f(z) exists. Then by Cauchy-Riemann’s equations, 


(63) m = g 


dv du . du 
dx dx dy 


du 1 du \ 
dx i dy) 
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Since the function u(x,y) is given, (63) shows that f'(z) can be computed directly as an analytic 
function, so the right hand side of (63) must be expressed as a function in z (= x + iy) alone. Finally, 
we get f[z ) by an indefinite integration, where the arbitrary constant is determined, such that ^Rf(z) = 
u(x, y) in fi. 


Example 2.1.3 We shall reconsider the function u(x,y) = x 3 — 3 xy 2 from Example 2.1.2. We can 
now prove that u(x, y) is harmonic without inserting into the Laplace equation A u = 0. In fact, if u 
were harmonic, then there would exist an analytic function f(z), such that 

f'(z ) = ^ ^ = 3x 2 - 3y 2 - i{—6xy} = 3z 2 . 

ox dy 

Since f'(z) = 3 z 2 has the indefinite integral z 3 + C , it is straightforward to check that 
5iz 3 = a ; 3 — 3 xy 2 = u(x, y), 
so u is an harmonic function. 


In the present case one cannot claim that Example 2.1.3 is easier than Example 2.1.2. However, it is 
indeed very easy to provide examples, where 


d 2 u 
dx 2 


and 


d 2 u 

dy 2 


become very cumbersome to compute, while 


du 

dx 


and 


du 

dy 


occurring in f'(z ) given by (63) are easy to find, so it becomes a simple task to find f'(z ) as a function 
of 2 alone. 0 


2.2 The maximum principle for harmonic functions. 

Given an harmonic function u(x, y) in an open domain O, and assume that u has a local maximum (or 
a local minimum) at an inner point zo ~ (xq, yo ) £ There exists an r > 0, such that B ( zq , r) C f l. 
Since B (zo,r) is simply connected, we can find an analytic function f(z), such that 3?/(z) = u(x,y) 
locally in B (z 0 , r). 


Since 



= e 3 ^ 2 ) = e{x,y) > 0, 


we conclude that |e-^ 2 ) | has a maximum (a minimum) at zq. The function e-h z ) is analytic, hence e-h z ) 
is a constant by the maximum (the minimum) principle, and u(x,y) is constant in B (zo,r). Since u 
is continuous, it is constant in every simply connected subdomain u> C P, hence also in f 1 itself, and 
it follows by contraposition that we have proved 


Theorem 2.2.1 The maximum (minimum) principle for harmonic functions. Letu(x,y) be harmonic 
and not a constant in the open domain fb Then u(x,y) has neither a maximum nor a minimum in 


fi. 
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Let K C Cl be a compact set, and let u{x,y) be harmonic in f l. Since u(x,y) only has real values, 
it follows by the Main Theorem for the Continuous Functions that u{x, y ) has a maximum and a 
minimum in I\. It follows from the above that these values are always obtained on the boundary of 
K , no matter if u(x,y) is a constant or not in f 1. It therefore immediately follows that we have 


Corollary 2.2.1 Assume that C1 is an open and bounded domain. Let u{x,y) be continuous on the 
closure f l and harmonic in Cl itself. Then u(x,y) has its maximum and minimum on the boundary 

dci = ci\ci of Cl. 


Proof. Just notice that Cl is compact and apply the previous argument. □ 

In Ventus: Complex Functions Theory a-1 we proved the Mean Value Theorem for Analytic Functions, 

i.e. 

(64) f{zo) = ^ J 0 f( z o + rC e )dQ, 

assuming that f(z) is analytic in a neighbourhood of the closed disc B [ 2 0 ,r] of radius r > 0. 

Clearly, we obtain a similar Mean Value Theorem for Harmonic Functions by simply taking the real 
part of (64). We shall more precisely state this as a theorem, 


Theorem 2.2.2 Mean Value Theorem for Harmonic Functions. Let u{x,y ) be harmonic in an open 
domain Cl containing the closed disc B [zo,r] of radius r > 0. Then the value u(zq) = u(xo,yo) at the 
centre of the disc is equal to the mean value of u{x,y) over the circle of centre (cco, J/o) an d radius r, 

i.e. 


(65) u (z 0 ) = — 


p2tv 

/ u (zo + re 10 ) dO. 

J o 


Finally we prove that if Cl is an open bounded domain, and h{x,y) is a continuous function on the 
boundary dCl of Cl, where dCl is composed of continuous and piecewise differentiable curves, then the 
boundary value problem 


f rflu 3 2 U 

Au =Jv + Jv = 0 ’ for (x, y) G Cl, 

dx °y 

u(x, y) = h(x, y), for (a;, y) G dCl, 

has at most one solution. 


Theorem 2.2.3 Assume that Cl is an open and bounded domain. Let u and v be continuous functions 
on the closure U and harmonic in Cl itself. If u = v on the boundary dCl, then u = v in all of Cl. 

Proof. Put y = u—v. Then ip is harmonic in Cl and varphi = 0 on dCl. It follows from Corollary 2.2.1 
that since ip has its maximum and minimum (both = 0 ) on the boundary, we must have p = u — v = 0 
in all of Cl, hence u = v in Cl. □ 
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Example 2.2.1 We shall here show that we cannot exclude the assumption that Cl is bounded in 
Theorem 2.2.3. 

It follows from the obvious fact 
d 2 d 2 

A(xy) = dx 2 ^ + dy 2 '^ XV) = ° 

that u(x, y) = xy = ^ 3z 2 is harmonic in M 2 . Let Cl = {(re, y) G I 2 | y > 0} denote the upper half 
plane. Then the boundary dCl is the X-axis, and we have obviously zi(x,0) = 0 on the X-axis. 

Clearly, u{x,y) has neither a maximum nor a minimum in Cl, because its range is R. 

Since v(x, y) = 0 is another harmonic function, for which v(x, 0) = 0, this example illustrates that if 
Cl is not bounded, then the solution of the boundary value problem (66) is not unique. 

We shall later in Section 2.4 explicitly solve (66) in the special case, when Cl is an open disc. 0 


2.3 The biharmonic equation 

In problems from the two-dimensional elasticity theory one has to deal with the biharmonic equation 


( 67 ) AA4>=|^+2 


<9 4 <t> <9 4 d> 


= 0 . 


etc 4 dx 2 dy 2 dy 4 

We shall in the following solve it in simply connected domains. 


Remark 2.3.1 The conventional name “biharmonic equation” is misleading, because there are bi¬ 
harmonic functions which are not harmonic, while all harmonic functions trivially are biharmonic. A 
better name would therefore be “semiharmonic equation”. It has, however, become customary to call 
it the biharmonic equation, so we shall stick to this notation. 0 

Theorem 2.3.1 Let Cl be an open simply connected domain. Every solution of the biharmonic equa¬ 
tion AA4> = 0 can be written in the form 

(68) <S> = M{zf(z)+g(z)}, 

where f and g are analytic functions in Cl. 


Proof. As mentioned above, every harmonic function is also biharmonic. Since Hlg(z) is harmonic, 
we shall only prove that ft{z f(z)} is biharmonic to conclude that every function <I> of the form (68) 
is biharmonic. 

Put / = u + iv. Then 

5ft {z f(z)} = 5ft{(x — iy)(u + iu)} = xu + yv. 
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When (68) is inserted into (67), we get 
AA<f> = A 


d 2 d 2 

(xu + yv ) + --—(xu + yv) 
dy 2 


= A 


dx 2 
d 

dx 


u + x 


du 

dx 


d 2 i 


d 2 i 


^ dx 2 dy 2 dy 


v + y 


dv 

dy 


. i n du d 2 u d 2 v d 2 u dv d 2 v 

A { 2 —— \- x —- j + y 2 + x 71 2 ^ o—I" y 

dx dx z dx z dy z dy 


dy 2 


5u 


5a; 


dv 


= 2A < — > + 2A< — > + AjxArt + y Aw} 


5y 


O A 

= 2 — {Au} + 2 — {Au} + A(x ■ 0 + y ■ 0) = 0, 
dx dy 

and the claim is proved. 

Then assume that $ is biharmonic in 12, i.e. AA<f> = 0. We shall prove that $ can be written in the 
form (68) for some analytic functions / and g. 
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It follows readily from AAd> = A(Ad>) = 0 that Ad) is harmonic. Thus there exists an analytic 
function h in the simply connected domain Cl, such that Ad* = dlh. Using once more that Cl is simply 

connected we can define the indefinite integral / of the analytic function - h in Cl, thus f = —h. 
Then by a similar computation as above, 

= $th = Ad>, 


so 


A {d>-5i (*/(*))} = 0, 

and we have proved that d> — z f(z) is harmonic in Cl. Thus there exists an analytic function g on Cl, 
such that 

d >-3t{zf(z)} = 3tg(z), 
and (68) is obtained by a rearrangement. □ 
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2.4 Poisson’s Integral Formula 

Assuming that is an open bounded domain and h{x , y) is a continuous function on the boundary 
<9fi, we proved in Theorem 2.2.3 that the boundary value problem 

| AU =S + S= 0 ’ 

[ u{x,y) = h(x,y), for (x, y) £ dCl, 

has at most one solution. 

We shall in the following prove that in the special case of S2 = B( 0, R) there exists a solution, and we 
shall derive a solution formula. For practical reasons we shall in the following sometimes write u{z) 
instead of u(x, y ), where as usual z = x + iy. 


Theorem 2.4.1 Poisson’s integral formula (1820). Let f l be an open domain, containing the closed 
disc 1?[0,.R] of centre (0,0) and radius R > 0, and let f(z) be analytic in f l. Then we have for any 
point zq = xq + iyo = re* e £ B( 0, R), thus 0 < r < R, 


(69) u (z 0 ) = ^ J q 
and analogously, 


R 2 — r 2 


o R 2 +r 2 — 2 Rr cos(0 — t) 


u[Re lt ) dt, 


f 2-7T 


' ( 0 o) = dT 


R 2 — r 2 


27r J o R 2 + r 2 — 2 Rr cos(0 — t) 


v(Re lt ) dt. 


If u is harmonic in f l, then every harmonic conjugated function v(x,y) of u{x,y) is given by the 
formula 


1 


(.27T 


2 Rr sin(0 — t) 


(70) v (z 0 ) ^ R2 + r2 _ 2Rr cog(0 _ ^ 


u(Re lt ) dt + w(0). 


Finally, 


z+ z 0 


0 z ~ z 0 


i(z)dt + iv( 0), where z = Re lt . 


Proof. First, by Cauchy’s integral formula, 


(72) f(z 0 ) = 


1 

27T7 


I AA d2 . 

J\z\=R Z~ z 0 


If zq / 0, then z\ := R 2 /zq is a point outside |z| = R, so f(z)/{z — z±} is analytic in an open set 
containing the closed disc B[0,R]. Hence by Cauchy’s integral theorem , 

(73) 0 = —(f d z, where z\Zq = R 2 . 

27r * J\z\=rZ-Zi 

Notice that Arg zq = Arg z\ = 0 O - 
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Let z = z(t) = Re lt , t £ [0, 2i r], be a parametric description of the circle \z\ = R. Then dz = iRe lt d t = 
izdt. Writing f(z) = u + iv we get by insertion into (72), 


P4) /(*>)= A 


f27T 


u + iv 1 

*2 dt = — 


r 2ir 


(u + iv) dt. 


2ni J o z — Zq 2n J q z — Zq 

Apply the same substitution in (73) and conjugate the result. Using furthermore that 

R 2 


zz = R 
we get 


and 


2 2 

zi = — = —, 

ZQ Zo 


0 = 27T 


1 fz 


z — 


z z 
Zo 


(it — iv) dt = — 


1 / Zo 


thus when we multiply by —1, 

1 r 27r 

(75) 0 = — / ~° (u - iv) dt, 

^ 7r Jo Z — Zo 


2 tt Jo z 0 - z 


z = Re lt . 


(u — iv) dt, 


Formula (75) was proved, assuming that zq ^ 0, and it is trivial for Zq = 0. Hence we get by adding 
(74) and (75), 


/(2o) = kl {AA“ + ” ) + 


dt 


(76) 


' Z + Zo , , , * f . , , 

2 nj 0 J^o U{z)dt+ 2nl V W dt ‘ 


Now, z = Re lt and z 0 = re , so be get by (45) in Section 1.6 that 

, , - 1 + - 1 1 ' 1 2 
z + *o _ R _ 

r ,, a ^ 

1 + 


Ms)' 


2 — sin(0 — t) 

R 


z - z 0 


1 — — 
R 


— \ -2 — 
.R) R 


(77) 


R 2 - r 2 


+ i ■ 


cos(0 — t) 1 + { 


2i?r sin(0 — t) 


— ) -2 — 
R) R 


cos(0 — t) 


R? _|_ r 2 _ 2 Rr cos(0 — t) R 2 + r 2 — 2Rr cos(0 — t) 
Finally, it follows from the Mean Value theorem for Harmonic Functions that 


-| /»27T -| /»27T 

2ii v{z)it= ^L »(Se“)d< = «( 0), 


and (71) follows from (76). If (77) is put into (76), we get (69) and (70) by splitting into the real and 
the imaginary parts. □ 
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Example 2.4.1 Let = B( 0,1), and let h(z) be a continuous function on the unit circle \z\ = 1. We 
shall solve the boundary value problem , 


(78) 


A u = 0, for (x, y) £ f 1, 

u(x, y ) = h(x, y) for x 2 + y 2 = 1. 

First define a continuous function ip(t) on [0, 27 t] by 
<p(t) = h(e lt ) , t £ [0, 2t r], </?(0) = ip(2n). 

Then by the classical Theory of Fourier Series , y>(f) has a Fourier series expansion, 

+oo 


| i 

<*>(*) ~ x »o + E {a n cos nt + b n sin nt} , 


where 


a n = ~ 

7T 


[ <p(t) cosntdt and b n = — [ 

Jo 77 Jo 


ip(t) sinnf dt. 


The solution of (78) is given by (69), 


u(re 


1©^ _ 


1 


1 — r 2 


2-7T J 0 1 + r 2 — 2r cos(0 — t) 
We have in Section 1.6 proved (45), i.e. 

2 +°° 


V5(t)df, r £ [0,1[. 


(79) 


1 — r 2 


1 + r 2 — 2r cos(0 — t) 


— = 1 + 2 rn cos(n{0 - t}). 


71=1 


When 0 < r < 1 is fixed, then the series in (79) is uniformly convergent, so we can interchange 
summation and integration in the computation below, when (79) is inserted into the expression of 
u(re l& ), 


u(re' e ) = 


J, 


+oo 


dt + ‘ 


1 


<p(t) cos(nO — nt) d t. 


71=1 


Here, 


r 2n i 

J ^ <p(t)dt=-a 0 , 

ind 

/»271 /»27T /»27T 

- / ip(t) cos(nQ—nt) dt = — / (p(t) cosnf dt-cosn0+ — / <p(t) sinnf dt-sinn0 = a n cos n<d+b n sin n0, 

71 Jo 71 J o 7r Jo 


27T 


hence 


+oo 


(80) n(re* e ) = - oq + ^ r 11 {a„cosn0 + &„sinn0} , 


71=1 


and we have proved that the unique solution of the boundary value problem (78) is obtained from the 
Fourier series expansion of ip(t) by multiplying the n-th term of the series by r n , r £ [0,1]. 0 
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2.5 Electrostatic fields 

A planar electrostatic field is a vector field in the complex plane, which at each point of its domain 
is given by the vector £ ~ E x + iE y , corresponding to the force, which the field would exercise on a 
unit charge at the point. 



Figure 6: The tangent field and the normal field of a closed curve C in the plane. 


The flux through a closed curve C of the field is defined by 
N := j) £ ■ nds = 47re, 

(cf. Figure 6), where e is the sum of all charges inside C , and n is the outward normal vector field of 
the curve C, and s is the natural parameter of C defined by the curve length. At each point z, 

m tog ** e , 

where the curve C in some sense shrinks to z (choose e.g. C r as the circle of centre z and radius r, 
and let r —> 0+), and where S is the area of the bounded domain inside C 1 and finally, q is the density 
of the charge at the point z. 

The circulation of £ along C is equal to the work 
W := j) £ ■ tds, 

where t denotes the unit tangent field of C. When the work is 0 along every simple closed curve C, 
we get by Stokes’s theorem in two dimensions that 

8E V dE x 

rot £ = - —- = 0. 


dx 


dy 


This shows that the differential form E x da; + E y dy is closed, thus exact in (at least) every simply 
connected domain. Hence, there exists a potential v, such that 

E x da; + E y dy = —daa, 
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(the minus sign is due to convention), so 


dv . dv 
dx dy 


—grad v. 


If an open and simply connected domain Cl does not contain any charge, then it follows from (81) that 


div £ = 


dE x . dE, 


dx 


dy 


= 0 , 


so the differential form —E y dx + E x dy = du is also exact. Then it is easy to prove the the level 
curves of the function u(x , y) are the field lines , where the tangents define the direction of the field. 


Under the given assumptions above we have constructed u and v, such that 


du 

dx 


= -E. 


y> 


du 

dy 


— E x 


°1--E ---E 

dx~ x ’ dy~ v ' 


It follows that u and v satisfy Cauchy-Riemann’s equations , so they are harmonic functions in U, and 
the function 


f{z) = u(x, y) + i v(x, y) 


is analytic in Cl. 
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The function f(z) is called the complex potential of the field, and the electrostatic field £ can be 
represented by 


dv . du 
dx 1 dx 


-if'(z)- 


The field lines and the potential curves are level curves of harmonic functions. This shows that we also 
for (two-dimensional) electrostatic fields are interested in the solution of the boundary value problem 
( 66 ). 


2.6 Static temperature fields 

Let u(x,y,t) be a planar temperature field, in a domain il. If f l does not contain any source of heat, 
then it can be proved that u(x, y , t ) satisfies the heat equation 


<»>!= 


d 2 u d 2 u 'I 
dx 2 dy 2 / ’ 


where t denotes the time, and a 2 is a positive constant. We shall not here solve (82) in general. We 
only note that if the temperature field does not depend on time (corresponding to an equilibrium state 
of the temperature), then u(x,y,t) = u(x,y), and (82) reduces to the well-known Laplace equation 
A u = 0, and u(x, y) is an harmonic function. 
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f(z) 

u = 3t/(z) 

v = Qf(z) 

if(z) 

-v = -9/(z) 


u = 5R/(z) 

z 

X 


y 

z 2 

2 2 

2 T — y z 


2 xy 

z 3 

x 3 — 3 xy 2 


3 x 2 y — y 3 

1 

X 


y 

z 

x 2 + y 2 


(N 

+ 

(N 

e z 

e x cos y 


e x sin y 

sin z 

sin x • cosh y 


cos x ■ sinh y 

COS 2 

cos x • cosh y 


— sin x ■ sinh y 

sinh 2 

sinh x ■ cos y 


cosh x ■ sin y 

coshz 

cosh x ■ cos y 


sinh x ■ sin y 

tan z 

sin 2x 


sinh 2 y 

cos 2x + cosh 2 y 

cos 2x + cosh 2y 




Arg z = 0 

for 0 £ } — 7T, 7r[, 




A X 

Arccot — 

y 

for y > 0, 

Log z 

In | z| = In r = — In ( x 2 + y 2 ) 

< 

Arctan — 

X 

for x > 0, 




Arctan — 

X 

7r for a? < 0 and y < 0, 




A x 

Arccot - 

y 

7T for a; < 0 and y < 0. 


Table 1: Some elementary analytic and harmonic functions. 


70 


Download free eBooks at bookboon.com 













Calculus of Residua 


Laurent Series and Residua 


3 Laurent Series and Residua 

3.1 Laurent series 

We previously proved in Section 1.6 that if we have given the convergent power series 

+oo 

(83) f(z) = ^a n z n , for \z\ < i?, 

n= 0 

then we obtain a Fourier series by putting z = re 10 , where 0 < r < R and 0 € R. In fact, 

+oo +oo 

(84) ^(0) =f(r e 10 ) = £ a n r n e inQ = ^ c n e ™ 0 

n=0 n= 0 

This is a special case of a Fourier series, because we usually sum from —oo to +oo in the Theory of 
Fourier series. In other words, we are missing all terms of the form c_ n e -m0 , n £ N, in (84). From 

c- n e~ m& = c n r n • (re‘ e ) " = a_ n z~ n 

follows that we miss all terms of the form a n z~ n in (83). Thus, the general Fourier series force us to 
consider the more general functions of the form 

+oo 

(85) f(z)= ^2 a nZ n , 

n=— oo 

where we still have to discuss where this series is converging, and is representing an analytic function. 
Before we start on this project we sketch a useful application in the technical sciences of series of the 
form (85). 

Let u(t), t > 0, be a continuous function in time t, and assume that it is measured at the equidistant 
times t = nT , n£N. In this way we define a sequence ( a n ) by 

a n = u(nT), n C No- 

Such sequences are used in theoretical considerations in e.g. Cybernetics. One uses the so-called 
3 -transform, which is defined by 

+oo . +oo 1 

(86) 3t(u)(z) := ^ a n ■ — = ^ u{nT) ■ —, 

71—0 ^ 77=0 

in a domain, where ( 86 ) is convergent. In this case all exponents are non-positive. 

The 3 -transform may be considered as a discrete form of the Laplace transform. We shall in more 
details return to them in Ventus: Complex Functions Theory a~4- For the time being we have only 
used ( 86 ) as a motivation and excuse for also to consider series with negative exponents. 

We shall return to mathematics. 

First consider a function f(w), which is analytic in w for |w| < r, so f(w) has a convergent power 
series expansion 

+oo 

(87) f(w) = ^ b n w n , M < r. 

n— 0 
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Figure 7: A continuous function in time, which is measured at the equidistant times nT = n, n £ No, 
where we have normalized, such that T = 1. 


Let zq £ C be fixed, and put w = g(z ) = -, thus |u>| = T - T < r for \z — zq\ > -, and g(z) 

z-z 0 \z-z a \ r 

is analytic in this domain. Hence, the composite function (/ o g)(z) is analytic for \z — Zq\ > -, and 

1 r 
when w = - is put into (87), we get the convergent series expansion 


z- z 0 

(88) h(z) = (/ o g)(z) = y^b n -(z- z 0 )~ 


+oo 


for \z - zq | > -, 


n=0 


where the series is uniformly convergent in each closed subset 4C!1 = C\ B 
not necessary here to assume that A is compact, because g° - 1 (A) is compact. 


A>, 


Note that it is 


Since w = 0 corresponds to z = oo, it is natural to say that h(z) is analytic at oo, and we put 
h{ oo) = bo. 

If r = +oo, then the series ( 88 ) is convergent for all z zq. If r = 0, then the series is divergent for 
all z 7 ^ oo, i.e. it is only convergent at oo with the value bo¬ 
lt follows from the above that a series ( 88 ) with only non-positive exponents usually is convergent in 
the complement of a closed circle, C*\B zq, - , where C* = C U {oo} denotes the extended complex 

plane. If r = +oo, then the domain of convergence is C* \ {zo}, an( i if r = 0, then the series is only 
convergent at oo. 


Definition 3.1.1 A Laurent series expanded from zo G C is a series of the form 

+oo +oo +oo 

(89) ^2 a n (z - z 0 ) n := ^2 a n (z - z 0 ) n + ^2 a- n (z - z 0 )~ n ■ 

n =—oo n —0 n —1 

Its domain of convergence is the intersection of the two domains of convergence of the two series on 
the right hand side of (89). 
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If the “Taylor series” a « ( z ~ z o) U is convergent for \z — zq\ < R, and the series of terms with 

negative exponents a -n ( z ~ -o)~” is convergent for \z — zo\ > r, and if furthermore r < R, 

then the domain of convergence for the series (89) is given by the annulus r < \z — zq\ < R. 



Figure 8: The domain of convergence of a Laurent series is an annulus. 


Note that if r = 0, then the domain is B ( zq,R ) \ {-So}, where only the centre z o has been removed 
from the disc (a deleted neighbourhood of £o), and if R = +oo, then the domain is the complement 
of a closed disc, C\B [z 0 , r}. If r = 0 and R = +oo, then the domain is of course the deleted complex 
plane C \ {^o}. 

If r > R, then the Laurent series is divergent , and (89) does not represent a(n analytic) function. 

We shall prove a theorem which is analogous to Theorem 1.3.1. 
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Theorem 3.1.1 Laurent’s theorem. Assume that f(z) is analytic in the open annulus 
Ll = {z £ C | ri < \z — zq\ < r 2 } 1 where 0 < r\ < r 2 . 

Then f is uniquely determined in Ll by its Laurent series 

+ OO 

(90) f(z) = ^ a n {z- z 0 ) n , 

n— — 00 

where 

(91) a n = (f -— +1 dz, n £ Z. 

2 m Jc (z - z 0 ) 

Here, C is any simple and closed curve in Ll which separates the two boundaries \z — 2 ol = r i an d 
\z - z 0 1 = r 2 . 

The series (90) is uniformly convergent in every compact subset of LI. 

PROOF. Let C(r) denote the circle \z — zq\ = r. If rl < R\ < R 2 < r 2 and n € Z, then it follows 
from Cauchy's integral theorem for multiply connected domains that 

i m n+1 d z=i dz. 

JC{RP) {z - Zo) + Jc(R 2 ) (z - Zo) + 

so it suffices only to consider C = C(r) in (91). 

For every n £ Z the constant a n is uniquely determined by 

■j^ /» j? 

a n = <f> ~— -“wry dz, where ri < r < r 2 . 

2' 7r * Jc(r) (Z - Zo) 


We shall prove that the corresponding series Y^n =-00 (~ — 2 o)™ is convergent for every z € Cl, and 

that its sum function is f(z). 

Fix z C LI, and choose R± and R 2 , such that 

ri < Ri < \z — zq\ < R2 < i~2- 

Then by Cauchy’s integral formula, 


(92) f(z) 


1 

27 ri 


i m. 

Jc(R 2 ) C - 2 


dC 


1 

2ni 


/ m 

JC(Ri) C Z 


dO 


If C S C (Rf), then |£ 


(93) 


1 


1 

C - Zo 


— Zq\ = R 2 - From \z — Zq\ < R 2 follows that 


1 

z - Zo 

C - z 0 
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Figure 9: The paths of integration in the proof of Laurent’s theorem. 


where the series is uniformly convergent for £ £ C (R^)- Hence, (93) can be put into the former 
integral on the right hand side of (92), and we are allowed to interchange summation and integration, 
thus 


+oo 


(94) h{z) = ^ S dC = £- zo)n '^ 


/(C) 


+oo 


Jc(r 2 ) (C - z 0 ) 


n+1 


dC = 51 a n (z - zo) n ■ 


n —0 


'G(R 2 ) C z n=0 

The latter integral of (92) is treated similarly. If f £ C (Ri), then |£ — zq\ < \z — zq\, hence 

+oo 


(95) - 


1 


1 


1 


— =—y 2 \— 

— Zr\ 7 — 7.r\ • ^ I 7. — 


C - z z - z 0 1 _ C - z o Z - z 0 n=o 


C - z 0 


Z- z 0 


Z- z 0 

which is uniformly convergent for f £ C (R\). Then 

+oo 


1 


h(z) = /^-dC = g--(C - *o) n /(C) dC 

2tt* Jc{r{) C - z n=Q (z - z 0 ) + 27 n Jc^) 


(96) 


= £ (z - z 0 ) n ■ — 


/(C) 


2tt* JCiRi) (C- z 0 ) 


n+1 


dC = 5Z Un ( z ~ Z °Y 


We get (90), when (94) and (96) are put into (92). 


The series (94) and (96) are both uniformly convergent on every compact subset of fl, so we have 
proved the theorem. □ 

It is in particular easy to find the Laurent series of rational functions. We shall in the following give 
some examples which show the technique. 
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Example 3.1.1 The function f(z) = -is analytic in C\{2}. Seen from zq = 0 the function f(z) 

z 2 

is analytic in the open disc fli = B( 0, 2), and in the complement = C \ 1?[0, 2] of its closure. We 
shall find the Laurent series of f(z ) in fli and in fl 2 . In both cases we apply the geometric series in 
an essential way. 



3- 

_—a- 

Omega_2 


/ Omega_1 


-3 “I 

-1 

-j- ^ 3 


'"=2- 



-3- 



Figure 10: The two domains and 12 2 , in which we have a Laurent series expansion of f(z) = 
from zq = 0. 


First consider z £ fli = B( 0,2). Then |^:| < | — 2| = 2, so 

-1 +oo +oo 1 


< 1, and we get 


(97) 


1 


z-2 


1 




2 :{|}~£ 


2 n+1 


z , 


n —0 n —0 

and f(z) is in described by its Taylor series expanded from zq = 0. 

2 


Then let z £ fl 2 . In this case, |z| > | — 2| = 2, so 


< 1, and 


(98) 



2 +oo 


n —0 


+oo 1 

-!• =Sf 2 ” _ 1 -F- 

n—\ Z 


corresponding to that f(z) is given by a Laurent series in fl 2 . Since the Laurent series expansion is 
unique according to the theorem above in each of the domains, the problem is solved, and we have 
the descriptions given by (97) and (98), i.e. 


f(z) = | 

The Laurent 


E +°° 1 , r n 

™=o 2«+i ' " ’ 


E + OO iyn—1 
n=1 Z 


series expansion from zq 


for \z\ < 2, 
for |z| > 2. 

0 does not make sense on the circle \z\ = 2. 0 
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Example 3.1.2 The function f(z) = 


is defined and analytic in C\{1,2}. If we choose 


(z-l)(z-2) 

z o = 0, then the Laurent domains are the disc fii = B( 0,1), the annulus fl 2 = {0 € C | 1 < \z\ < 2}, 
and the complement fl 3 = C \ B[ 0,2] of a closed disc. 



Omega 3 






Omega_2 


/ /-'I- 



/ 

)mega_1 \ 


4 4 z - 

’ ,1 

' i 

j 

~ “2- 




Figure 11: The three Laurent domains for f(z) 


1 

(z-l)(z~ 2 ) 


and zq = 0. 


First decompose 


m = 


( 0 — 1)(0 —2) z — 2 z — 1 


Notice that we have already found the Laurent series of-— in Example 3.1.1. 

If z e fli = B( 0,1), then \z\ < 1, so 

m = 


m = 


m = 


l 

1 


1 

0-2 

z — 

1 

2 1 

i.e. 1 < 

1*1 < 

2, then 

1 

1 


1 

CM 

1 

z — 

1 

2 1 

then 0 

>2, 

and 

we get 

1 

1 


1 1 


+oo 




1-0 


n— 0 


2 n+1 


11 +°° i 

_I . = _ V — 0" 

£ 0 , 1 ^ 11 

r\ 1 — — n —0 

z Z 


+oo 


2 n+1 


- V-. 

Z_^ 


+oo 


+oo 


0-2 0-1 0 i_ 2 


'Er'-'iv'Er 1 - 1 }'? « 


1 — — n—1 
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Example 3.1.3 The function f(z) = 


1 


is analytic in C \ {1}. The denominator (1 — z ) 2 has 


(1 . -! 2 

the root z = 1 oi multiplicity 2, so the direct determination of the laurent expansions from zo = 0 
in fli = B( 0,1) and Q 2 = C \ B[ 0,1] becomes more difficult than in the previous two examples. The 

trick is instead first to find the Laurent series of g(z) = -, where the denominator only has a 

simple root, and then find the Laurent series of f(z ) by termwise differentiation. 

We find in Lli = B( 0,1), 


ff(*) = 


1 


+ OO 

= E ; 


n—0 


l-z 

and in fi 2 = C \ B[ 0,1], 


N <1, 


9(z) = 


1 

l-z 


+OO 


+ OO 


1 - 


= -E^ = -E' 


Izl > 1. 



360 ° 

thinking. 


Deloitte. 
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Clearly, f(z) = 


1 


+oo +oo 


= g'(z ), so we get by termwise differentiation, 


/(z) = ^rt/- 1 = ^(n+ l)z n for zgSl 


l) 


and 


n —0 


+oo +oo 

„— n — 


f(z) = J2n-z- n - 1 = J2{n-l)z- n for z e Sl 2 . 


n • z 

n =1 n—2 


This method can of course be generalized to h„.(z) = ^- , in which case we apply that /i(z) = 


1 

n! 


(i - z y 


9 {n) (z). 0 


3.2 Fourier series II 

Assume that f(z) is analytic in an open annulus r± < |^ I < J" 2 , and let f(z) in this annulus be 
represented by its Laurent series 

+oo 


f(z) = a n(z~Zo) n , ri < \z — z 0 | < r 2 . 


Let ri < r < r 2 and 0 € R, and put z = zo + r e l& . Then 

+oo 

(99) f( Zo + re 10 ) = £ a n r n e in& , 


n =—oo 


where the series in (99) is uniformly convergent in 0 for fixed r. 

We see that for given zq and r the series expansion (99) is the Fourier series of the function 
¥>(©) := / (zo + re’ e ) . 

This is in agreement with Laurent’s theorem, because 


dn, — 


1 C/(*. + -») ire ,e dei 


2m J o r n+1 e*(” +1 ) e 
from which 

1 r 2v 1 r 2 

/(=» + ”' e )«-” ede = 2 rj 0 


ip(Q)e~ in& d0 


which is the usual formula for the n-th Fourier coefficient of <p(0). 

We proved in Section 1.6 ParsevaVs formula in a special case, necessary for the proof of the maximum 
principle. We can now easily prove the general ParsevaVs formula. 
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Theorem 3.2.1 Parseval’s formula. Let 

+oo +oo 

fi z ) = ^2 a « and g ( z ') = hn zH 

n= — oo n =—oo 


be analytic in the annulus rq < |z| < V 2 - Then for every r £ ]ri,?" 2 [, 

/.2-n- _ +oo 

(100) — J f (re 10 ) g(re ie ) d<3 = ^ a n b n r 2n . 


n ——oo 


Proof. The only difference from the proof in Section 1.6 is that we here sum from — oo to +oo 
instead of from 0 to +oo. □ 


3.3 Solution of a linear differential equation by means of Laurent series 

We gave in Section 1.4 a solution procedure for a linear differential equation of order n with analytic 
coefficients in an open domain 12, 

slTL •£ 777 ,— 1 p 

( 101 ) a 0 (z) + ai{z) d ^ n _\ H-b a n -i(z) f(z) + a n (z) f(z ) = g(z), z £ i 2, 

where the coefficient ao(z) of the highest order term is not identically zero. If a 0 (zq) yf 0 for zq £ 12, 
then we can apply Theorem 1.4.1. However, if z 0 £ 12 is a singular point for the differential equation 

(101) , i.e. a 0 ( 20 ) = 0, then we do not have a result, which guarantees that there exists an analytic 
solution in a (deleted) neighbourhood of zo ■ 

One may assume that (101) in the singular case has a Laurent series solution of the form 

+ OO 

( 102 ) f(z) = ^2 a n(z- ~ 0 )” , 


where the Laurent series expansion of f^{z) is formally found by termwise differentiation, and where 
Zq clearly does not belong to the domain of convergence of /(z), which is still to be found. In particular, 

+OO 

(103) f(z)= ^2 na n (z - z 0 ) n ~ l 

n =—00 

in the possible domain of convergence. Since zq is excluded from this domain, we do not have an 
exceptional case in (103) for n = 0 as we did in Section 1.4, where we only dealt with power series. 
Therefore, we shall no longer be careful with the domain of summation, which always may be put 
equal to Z. Hence, if we are looking for a Laurent series solution, there is no need to specify the 
bounds of summation, and we just write ^ instead of 

We use the same method as in Section 1.4 to set up a recursion formula for the coefficients a n in 

(102). The new aspect here is that this recursion formula must be used to find both the coefficients a n 
for n > 0, and the coefficients o_ n for n > 0. Since the recursion formula for a_ n of negative indices is 
used in the “opposite direction” of the customary one, it is highly recommended for the novice of this 
theory (it is, however, not mandatory) to put b n = a- n for n > 0 and then find a (usual) recursion 
formula for the b n . 
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A small and useful, though often neglected detail here is to start by finding the values of n £ Z, 
for which at least one term in the recursion formula becomes zero. Such an observation will give us 
some information of which coefficients are zero, and which coefficients are arbitrary. Then the formal 
Laurent series is expressed by these arbitrary coefficients. 

Finally, when we have found all the formal Laurent series (102) which satisfy the equation (101), we 
must find its domain of convergence!!! This point is of paramount importance and must never be 
forgotten, because otherwise the computations do not make sense. 

This point can best be illustrated by the fact that it is not hard to construct a linear differential 
equation of second order (101) for which the formal Laurent series solution (102) contains three 
arbitrary constants, which is contradicting the theory, unless the arbitrary constants are sorted out 
by a discussion of the domains of convergence, which always will force at least one of the arbitrary 
constants to be zero in any given subdomain. 

To be more precise, when we find Laurent series solutions expanded from a singular point zq of a linear 
differential equation of order re, then there are at most n linear independent Laurent series solutions, 
which are convergent in the same domain, while there may be more formal Laurent series solutions. 


Example 3.3.1 We shall give one example, showing the method. We tried previously in Exam¬ 
ple 1.4.3 to solve the linear differential equation 

(104) z 2 f\z) - f(z) = -z. 

We shall here try to solve it by means of Laurent series. 

We first solve the corresponding homogeneous equation 

(105) z 2 f\z) - f(z) = 0. 


First note that it easily follows from the rearrangement J for f(z) ^ 0 and a check that the 

/(-) -2 

solutions are 


(106) f(z) = c-exp - 


z 7 ! 0 and c £ C arbitrary. 


We shall now prove (106) by using formal Laurent series instead. So assume that f(z) = )T) a : n z n 
is given by a formal Laurent series. Then by insertion into a reversed (105) 


0 = z 2 f(z)-f(z) = z 2 J2 

(107) 


na n z 


n— 1 


-E«^ n = E 


na n z 


n+1 


-E 


a n z 


= ^2 na nZ n+1 - a n +iz n+1 = - a n + 1 ) z n+1 , 

so we get the recursion formula 
(108) a 7 j + i = na n for n € Z. 

An obvious “zero” of (108) is n = 0, for which value 


ai = 0 • ao = 0, 
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no matter how a o is chosen. Hence, a o is for the time being an arbitrary constant, while a i = 0. 

We shall exclude n = 0 in the following. This case was discussed above. The investigation is then 
split into the two cases, n > 0 and n < 0. 

When n > 0 and a\ = 0, it follows by recursion of (108), or by a division by n!, from which 

1 1 1 
T a n+l 7 7TT a n ' ' ' — 77 a l — 0, 

n\ (n — 1)! 1! 

that a n = 0 for all n £ N. 

Then assume that n < 0 and put m = —n > 0 and b m = a_ n . Then (108) becomes 

a n +l — ft—m+l = ^—{m— 1) = 1 = ri d n — 1TI d—m = 

so we get the recursion formula 
(109) 6 m _i = —mb m for m £ N. 

Notice that m = 1 corresponds to — b\ = bo = ao, where ag is arbitrary. 
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If (109) is multiplied by (—l) m 1 (m — 1)! ^ 0 and we read from the right to the left, then we get by 
recursion 

(-1 ) m rn\b m = (-1 ) m ~\rn - 1 )!6 m _! = • • • = (-l)°0!6 o = «o, 
hence 


(-l) m 

Qj—m bm — | flO) rn £ N, 

m\ 

which of course also holds for m = 0. 


Thus we have derived the formal Laurent series solution of (105), which is given by 


+oo 


0 + a_„ 


m —0 


+°° / IV 

_ \ ' W 1 ! 

< ?7,! 

n =0 


n\ 


a 0 


+oo 

= a,, V' 


= ao 


^ n\ 

n —0 


= a 0 exp - 


where we have recognized the coefficients of the series expansion of the exponential function. The 

latter is of course convergent for all ;€C, for which-eC, thus for 2 £ C\ {0}, and the domain of 

z 

convergence is C\ {0}. There is no need to use any other procedure of finding the radii of convergence 
(r = 0 and R = +oo). 


We have seen that equation (105) could be solved by guessing a Laurent series solution. However, if 
we try the same method on the inhomogeneous equation (104), then it follows by a modification of 
the computation of (107) that 

-z = z 2 f'{z) - f(z) = ^2(na n - a n+ i) z n+1 . 

Here, n — 0 corresponds to z n+1 — z, so we get the recursion formula 

{ n-a n = a n+ i, for n ± 0 , 

—a\ = — 1 , for n = 0 . 

The only change from (108) is that now ai = 1. This change, however becomes disastrous, because 
then it follows from the first recursion formula above that 

a n +i = na n = n(n — l)a„_i = ■ • • = n\ oi = n! for n £ N, 

and trivially ai = 1 = 0!. Hence by a translation of the index, a n = (n — 1)! for n £ N. Since there 
is no change in the recursion formula for n < 0 , we again get the previous solution, so in this case all 
formal Laurent series solutions are given by 

+oo +oo 1 

X>-i)n" + a„£-- i (-irv” 

n= 1 n =0 

The latter series is of course convergent for z / 0 with the sum function ao exp ^^, but the former 

series has 0 as radius of convergence, so it is only convergent for z = 0 (where the latter series is not 
convergent, unless a o = 0 ). 


We conclude that there does not exist any convergent Laurent series solution of (104), when expanded 
from zq = 0. However, if we move the expansion point Zq away from 0, then it follows from Theo¬ 
rem 1.4.1 that there will always exist even power series solutions in the open disc B(zo, |-o|), z o 7^ 0. 
0 
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3.4 Isolated boundary points 

We have seen that we may have Laurent series which are convergent in an annulus “far away” from 
the expansion point zq. It should, however, not be of any surprise that the most important case is 
when the domain of convergence is a deleted disc D(zq, R ) = B(z 0 , R) \ {z 0 }, in which case the centre 
Zq is an isolated boundary point of the domain of convergence. In fact, it will lead us to the important 
Calculus of Residua. 

Let / : Cl —» C be analytic in an open domain Cl, and let zo G C be an isolated boundary point of Cl. 
This means that zq Cl and that we can find R > 0, such that the deleted disc 

D(z 0 ,R) :=B(z 0 ,R)\{z 0 }CCl. 

In this case we can apply Laurent’s theorem to get 


+OO 

(110) f(z) = ^ a n (z- z 0 ) n , for z G D(z 0 , R ), 

n =—oo 

where 

(111) a n = - — ~^ n+1 dz, for all r G ]0, R[, 

JC{z 0 ,r) [Z - Zq) 

where C (zo, r) denotes the circle |z — zo| = r of centre zq and radius r, traversed in the positive sense 
of the plane. 

We shall consider three cases. 


3.4.1 Case I, where a n = 0 for all negative n. 

In this case /(z) is given by an ordinary power series 

+oo 

(112) f(z) = ^2 a n (z - z 0 ) n , z G D(z 0 , r) . 

n =0 

The power series is clearly convergent and analytic in the full disc B(zo, r), so we extend /(z) by (112) 
to all of B(zo,r), i.e. we add /(zo) := oo to the definition. 

In this case we call zq a removable singularity. 

We have 


Theorem 3.4.1 If the analytic function f is bounded in a deleted disc D(zo,g) = B(zo,g) \ {zo}, 
then zq is a removable singularity. 


Proof. Assume that |/(z)| < M for all z G D(zq,q). 
r G ]0, o[ and every oGZ, 


Then by an application 


of (111) for any 


m 


2?r* Ic(z 0 ,r) (Z - Zq) 


n+1 


dz 


1 

< — 


M 


2n r n+1 


2tt r 


M 
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It follows from this estimate for n < 0 by letting r —> 0 + that a n = 0 for all n < 0 . □ 


sin 2 


Example 3.4.1 The function f{z) = - for 2 £ C \ {0} has a removable singularity at 2 = 0. In 

2 

fact, 


sin 2 lpZ (-1) 


+oo 


n=0 


2 2 2 4 


z t"A 2n + l )- 


2n+l _ 'y ^ J2n __ ^ __ ^_|_ “_ 


—' (2n + 1)! 


n —0 


3 ! 5 ! 


for 2/0, 


and the series is convergent for all 2 £ C, so we can put /( 0 ) := 1 . 0 


3.4.2 Case II, where a n 7 ^ 0 for a finite number of negative n. 

In this case there is a q € N, such that a_ q yf 0 and a n = 0 for all n < —q. Then in D(zo, Q ), 

+OO 

f ( z ) = V] a n (2 - Z 0 ) n = . a ~ q 4 -h a-1 + a 0 + ai (2 - z 0 ) 4 -. 

( z - * 0 ) * - *0 

The analytic function 
g{z) := (z- z 0 ) q f{z) 

has clearly a removable singularity at zq , so it can be extended to the whole disc B(zq,q) by adding 
the value g (20) := a_ g to its definition. 

In this case we say that f(z) has a pole of order q at 20, and we put / (20) = 00 (the complex infinity). 
We have 


Theorem 3.4.2 If f(z) —> 00 for z —> Zq in D(zo,g), then f(z) has a pole at z 0 of some (finite) 
order q £ N. 


Proof. This proof is not so easy as the proof of Theorem 3 . 4 . 1 , because we shall show that the 
singularity has a finite order q £ N. 


From the assumption f(z) —> 00 for 2 —> 20 follows that there exists a deleted disc D(zq,q) = 
B{z 0 , g) \ {2 0 i, such that f(z) ^ 0 , and such that the reciprocal hiz) = ——— is bounded in D(z 0 , g). 

fW 

Since h(z) —> 0 for 2 —> 20, the point 20 is a removable singularity of h(z), and h{z) can be extended 
to an analytic function in the full disc B(zq, g) by adding the value h (20) = 0. 


Since h{z) is not identically zero in B(zo, g ), we may apply Theorem 1 . 5.2 to conclude that the zero 
2 q is isolated. Thus by Theorem 1 . 5.1 the zero zq has a finite order q £ N, so 


( 113 ) h(z) = (2 - 2 0 ) 9 {b Q + bi (2 - 2 0 ) H-} 


in a possibly smaller open disc B(z 0 , gf), gi £ ] 0 , o[. Note that b 0 = 
from the continuity that g\ can be chosen, such that the latter factor 


h^ (2 0 ) ^ 0, hence it follows 
{bo + 61 (2 — 2q) + • ■ ■ } in ( 113 ) 
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is analytic and ^ 0 in B(zo, Q\). Then its reciprocal {6 0 + bi(z — 2o) + • • • } 1 is also analytic in 
B(zq,qi), so it even has a Taylor series expansion 

+ OO 

{b 0 + h (z - z 0 )-\ - y = a n (z - z 0 ) n for z £ B(zq, £?i) , 

n =0 

where a o ^ 0. Thus we conclude from (113) that 


+oo 


+oo 


/(*) = 


^2 a n(z- Z 0 ) n = ^2 a n+q (z - zo) n , z e D(z 0 , pi) = B(z 0 , £>i)\{^o} • 


h (z ) (z - z 0 ) q 

and we conclude that zq is indeed a pole of order q of f(z). □ 


z-2 


Example 3.4.2 The rational function - „ 

(z 2 + 1 ){z- l) 3 

poles (i.e. poles of order 1) at z = i and z = —i. 0 


has a pole of order 3 at x = 1, and simple 
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3.4.3 Case III, where a n 0 for infinitely many negative n. 

In this case we call zq an essential singularity for f(z). 

We mention without proof (because the full proof is very difficult) the following theorem. 


Theorem 3.4.3 Picard’s theorem (1879). If zo is an (isolated) essential singularity for f(z), then 
the image f(D(zo, g)) is either all of C, or C \ {u>o} with just one exception point wo, where D(zo, g) 
is any deleted disc contained in the domain of f. 


Picard’s theorem shows that the behaviour of an analytic function / in a deleted neighbourhood of 
an essential singularity is very wild. 

To give some understanding of Picard’s theorem we shall below prove the following weaker result. 


Theorem 3.4.4 Casorati-WeierstraB’s theorem. Let z o he an (isolated) essential singularity of an 
analytic function f. Given any deleted disc D(zq, g) = B(zq, g) \ {^o} Q contained in the domain 
C 1 of f, the image f(D(zo, g)) is a dense set in C. 


Proof. Contrariwise. Assume that we can find g > 0, such that D(zo,g) C Cl, and such that the 
image f(D{zo , g)) is not dense everywhere in C. Then we can find wq € C and S > 0, such that 

f(D(z 0 , g)) D B [w 0 , <5] = 0, 

thus | f(z) — u>o| > 0 for every z € D(zq, g). In particular, g(z) = - is bounded and analytic 

in D(zo, £0> an d ~o is a removable singularity of g(z). 



Figure 12: Proof of Casorati-WeierstraB’s theorem. 


There are two possibilities: If lim^ Zo g{z) 
that 


= x 0, then it follows from the rules of computation 


f{z) - w 0 


1 


thus f(z) —> wq + 


1 


for z —> zq, 
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and Zo is a removable singularity of f{z). 

If instead lim z ^ Zo g{z) = 0, then /(z) — wo = oo, hence /(z) —> oo for z —> zo, and Zq is according to 
Theorem 3.4.2 a pole for /(z). 

In neither of the two possible cases above zo is an essential singularity. By contraposition of the proof 
we get by negating the first sentence in this proof that for every D(zq, g) C Q the image f(D(zo , g)) 
is dense everywhere in C, and the theorem is proved. □ 


We mention two ways to show that an isolated singularity zo £ Q for an analytic function / : —> C 

is an essential singularity. 

1) If the Laurent series expansion /(z) = a n (z — zo)" expanded from zq is convergent in a deleted 
disc D(zo, g) = B(zq, g) \ {zo}, and infinitely many of the coefficients a n for n < 0 are not zero, 
then zo is an essential singularity for /(z). 

Notice that we do not require that all a n ^ 0 for n < 0, just infinitely many of the a_i, a_2, 

2) If we can find a sequence (z„) C f2, such that z n —> zo, while the limit of (/(z n )) does not exist 
for n —> Too, or “ alternatively ” (it is actually a variant of the same) if we can find two sequences 
(~n) ~ i * z o and (z") —» zo in for n —> Too, such that 

lim / (4) ^ I™ / (4) > 

n —>+oo n —>+oo 


then zo is an essential singularity for /. 

We here allow the complex infinity oo to be a possible limit. 


Example 3.4.3 The function /(z) = exp( ), z G C \ {0}, has an essential singularity at Zo = 0. 
In fact, the Laurent expansion of / is 


exp 

where 


+oo 


= E 


i r i 


n —0 


n! 1 z 2 


n +oo ^ ^ 

= E n for z C C \ {0}, 

n =0 


a n = 0 for n > 0, a_2n+i =0 for n > 0 and a_2n = —for n > 0. 

n! 


We see that a n ^ 0 for infinitely many negative n (all even negative numbers), and a n = 0 also for 
infinitely many negative n (all odd negative numbers). 


i % 

An alternative proof is to choose z'= — and z" = —. Then clearly z' —> 0 and z" —> 0 for n —> Too, 

n n 

while 


/ ( 4 ) = exp(n 2 ) —> Too and 


/( 4 ) = exp(—n 2 ) —> 0 for n —> Too. 


Finally, we mention that we cannot conclude that 0 is an essential singularity, if we instead of z” have 
chosen 

1 

z„ = . = —> (J for n —» Too, 

yn 2 T *7r 
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where the square root is any one of the two possible definitions. In this case we get 
f(z n ) = exp (?r 2 + in) = — exp(n 2 ) —> —oo for n —> +oo. 

Then 

f(z' n ) = exp(?r 2 ) —► Too and /(z„) = — exp(?r 2 ) —► — oo for n —> Too. 

This is, however, not sufficient to conclude that zq = 0 is an essential singularity (what it is!), because 
the two real infinities, — oo and Too, in the complex plane both are identified as oo, so we have not 
by this unfortunate choice excluded the possibility of a pole at 0. Therefore, we need the sequence 
(z"), where the limit of (/(z")) is finite. 

It was proved in Ventus: Analytic Functions Theory a-1 that exp(C) = C\ {0}. Hence, 
f(D(0, g)) C C \ {0} 

for every deleted disc D(0,g). However, according to Picard’s theorem the image has at most one 
exception point, so we conclude that 

f(D( 0, q)) = C \ {0} for every g > 0. 0 

3.5 Infinity oo as an isolated boundary point 

We shall in the following sections also need to consider the case, where oo is an isolated boundary 
point. This case is just as easy as the finite case in Section 3.4. There is, however, a psychological 
obstacle here, because what may be obvious in the finite case if often difficult in the infinite case, and 
vice versa. We shall later benefit from the results in this section, and the gain is indeed much bigger 
than the effort we must use here to understand, what is going on “around oo”. 

Let / : fl —> C be analytic in an open domain H, which has oo as an isolated boundary point. This 
means more precisely that there exists an a > 0, such that C \ B[0, a] C H, so every z £ C, for which 
|z| > a, belongs to fb 



Figure 13: A deleted neighbourhood of oo is given by C \ B[0, a}. 


89 


Download free eBooks at bookboon.com 





Calculus of Residua 


Laurent Series and Residua 


Let r > a, and let C{r) denote the circle {z £ C | \z\ = r}. Then by Laurent’s theorem, 

+OO 

(114) f(z ) = ^ a nZ n for M > 
where 


n ——oo 


27T?’ 


C( r ) 


m 

?n+l 


d z, for r > a. 


If we define an analytic function g by 
9(w) :=/(—'] for — G 12, 


then we can use the discussion of Section 3.4. In fact, the domain 
deleted disc 


€ 12 


w 


of g{w) contains the 


D 



B (o,i)\(0}. 
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It follows from (114) that g(w) in D ( 0, — ) is given by the convergent Laurent series 


+ OO 


n— — oo 


3 m= a M-f = S a - 


n +oo 


n ——oo 


for |u>| < —, 
a 


and z = oo corresponds to w = 0. 

We consider as in Section 3.4 three different cases: 

3.5.1 Case I*, where a n = 0 for all positive n. 

In this case we get for \z\ > a, 

0 +oo f 1 1 n +°° 1 

f(z) = ^2 an = a ~ n 1 “ [ =^2 a -nW n , where w = 

n=— oo n —0 ^ * n —0 

When we extend f(z) to Q U {00} by adding the value /(oo) = ao, it follows that / is continuous in 

U {00}, because z = 00 corresponds to w = 0, where g(w) = f( — ) is analytic. 

\wj 

In this case it is natural to say that the extended function is analytic in the set f l U {00}, and we say 
that 00 is a removable singularity of f(z). 

When we consider the 3-transform, we get precisely series of this type. 


Example 3.5.1 It will be proved in Ventus: Complex Functions Theory a~4 that the function 

z 

f(z) := -- is the 3-transform of the constant function ip(t) = 1. It is of course analytic for 

zeC\{i}. 

We get for \z\ > 1 the Laurent series expansion, 

. +00 1 +00 

f(z) = ^ = V - = T,z~ n . 

J w z- 1 1 ^ z n ^ 

1 — — n=0 n—0 


The extension is given by /(oo) = 1, corresponding to 

z 


lim 

z^oo Z — 1 


= 1. 0 


3.5.2 Case II*, where a n 0 for finitely many positive n. 

If a n = 0 for all n > q (> 0) and a q ^ 0, then we get in C \ fl[0, a] that 

f(z) = a q z q -\ -1- ai z + a 0 + —- H-, 

z 

and it follows that f(z) —> 00 for z —» 00. More precisely, z~ q f(z) —> a q for z —> 00 . 

In this case we say that the analytic function f(z) has a pole of order q at 00, and we put /(oo) = 00. 
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Example 3.5.2 Every polynomial of degree n 
P n (z) = tin P 1 + ' ' ’ + 0,1 Z + do, Z £ C, 
where a n ^ 0, has a pole of order n at oo, and we put P„(oo) = oo. 0 


3.5.3 Case III*, where a n ^ 0 for infinitely many positive n. 

In this case we say that / has an essential singularity at oo. 

We clearly have an analogous result as in Section 3.4: 

Theorem 3.5.1 Picard’s theorem. Let f : LI —> <C be analytic with an essential singularity at oo. For 
every C \ U[0, a] C Ll the image /(C \ B[ 0, a]) is either C or C \ {w} with precisely one exceptional 
point w £ C. 


Example 3.5.3 Important. Every “nice” transcendent function like e.g. 

exp z, sin z, cos z, sinh z, cosh z, 
have all an essential singularity at oo. 

It follows in particular from Picard’s theorem that none of them has a well-defined limit for z —> oo 
in C. They may, however, have well-defined limits on certain (one-dimensional) curves extending to 
oo. For instance, expz —» 0 for z = x — oo along R_. 

On the other hand, Log: Ll — > C, where LI = C \ {R_ U {0}} does not have an essential singularity 
at oo. The reason is that the domain LI does not contain the complement of a disc, i.e. C \ B[ 0, a] of 
centre 0, so the image “Log(C \ B[ 0, a])” is not defined. 0 


Example 3 . 5.4 The function f(z) = - is defined on the oo-connected domain 

sin 0 

LI = C \ {p7T | p £ Z}. 

The denominator sin z has simple zeros for z = pn, p £ Z, because 
lim — sin z = lim cosz: = (—l) p ^ 0. 

z—>pir dz z >pTT 

Hence, f(z) has simple poles at the same points. 

7r 

If we put z' n = — + 2mr —> oo for n — > +00, we get 

/CO = 1. 

7T 

If instead = — — + 2 nn —> 00 for n —> +00, then 

/(O - -1. 

Intuitively we would then say that 00 is an essential singularity. This is not true , because LI does not 
contain any complement of a disc of centre, hence 00 is not an isolated boundary point. 0 
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We may repair Example 3.5.4 by introducing the following definition. 

Definition 3 . 5.1 Let f : SI —> C be an analytic function in an open domain, and let 
Zo£C*:=Cll {oo}. 

If there is a sequence (z n ) of either poles or essential singularities for f, such that z n —> zq for 
n —> +oo, then we say that Zo is a non-isolated essential singularity. 

Example 3 . 5.5 It follows readily from Definition 3.5.1 that 
tan z, cot z, tanli z, cothz, 
all have a non-isolated essential singularity at z$ = oo, and that the function 
1 

r 

sin - 

has a non-isolated essential singularity at zq = 0. 0 

From time to time we shall meet problems containing non-isolated essential singularities. This is the 
reason why we have given them a name in Definition 3.5.1, so we can identify them. Then in practical 
computations the rule of thumb is always to avoid this type of singularity, and instead find another 
method than just considering such a nasty singularity. We shall in the following only consider isolated 
singularities, i.e. removable singularities, poles or essential singularities. 


3.6 Residua 

Using the previous sections on isolated singularities we shall now introduce the important concept 
of a residuum of an analytic function at such points. The powerful applications of this theory is 
for practical reasons postponed to Chapter 4. In the remainder of this chapter we shall define the 
residuum and derive some easy rules of computations of it. 

Let / : 12 —> C be analytic in an open domain 12, and let zq G C be an isolated boundary point of fl, 
i.e. a singularity of /. Let C be any simple closed curve in 12 surrounding zq and no other singularity 
of /. Then it follows from Cauchy’s integral theorem that the value of the integral j> c f(z)dz is the 
same for all such closed curves around Zq. This shows that the following definition makes sense. 

Definition 3 . 6.1 The residuum of the complex differential form f(z) dz at zq is defined as 
(115) res (f(z) dz-, z 0 ) := j) f(z) dz, 

where C is any simple closed curve in 12 surrounding zq an d no other boundary point of 12. 

In general, the correct notation res (f(z) dz; zq) is too clumsy, so we shall usually incorrectly write 
res (/; zq) instead, when there is no risk of misunderstanding. One may in the literature also find other 
notations, like e.g. Res[/;zo] and res/(zo). It should, however, be emphasized that the residuum is 
a number, which is connected with the differential form f{z) dz and not the analytic function f(z) 
itself. We shall later explain this point in more details. 
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Figure 14: A simple closed curve C in surrounding just one isolated singularity zq of f(z). 


Remark 3.6.1 “Residuum”, or in some texts “residue”, is a Latin word meaning “remaining”, i.e. 
what is left after a part is taken away, namely the often complicated process of integrating the 
differential form f(z) dz along the simple closed curve C around zq. We shall in the present text use 
the Latin plural “residua”, but one may of course also use the English plural “residues”. 0 
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The first result is very simple. 

Theorem 3.6.1 Let f(z) have the convergent Laurent series expansion 

+ OO 

(116 ) f{z)= ^2 a n (z — Zq) u , forO < \z - z 0 \ < g, 

n =—oo 

in a deleted disc D(zo , g) = B(zq, g) \ {-o}- Then 

(117) res(f;zo) = a-!. 

Proof. The Laurent series (116) is uniformly convergent on the circle C (zo,r) : fz — Zo\ = r for all 
fixed r £ ]0, g[. Thus, it can be integrated termwise, 


res 


+oo 


(/: Z 0 ) = <f f(z) dz= a n-7Z~ 

2m JC(z „,r) nt^oo J C{z 0 ,r) 


(z - Z 0 ) n dz = a- 1, 


where we have used the result 


0 

(z — zo) n dz = (f> z n dz = 

C(z 0 ,r) Jc( 0,r) | 27 xi 


proved in Ventus: Complex Functions Theory a-1. □ 


for n£Z \ {—1}, 
for n = — 1, 


Theorem 3.6.1 above is in particular applied when we shall find the residuum at an essential singularity. 
It may of course be used in general, but it is often easier to apply other methods, when the singularity 
is a pole. 


Example 3.6.1 We shall find 

“(^j 

The Laurent series expansion from z$ = 0 is in the deleted disc 0 < \z\ < 1 given by 

+oo 


1 


1 1 


Z 2 (Z - 1) Z 2 1 - 3 

thus a_i = —1, and we get 

1 


_ * t ;2 „2 * 


■ • .z — 


n=0 


res 


o = -i. 


^ 2 (*-i)’ 

As a consequence - cf. (115) - we see that also 


■dz = 27 n res ( ——; 0 I = — 27ri. 


I*H 


2 (^-l) 


z 2 (^l)’ 


Clearly, the traditional computation of the line integral in the left hand side of this equation becomes 
very difficult, so we have indeed derived an easier method of computation in this case. 0 
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Example 3.6.2 It is obvious that the Laurent series expansion of 
centre 0 only contains terms of even exponents , so a_i = 0, and 
1 


^ 2 - 1) 


in a deleted disc of 


res 


z 2 (z 2 -l) : 


; o =o, 


in which case it is not even necessary to find the explicit Laurent series expansion. We just use a 
trivial inspection. 0 


An important observation is that if zo is a removable singularity , then trivially a_i = 0, so 
res(/; Zq) = 0, if Zo is a removable singularity. 

This implies that if we occasionally include removable singularities in the Calculus of Residua , this 
does not matter much because such removable singularities will not contribute to the final result. This 
is convenient because it in general in many cases suffices to notice that zq is a pole of at most some 
order q. 



Let C be a simple closed curve in the domain f 1 of /, such that / is analytic inside C, with the 
exception of a finite number of isolated boundary points z i, ..., Zk of Cl, cf. Figure 15, where k = 2. 

“Inside” means here the bounded set surrounded by C, and the direction of C is in the positive sense 
of the complex plane, which also means that the the bounded component lies to the left of C seen in 
the direction of C. 

It follows from Cauchy’s integral theorem that 

I f(z)dz= (j) f(z)dz+--- + <f> f(z)dz, 

J C J Ci J Ck 

where each simple closed curve Cj (typically a small circle of centre Zj) only surrounds one singularity 
Zj, j = 1,, k. Using that 


P f{z)dz 
Cj 


2ni ■ res (f (z)-, Zj), 
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we immediately get the following main theorem: 


Theorem 3.6.2 Cauchy’s residuum theorem. Assume that /(z) is analytic in an open domain Cl, 
and let C be a simple closed curve in Cl oriented in the positive sense of the complex plane and with 
only a finite number of isolated boundary points Z\,...,Zk of Cl inside C (i.e. to the left of the curve), 
and analytic at all other points inside C. Then 

1 f 

(118) — f{ z ) dz = res(f\z i)H- res(f;z k ) = Y^res(f;Zj). 

The importance of this main theorem will be made clear in Chapter 4. Before this chapter we shall 
in the next section derive some simple rules of computation of the residuum of a function / at a pole. 


3.7 Simple rules of computation of the residuum at a (finite) pole 

It will in the following be convenient to consider a removable singularity as a pole of order 0. We shall 
in this section only consider residua of finite poles of order g £ No- 

If zo is an essential singularity of /, we either apply Theorem 3.6.1, or a technique which we shall 
develop in Section 3.8. 


Theorem 3.7.1 Assume that /(z) has the pole zq of order < q for some q £ N. Then 
1 < 7 9-1 

(119) res (/; z 0 ) = hrn {(z - z 0 ) q f(z)} . 


Proof. We have assumed that the order of the pole is at most q, hence a_ n = 0 for all n > q, and 
we have the Laurent series expansion, 


m = 


a-i 


+ cio + ai (z — Zq) + ' 


for 0 < \z - zq\ < q. 


(z - z 0 f z- z 0 

Notice that we do not assume that a _ 9 ^ 0. 

Multiply this equation by (z — Zo) q to get 

(z - z 0 ) q f{z) = a_ q H-ho _ 2 ( 2 - z 0 ) q ~ 2 + a_i (z - z 0 ) 9_1 + (z - z 0 ) q {a 0 H-} . 

By q—1 successive differentiations the polynomial of degree q — 2 disappears, so 

r i( z - z o) q f{z)} = (q~ 1 )! 0-1 + {z - z 0 ) {•■■}, 

where the dots are a shorthand for some analytic function. Then divide by (q— 1)! and take the limit 
z —> Zo to get 

1 d 9-1 

lim -—— {(z - z 0 ) 9 /(-)} = o_i + 0 = o_i = res (/; z 0 ) , 


(q — 1 )! 2 —> z o dz q 
and the theorem is proved. □ 
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Remark 3.7.1 The proof can of course be generalized to find all coefficients a-j, j = 1, ..., 1 (and 
of course also of the not so interesting coefficients a n , n £ No). By a modification of the proof above 
the reader easily verify that 


( 120 ) a-j = 


lim 


(q — j)\ z^z 0 dz q i 


d ^ {(z - z 0 ) q f(z)} . 


j = l,...,q. 


When f(z) is a rational function of multiple poles, this is just decomposition in a new way, because 


a-j is precisely the coefficient of {z — zq) 

0 


-1 


in the decomposition. The details are left to the reader. 


An important special case of Theorem 3.7.1 is the following: 

Theorem 3.7.2 Assume that zq is either a simple pole or a removable singularity of f(z). Then 
(121) res(f-,z 0 )= lim (z-z 0 )f(z). 

Z^Zo 

Proof. Put q = 1 into (119). □ 



AACSB 


ACCREDITED 


Excellent Economics and Business programmes at: 

"V &r 


university of 
groningen 




www.rug.nl/feb/education 


“The perfect start 
of a successful, 
international career.” 


CLICK HERE 

to discover why both socially 
and academically the University 
of Groningen is one ofthe best 
places for a student to be 


98 




Click on the ad to read more 


Download free eBooks at bookboon.com 






Calculus of Residua 


Laurent Series and Residua 


Example 3.7.1 Theorem 3.7.1 has the convenient property that one shall not find the exact order 
of the pole. Let us e.g. consider the function 


m = 


sin 2 2 


for 2 G C \ {0}. 


Since 


sin z 


has a removable singularity at zq = 0, and 


m = 


sin 2 z 


sin z 


we conclude that f(z) has a pole of exact order 3 at Zq = 0. However, if we put q = 3 into (119), we 
get 


res 


sin 2 z 


0|=± lim 

2! z— >o dz 2 


sin 2 2 


which clearly will give us some difficulties, if we continue this computation. 


We shall not do this, for if we instead use that zo = 0 is at most of order q = 5 (> 3), then we get 
from (119), 


res 


sin 2 z 



— lim ——r jsin2z} = — lim 2 3 {— cos 22;) 
24 z—>o dz 3 1 ; 24 z^o 1 J 


1 

3 


In other words, by choosing a higher order than the exact one for the pole we ease the computations, 
at least in this particular case. 0 


Example 3.7.2 We shall find the value of the line integral along the circle \z\ = 2, 
dz. 


|z|— 2 | z(z-l) 2 

The mixture of an exponential and a rational function in the integrand will make the usual method 
of inserting a parametric description of the integration curve very complicated, if successful at all. 
Instead we notice that we have inside |^| = 2 two isolated singularities, namely the simple pole z = 0 
and the double pole z = 1, i.e. of order 2. Then by Cauchy’s residuum theorem, 

dz = 27ri {res(/; 0) + res(/; 1)}. 


'|*|=2 z ( z - !) 2 

Since z = 0 is a simple pole , it follows from Theorem 3.7.2 that 

e z 


res(/; 0) = lim 2 f{z) = lim 

z—>0 z—>0 [z — l) z 


= 1 . 


Since z = 1 is a pole of order q = 2, it follows from Theorem 3.7.1 that 

res(/; 1} = (2^1)! - 1)2/(z) ^ = H Tz 

Finally, we get by insertion 
f\z \=2 z{z - l) 2 


= lim (z — 1) = 0. 
z — * 1 z 2 


dz = 27 xi {1 + 0} = 27ri. 0 
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Theorem 3.7.3 Assume that both A(z) and B(z) are analytic in a neighbourhood of zq, and assume 
that zq is a zero of order 1 for B(z), i.e. B(zq) = 0 and B' (zo) ^ 0. Then 


(122) res 



A{z 0 ) 
B' (*„)' 


A(z) 

Proof. Since ^ has a pole of at most order 1 at z 0 , it follows from Theorem 3.7.2 that 


res 


A(z) 

B{z) 


;zo 


= .‘i” - 2 »> ■ W) = 15. A{s) ' Bgi (*,) 


z — Z 0 


— A{zq) 


lim. 


1 _ = A(z 0 ) 

B(z) - B(z 0 ) B' (z 0 ) ' 
z- Zo 


□ 


Example 3.7.3 We shall find the value of 


1 * 1=2 


z e 
z 2 — 1 


dz. 


Put A(z) = ze z and B(z) = z 2 — 1. Then A{z) and B(z) are analytic in C, and B(z) has its simple 
zeros at z = ±1. Let zo be any of the zeros ±1. By Theorem 3.7.3, 

i , \ A(z 0 ) z 0 exp(z 0 ) 1 

™ S(/; 2o) = WW) = 2 sq = 2 exp(2 »>' 


hence 


1 * 1=2 


ze 

z 2 — 1 


dz = 2iri {res(/; 1) + res(/; —1)} = 2ni ■ 


e 1 + e 1 


= 2ni ■ cosh 1. 0 


Example 3.7.4 Theorem 3.7.3 is in particular applied when we shall find the residua at several simple 
poles. We have e.g. 

dz = 2m {res(/; 1) + res(/; -1) + res(/; i) + res(/; -*)}, 


1 * 1=2 24 - 1 


where all poles , 1, -1, i and —i, are simple. Let zq be any one of these. Then Zq = 1. 
Choose A(z) = z and B(z ) = z 4 — 1, i.e. B’(z ) = 4z 3 , so it follows from Theorem 3.7.3 that 


res(/; zo) = jzj 


A(z 0 ) zo 1 zl 1 


B'(z 0 ) 4zg A'4 4 Z °’ 


hence by insertion, 


1 * 1=2 


? ^ T d 2 = 2 p{C + (-i) 2 + i ! + H)»} = o. 0 
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In the following theorem we assume that B(z) has a zero of second order at zq. It is given here for 
completeness, but it should also be mentioned that Theorem 3.7.1 usually give smaller computations, 
and furthermore, it is very easy to make errors in the computations of Theorem 3.7.4 below, so the 
reader is warned against applying it uncritically. 


Theorem 3.7.4 Assume that A(z) and B{z) are analytic in a neighbourhood of zq . Furthermore, 
assume that B(z) has a zero of exactly second order at zq. Then 

oo\ (Mz) \ GA'(zo)B"(zo)-2A(zo)bW(z 0 ) 

J 3{£"(z 0 )} 2 

Proof. It follows from the assumptions that 

A(z) = a 0 + ai{z - z 0 ) -\ -, B(z) = b 2 (z - z 0 j 2 + b 3 (z - z 0 ) 3 H- for \z - z 0 \ < g, 

where 62 7 ^ 0, and B(z) 7 ^ 0 for 0 < \z — zo| < Q- We shall find the Laurent series expansion of f(z) 
in the deleted disc 0 < \z — zq\ < g. This is given by 

_ c _2 c-i _ A(z) _ 1 a 0 + qi (z - 00 ) -I- 

(Z — Zo ) 2 z — Zo B(z) (z — Zo ) 2 + b 3 {z — Zo) + • • • 

Since res (/; zq) = c_ 1 , we shall only find c_i of the equation B(z)f(z) = A(z), i.e. 


{c _2 + c_i (z - z 0 ) H-} ■ {b 2 + b 3 (z - z 0 ) H-} = a 0 + ai (z - z 0 ) H-. 

By multiplication, 

c_ 262 + (c- 2 b 3 + c-ib 2 ) (z - z 0 ) H-= a 0 + 01 (z - z 0 ) H-. 

It follows from the Identity theorem that c _2 = —, and 

b 2 

(124) res (/; z 0 ) = c_i = ^ W - c. 2 b 3 } = - ja, - —| ^-. 

Using the notation of the theorem we get 

A(zq) = a 0 , A'(zo) = ai, B''(z 0 ) = 2\b 2 = 2b 2 , U ( 3 ) (z 0 ) = 3!6 3 = 663 . 

We finally get by insertion into (124), 

,, , _ A'(zq) • | B"(z 0 ) - A(z 0 ) ■ I ff ( 3 ) (z 0 ) _ 6 A'(zo) B"{zo) - 2A{zo) B^{z 0 ) 

ieS[t ' Z0} \{B"{zo)} 2 3 {B"{zo)} 2 
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3.8 The residuum at oo 

Assume that the analytic function f(z) for \z\ > R is given by its convergent Laurent series expansion 

+ OO 

(125) f(z) — ^ a n z n for \z\ > R. 

n=—o o 

Definition 3.8.1 Assume that f(z) is analytic for \z\ > R. Then we define the residuum of the 
differential form f(z) dz by 

(126) res(f(z) dz-, oo) := (j) f(z) dz, 

J c 

where —C is any simple closed curve in z > R surrounding oo. 


We notice the minus sign on the right hand side of (126), and that — C denotes the curve C with 
reversed direction, i.e. opposite the orientation of the complex plane. That —C is surrounding oo 
means that oo lies to the left of this curve —C, seen in this reversed direction of C. 

We have earlier proved that the value of (126) is independent of the curve C , as long as it fulfils the 
conditions of Definition 3.8.1. 
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Choose in particular C — (7(0, r), i.e. |z| = r, where r > R. Since we have uniform convergence of the 
Laurent series expansion (125) on (7(0, r), we can interchange summation and integration, and we get 


res (f(z) dz; oo) 


1 

27 xi 


/ 

JC{ 0,r) 


f(z)dz 


+oo 


0(0,r) 


]a n z n dz = —a_i, 


so we have proved 


Theorem 3.8.1 If /(z) = Y2n°°-oo a « zn ^ le convergent Laurent series expansion for \z\ > R, then 

the residuum of f(z) dz at infinity is given by 

(127) res(f(z) dz; oo) = — a_i. 

It is easy to extend Cauchy’s residuum theorem to also include oo. 


Theorem 3.8.2 Cauchy’s residuum theorem for unbounded domains. Assume that /(z) is analytic 
in an open domain 12, which contains the set C\B[0,i2] for some R > 0. Let C be any simply closed 
curve in 12, such that outside C, i.e. to the right of C seen in its positive direction - which can also 
be described as the unbounded domain of C having C as its boundary - there are only a finite number 
of (necessarily isolated) boundary points Zi,..., Zk of 12. Then 

(128) f( z ) dz = res(f dz; z-f) H-h res(f dz; z k ) + res{fdz;oo). 



Figure 16: Cauchy’s residue theorem for unbounded domains. 


Proof. Choose r > 0, such that \zj\ < r for all j = 1,..., k, so all the finite singularities Zi,..., z k lie 
between C and (7(0, r). We get by adding and subtracting the line integral along (7(0, r) and applying 
Cauchy’s residuum theorem for finite singularities and Definition 3.8.1, cf. Figure 16, 


1 

27t i 


f(z)dz 


1 

27rz 


f(z)dz + 


1 

27 xi 


/ 

JC( 0,r) 


f(z)dz 


1 

2m 


/ 

JC( 0,r) 


f(z)dz 


= res(/dz; Zi) + • • • + res(/dz; Zfc) + res(/dz; oo). □ 


103 


Download free eBooks at bookboon.com 





Calculus of Residua 


Laurent Series and Residua 


Theorem 3.8.3 Assume that f : Q, —> C is analytic, where = C \ {zi,..., Zk}- Then the sum of all 
residua, including the residuum at oo, is equal to zero, 

(129) res(f d z; Zi) + • • • + res{f dz; Zk) + res(f dz; oo) = 0. 

Proof. Let C be any closed curve in f l. Then by the two versions of Cauchy’s residuum theorem, 

0 =--*- & f(z)dz+ (f /(z)dz = res(/dz;zi) 4-f- res(/dz; z k ) + res(/dz; oo). □ 

2m J c 2m J c 


Theorem 3.8.4 If f(z) has a zero at oo, then 
(130) res(fdz-oo) = — lim z f(z). 

z—>o o 

In particular, res(f cl z\ oo) = 0, if oo is a zero of order > 2 for f. 

PROOF. By the assumption, 

f(z) = —- + H- for | z | > R. 

z z 

Then by Theorem 3.8.1, 

res(/dz;oc) = — a_i = — lim zf(z). □ 


The residuum at oo is in particular applied for line integrals along simple closed curves C, for which 
the residuum theorem of Section 3.6 is either difficult to apply, or where its assumptions are not 
fulfilled at all. We shall in the following give some examples of these phenomena. Notice that we often 
can choose between various methods of computation. The art is then to choose the easiest one. 


Example 3.8.1 We computed in Example 3.7.4 the integral 

z 


W=2 


z 4 - 1 


dz 


by applying Theorem 3.7.3. It is, however, much easier to apply Theorem 3.8.4, 


1 * 1=2 ^ - 1 


dz = — < — 


1 * 1=2 


Z 4 - 1 


dz > = —2m ■ res 


z 4 - 1 


dz; oo I = 2ni ■ lim 


z 4 - 1 


= 0 , 


because all (finite) singularities of the integrand lie inside the circle \z\ = 2, and because the integrand 
has a zero of order 2 at oo. 0 
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Theorem 3.8.5 If f(z) is analytic for |z| > R, then 


(131) res(f(z) dz; oo) = -resl \ f(- ) dz;0 ) . 


Proof. The Laurent series expansion is 

+OO 

f( z ) = E for N > R ' 


SO 


+oo 


res (i/(c) d ^ ;0 ) = res ( E 2 


= a_i = —res(/dz; oo). □ 


\Tl— — OO 


Remark 3.8.1 Notice that 

-res (? / G) d * ; °) =res ( / (l) d (l) ; ^=°) ■ 

so (131) can be rewritten as 

(132) res(/(z) dz; oo) = res d 

which indicates why the residuum is linked to the differential form and not to the function /(z) itself. 
In fact, in this way the residuum becomes invariant under transformations of z. 0 



Example 3.8.2 Important! The analytic function /(z) = -, z yf 0, is trivially extended to oo by 

z 

putting /(oo) = 0. The Laurent series expansion of / is trivial, 

+oo 1 1 

f( z )= V a n z n =- = a _i • 

' z z 

n =—oo 

thus a_i = 1, and res(/dz; oo) = — 1 0 according to Theorem 3.8.1. This simple example shows that 
analyticity at oo does not imply that res(f dz\ oo) is zero! This is in contrast to the residuum at a finite 

removable singularity, in which the residuum is indeed zero. Many years ago the author experienced 

that even trained professors in Mathematics can make errors here. Therefore this warning. 0 
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Example 3.8.3 To show the power of this theory we shall show that we are now able to compute 
even a nasty integral as 



The standard procedure of computing a line integral, using a parametric description is clearly doomed 

to failure. The integrand has its singularities at the zeros of sin -, i.e. for z = — for n £ Z \ {0} (all 

z mr 

simple poles lying inside \z\ = 1), supplied with the non-isolated singularity at z = 0, so we cannot 
use the version of Cauchy’s residuum theorem given in Section 3.6. 


We then choose to apply the residuum at oo. This is done in the following way, 


1 

|*|=i sin % 


dz = - { - 


1 

|*|=i sin % 


dz > = —27ri ■ res 


1 


sin - 


dz; oo 


= — 2ni • res 


1 J 1 


sin z 


cl - ; 0 = 27ri • res 


dz; 0 I = 27ri • a_i, 
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where o_i is the coefficient of - in the Laurent series expansion of „- from zo = 0. Now, 

z z 2 sin z 

z 2 sin ^ clearly has a zero of order 3 at 0, so ;— has a pole of order 3 at 0. Using that z 2 sin x ^ 0 

z 2 sin z 

for 0 < \z\ < 7r we therefore have 


1 


a_3 a_2 , a _1 


for 0 < \z\ < 1 r, 


where the task is to find a_i. If this equation is multiplied by 


2 • 2 I x 3 

2 smz = z iz ~e z 




then we get 

1 O o ( z 2 

1 = — 3 {u -3 + CI— 2 Z + a-\Z 2 + ■••}• z <1 —— + • ■ 

= {u -3 + CI- 2 Z + a_i z 2 + •••}• ^1 — — z 2 + • • • 


= a-3 + a- 2 Z + ^ a_i — - a_ 3 f z 

where the dots everywhere indicate terms of higher order. When we identify the coefficients, we get 

a _3 = 1, a -2 = 0 and a_i = \ a_ 3 = 

6 6 


Hence by insertion, 
f 1 
' 1 * 1=1 ~ 


1 \ 7TZ 

dz = 2ni ■ res I -dz; 0 I = 2-jri • a_i = —. 0 

1 sin z V z 2 sin z 1 3 


Example 3.8.4 Since f(z) = sin - is analytic in C \ {0}, it follows from Theorem 3.8.3 that 

z 

res^sin-;oo^ = — res^sin -; 0 ^ = — 1 , 

because the Laurent series expansion of f(z) from z : 0 = 0 is given by 

f(z) = sin 1 = 1 - ^ ^ H-, for z e C \ {0}, 

Z Z 01 Z° 

from which a_i = 1 . 

Alternatively we apply Theorem 3.8.5, because 
1 c( l \ 1 • 1 / ^ 

/ - 1 = -ffSinz=-^<z- — 


thus 

res ( sin 00 ) = — res( \ f(~ ) ; 0 ) = —1. 0 
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Remark 3.8.2 It must be emphasized that oo is not 


-, —, —-—, tana;, cot z, tanhz, cothz, etc.. 

cos z smh z cosh z 
defined for these function. 0 


an isolated singularity of functions like -, 

Sill z 

This means that the residuum at oo is never 


3.9 Summary of the Calculus of Residua 

It is very important to be able to compute the residuum in Chapter 4, so we collect the basic properties 
in this section. 

1) The residuum of the complex differential form f(z) d z at Zq is defined as 

res (f(z) dz; z 0 ) := <j> f{z) dz, 

where C is any simple closed curve in f 1 surrounding zq and no other boundary point of Cl, cf. 
Definition 3.6.1. 

2) Let f(z) have the convergent Laurent series expansion 

+oo 

f(z)= ^2 a n (z-z 0 ) n , for 0 < \z - z 0 \ < g, 

n=— oo 

in a deleted disc D(z 0 ,g) = B(z 0 , g) \ {ar 0 }- Then 
res (/; z 0 ) = a_i. 

Cf. Theorem 3.6.1. 

3) If zq is a removable singularity, then 

res(/; z 0 ) = 0. 

4) Cauchy’s residuum theorem. Assume that f(z) is analytic in an open domain Cl, and let C be a 
simple closed curve in Cl oriented in the positive sense of the complex plane and with only a finite 
number of isolated boundary points zi,...,z k of Cl inside C (i.e. to the left of the curve), and 
analytic at all other points inside C. Then 

2 — j c f ( z ) dz = res (/; z i)d -b res (/; z k ) = res (/; Zj), 

cf. page 93. 

5) Assume that f(z) has the pole zq of order < q for some q £ N. Then 

i( z ~ z °y /(«)}» 

cf. Theorem 3.7.1. 


res (/; z 0 ) = 


(g- 1)! z 


lim 

20 
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6 ) Assume that zq is either a simple pole or a removable singularity of f(z). Then 


res (/; z 0 ) = lim (z - z 0 ) f(z), 

Z^-Zo 

cf. Theorem 3.7.2 

7) Assume that both A(z) and B(z) are analytic in a neighbourhood of zo, and assume that zq is a 
zero of order 1 for B(z), i.e. B(zq) = 0 and B' (c 0 ) if 0. Then 


res 




',zo = 


B{z) 

cf. Theorem 3.7.3 


A{z 0 ) 
B' {z 0 ) ’ 


8 ) (Use this result with care.) Assume that A(z) and B(z) are analytic in a neighbourhood of zq. 
Furthermore, assume that B{z) has a zero of exactly second order at zq- Then 


res 




;zo = 


B(z) 
cf. Theorem 3.7.4. 


&A' (z 0 ) B" {z 0 ) - 2A(^ 0 ) B& (z 0 ) 
3{B"(z 0 )} 2 


9) Assume that f(z) is analytic for \z\ > R. Then we define the residuum of the differential form 
f(z) d 2 by 


res (f(z) dz-, oo) := j> f(z) dz, 


where — C is any simple closed curve in z > R surrounding oo, cf. Definition 3.8.1 

10) If f(z) = Y^n°°-oa a n z n is the convergent Laurent series expansion for \z\ > R, then the residuum 
of f(z) dz at infinity is given by 

res(f(z)dz-, oo) = -o_i, 
c.f. Theorem 3.8.1. 

11) Assume that f(z) is analytic in an open domain U, which contains the set C \ B[0,R] for some 
R > 0. Let C be any simply closed curve in U, such that to the right of C seen in its positive 
direction there are only a finite number of (necessarily isolated) boundary points z\,... ,Zk of Q. 
Then 


f{z) dz = res(/ dz; zi) H-h res(/dz; z k ) + res(f dz; oo), 

cf. Theorem 3.8.2. 

12) Assume that / : Q —> C is analytic, where f2 = C \ {^i,..., z k }- Then the sum of all residua, 
including the residuum at oo, is equal to zero, 

res(/d^; zi) H-+ res(/dz; z k ) + res(f dz; oo) = 0, 

cf. Theorem 3.8.3. 
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13) If f(z) has a zero at oo, then 

res(/cb;oo) = — lim zf(z). 

z —>-oo 

In particular, res(/clz; oo) = 0, if oo is a zero of order > 2 for /, cf. Theorem 3.8.4. 

14) If f(z) is analytic for \z\ > R, then 

les (f( z ) dz; oo) = -res^ /(^) d^;0^ , 
cf. Theorem 3.8.5. 
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4 Applications of the Calculus of Residua 


4.1 Trigonometric integrals 

We shall start this chapter by demonstrating that some trigonometric integrals are easier to compute 
by using the Calculus of Residua than the traditional method from Real Calculus. We first prove the 
following theorem: 

Theorem 4.1.1 Let R(£,rj) denote a function in two real variables defined in a subset o/R 2 . If 


m = r 


z 2 - 1 z 2 + 1 


2 iz 2z 

is analytic in a domain O C C, containing the circle \z\ = 1, then 

' z 2 — 1 z 2 + 1 \ dz 


c 2"7T /» 

(133) / l?(sin0,cos0) d© = ® 

J o 11*1=1 


R 


2 iz 


2z 


iz 


Proof. When we apply the parametric description 2 = e* 0 , 0 G [0,27r], of \z\ = 1, we get 
d z = i e 10 d0 = iz cl0. Then 

E—^ = 1 t e ie -e" 10 } =sin0 and E±1 = I {e ie + e~ 10 } = cos©, 
and the result follows immediately by insertion into the right hand side of (133). □ 

Obviously, (133) should be applied from the left to the right, because then we can apply the Calculus 
of Residua. That this method is really powerful is demonstrated by the following example. 


Example 4.1.1 We shall compute e 2cos& d0. The tradition substitution t = 2cos0 with a 
discussion of its intervals in which it is monotone does not look promising. Instead note that -R(£, rf) = 
e 2r? , in which even £ is missing, and where the function 


R 


^+1\ 

2 iz ’ 2z J 


= exp 




= exp 




2 / 0, 


is analytic in C \ {0}. Hence, by Theorem 4.1.1, 


r 2n 


e 2 cos e d Q = 




d2 

iz 


2ni 

-res 

i 


1 

z 


exp 




The unpleasant fact is of course that the only two singularities, 0 and oo, are both essential. This 
means that the only possible method is to find a_i in the Laurent series expansion. We have 


1 

z 


exp 



1 1 

- exp z • exp - 
z z 


1 

z 


+oo 

E 


m =0 


l 

m! 


+oo 

•E 


71=0 


1 
n ! 



*ec\{0}, 


so a_i is by Cauchy multiplication equal to the sum of the coefficients, which correspond to m = n 

-v+OO 1 


i.e. by a summation, a_i = X),i=o 


("!) 


2 ' 


Hence, 



+oo 

e 2cose d0 = 27r ^- 
71=0 


1 

P) 2 ' 
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This series is clearly rapidly convergent. It should here be mentioned that it can be proved that 

+oo 1 

X] t n \\2 = 

n=0 V 

where Jq{z) denotes the (complex) Bessel function of order 0. 0 


The value of Theorem 4.1.1 lies in the fact that it can be applied when the integrand is not a rational 
function in cos© and sin©. In many textbooks this theorem is however only formulated for such 
rational functions. In some sense this is an overkill, because the traditional method known from Real 
Calculus is often easier to apply. We shall demonstrate this in Example 4.1.2 below. 


Example 4.1.2 We shall compute the integral J^ n ^ " ■ q ’ ^irst we by Theorem 4.1.1, 


r 2-K 


d© 

2 + cos 0 


1 


dz 


-2 i 


UM 2 , * 2 + 1 ™ T\z\=i z 2 +4z + l 
2 z 


dz 


= (—2 i) ■ 2ni ■ res 


1 


z 2 +dz+l 


—2 + \/3 ) = 47t lim 


2tt 


—2+\/3 z T 2 T "\/3 


where we have used that z 2 + dz + 1 has the two simple roots —2 ± \/3, of which only —2 + a/ 3 lies 
inside \z\ = 1, and then used Theorem 3.7.2. 

Alternatively , we get by the traditional substitution f = ~ and u = tant, that 


d© 


2 + cos 0 


d© 


3 cos 2 


sin 


0 


= 2-2 


4 [ + °° _du 


>o 


1 + lu 2 


4-\/3 ir 2i A 
= 0 


dt 


/o 3 cos 2 t + sin 2 t 
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4.2 Improper integrals 

Up to this point we have only considered bounded line integrals. We shall here prove that the real axis 
R. in some cases - depending on the integrand - may be considered as a simple curve which is closed 
by adding oo, where we identify +oo and —oo in the complex plane. Hence under some additional 
assumptions we shall expand Cauchy’s residuum theorem to half planes. 


Theorem 4.2.1 Let f : H —> C be analytic in an open domain 12, which contains the closed upper 
half plane ?sz > 0 , except for a finite number of points z\,... ,z n , all of which lying in the open upper 
half plane ‘isZj > 0, j = 1,..., n. Thus, 


12 U {zi, ..., z n } D {z £ C | isz > 0}. 

Assume that there exist constants R > 0, c > 0 and a > 1, such that 

(134) |/( 2 )| < -j— 7 for | 2 | > R and 32 > 0. 

\ z \ 

Then the improper integral f _^ /( x) dx along the real axis is well-defined, and its value is given by 
the following residuum formula, 



Proof. First note that R + i’Ocd and that / is analytic in 12, so the restriction of / to R must be 
continuous. From a > 1 and (134) follows that /(x) in the real has an integrable majoring function, 

e.g. 


— for |x| > R , 


g(x) = 


1/0*01 for \x\ < R, 

and we conclude that the improper integral f(x ) dx exists, and its value is given by 


/*+0O 


f(x) dx = lim 

r 1, r 2 —*+00 


f(x) dx = lim / /(x) dx. 

r * • ot /_ 


We then exploit that f(z) is analytic in the upper half plane with the exception of only a finite number 
of singularities. 


It follows from (134) that all singularities in the upper half plane must lie in the disc B(0,R). 

Let r > R, and let C' r denote the circular arc of the parametric description 2 = re l& , 0 £ [0,7r], in 
the upper half plane. Let furthermore C r be the simple closed curve which is obtained by joining C' r 
and the interval [—r, r] on the real axis. Then by Cauchy’s residuum theorem , 


(136) 


[ f{x) dx + [ f(z) dz = (f f(z)dz = 2iri V res (f;Zj). 
J-r JCt . Jc r Sj2,->0 
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Figure 17: The path of integration in the proof of Theorem 4.2.1. 


The right hand side of (136) is constant for all r > R. The first integral on the left hand side of (136) 
converges towards J f[x ) da; for r —> +oo. Therefore, we shall only prove that 

lim [ f{z) dz = 0. 

r ^ +00 j c ^ 


Using that \z\ = r > R for 2 £ C' r of length 


f(z)dz 


< ^HC'r) = C7TT 1 


£ (C' r ) = irr, we get from (134) the estimate 
0 for r —> +oo, 


because a > 1 by assumption. □ 


Example 4.2.1 We choose arbitrarily 


/(*) = 


1 


z 2 + 1 


exp 


1 


z — i 


z € C \ {—*, i}, 


where it is more or less obvious that 


f*+oo 


X 2 + 1 


exp 


x + i 
x 2 + 1 


/» + 00 


dx = 


x 2 + 1 


exp 


x 2 + 1 


cos 


X 2 1 


i sm 


x 2 + 1 


dx 


cannot be computed by only using traditional methods from Real Calculus. 


Since 


1 


z — i 

an R > 1, such that 


0 for z —> oo, and thus exp 


1 


z — i 


e° = 1 for z —> oo, it is obvious that there exists 


\f(z)\ < —ry for \z\ > R, 


and it follows from Theorem 4.2.1 that 

r +oo 


X 2 + 1 


exp 


x +1 
x 2 + 1 


dx = 2ni- res 


z 2 + 1 


exp 


i = 2 tti res 


1 


w 2 + 2 iw 


— exp - ; 0 , 
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where we have applied the linear transform w = z — i. 

Obviously, wo = 0 is an essential singularity of —^——— exp( — V so we shall use Cauchy multipli- 

w z + 2 iw \w J 

cation in the Laurent series expansion from wq = 0 to find a_i by collecting all coefficients of —. We 

w 

get for 0 < |u>| < 2, 


1 1 


111 


_ . -. ex p — = — ■ — 

w 2i + w w 2 1 w 


1+ 2 i 


111 
r • exp - = — • - 
1 W Zl w 


+oo 


+oo 




m =0 


n=0 


l l 

w n 


Due to the factor —, the coefficient a_i is given by — times the constant term of the product of the 
w 2 i 

two series, i.e. for m = n. This gives 


1 ^ i 


9 


2 i ' n\ 1 2 

n —0 


= 2i 6XP 


= 2f ' C ° S 


• i sm - 
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Thus, 


p+oo 


X 2 + 1 


exp 


p + OO 


X 2 + 1 


X 2 + 1 


exp 


cos 


X 2 H - 1 


Z Sill 


X 2 + 1 


dx 


X 2 + 1 


dx = 27rz a_i = 7r <! cos - + z sin - 


Finally we get by splitting into the real and imaginary parts, 


+oo 


—oo 
+oo 


^I exp (^) c “(?TT) d " 


= 7r ■ cos 


1 

2’ 


x 2 + 1 


exp 


x 2 + 1 


sin 


x 2 + 1 


dx = 7r • sin -. 

2 


This example is of course only meant to demonstrate the power of the method, because it is most 
unlikely that one inn practical applications ever will need to find the exact value of these integrals. 

An alternative method of computation is the following: 


We adopt from the above the already derived formula 


r*+oo 


1 

x 2 + 1 


exp 


J —OO 

The function 

g(w) = 


f x + i\ 

Vx 2 + lj 


dx = 27 tz • res 


w 2 + 2 iw 


exp 


;° ■ 


1 


w 2 + 2 iw 


exp - = 


w J w(w + 2 i) 


exp - 
w 


is analytic in C \ {0, —2ij with only two singularities, so it follows from Theorem 3.8.3 that 


p+oo 


exp 


x 2 + 1 


dx = 2ni ■ res(g(w)\0) = 2ni{—res(g(w); —2i) — res {g{w)-, oo)}. 


Now w — —2i is a simple pole, so we apply Theorem 3.7.2, 


— res(g(w); —2 i) = — lim — exp 

z—> 2 i w 


(^) 


= 2} eXP 


Furthermore, lim^^c 

sM = -4 • 1 


, exp — = exp 0 = 1, so w = oo is a zero of order 2 for 

\w 


1 + 


2 i 


■ exp 


w 


and we get from Theorem 3.8.4 that 
—res (g(w); oo) = 0. 

Summing up we get 


x 2 + 1 


exp 


/ x + * \ 

Vx 2 + l) 


dx = 27ri{—res(g(w);—2i) — res(g(w); oo)} 


= 7t • exp 


1 1 
= 7r ^ cos - + * sin - 


and the results follow by taking the real and imaginary parts of this result. 0 
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Example 4.2.1 showed that Theorem 4.2.1 can be successfully applied even in very complicated cases. 
However, the most important case is of course when f(z) is a rational function , because the present 
method from Calculus of Residua is easier to apply than the ordinary decomposition method. We 
therefore explicitly formulate 


Corollary 4.2.1 Let f(z) 


P{z) 

Q{z) 


be a quotient of two polynomials, where the denominator Q(x) ^ 0 


for all i£l. 

If the degree of the denominator Q{z) is at least 2 larger than the degree of the numerator P{z), i.e. if 
f{z) has a zero of at least order 2 at oo, then the improper integral f_ /( x) da; exists, and its value 
is given by (135), i.e. 



f(x) da; = 2ni res(f\ Zj). 

Qzj >0 


Proof. Since degree Q(x)— degree P(z) > 2, we can choose a = 2 and then constants c, R > 0 to 
fulfil the assumptions of Theorem 134, and the corollary follows. □ 


Example 4.2.2 We shall find the value of the improper integral 
r+ °° da; 


/ 

J — C 


a; 4 + 1 


The integrand 


a; 4 + 1 


is a rational function, where the denominator a; 4 + 1 > 1 > 0 has degree 4, and 


the numerator is a constant. Hence, the assumptions of Corollary 4.2.1 are fulfilled. 

( X p 7T \ 

^ J , p = 1, 3, 5, 7, and they are all simple. Hence, the 
integrand has the same points as simple poles, of which only 

and 


exp (r) = 7i (1+i) 


f 3in\ 1 

ex TTj = ^ ( - 1+,) 


lie in the upper half plane. 

Then apply Theorem 3.7.3 with 

A(z) = 1, B(z) = z 4 + l and B\z) = 4 z 3 , 
and use that zfi = — 1 for all poles, to get 

1 

6o -a, 4 -- 


,, x A(z 0 ) 
re s(/; zo) = ^7 


1 _ z 0 _ 

B' (z 0 ) ~ 4Zg ~ 4Zg _ 4 Z °' 


Finally, by Corollary 4.2.1, 

r +co da; 

J- oo 1 + a; 4 


= 27rz< res 


1 + i 


= 27 ii < — - 


1 + z 4 ’ V2 
1] fl + i — 1 + i 


res 


1 -1 + i 


1 + z 4 


V2 


y/2 y/2 


2ni 2 i 7r 

4 ' y/2 - V2.' 
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Alternatively we sketch the following clever real decomposition , where where we first notice that 
x A + 1 = a: 4 + 2x 2 + 1 — 2x 2 = ( x 2 + l) 2 — ^%/2a/j = ^ x 2 + \plx + 1^ (x 2 — V2x + 1^ . 
There exist four uniquely determined constants A, B, C, D £ R, such that 


Ax + B 


Cx + D 


+ 1 (x 2 + y/2x + l) (x 2 — \/2x + l) x 2 + V2x + l x 2 — \/2x+l 


We get by some very tedious computations, 


A = 


1 


2^2’ 


B=~, 

2 ’ 


C = — 


from which an indefinite integral of 


2y/2 

1 

x 4 + 1 


and 


D = -, 
2 ’ 


can be found, and then the value of the improper integral 


by taking the limits. The details are left to the reader. 


A simpler alternative method is applying complex decomposition , where 

1 = res(/;zi) + res(/;z 2 ) + res(/; z 3 ) res(/; z 4 ) 

Z 4 + 1 Z — Z\ Z — Z2 z — Z3 Z — Zi ’ 

and res(/: Zj) = from the above. Then pair the results to get the real decomposition , which 

then is integrated in the usual way. 0 


Obviously, Example 4.2.2 shows that residuum formulas may be easier to apply than a straightforward 
real decomposition followed by an integration. Furthermore, Example 4.2.1 showed that we also can 
compute improper integrals, which could not be found by traditional real methods. However, the 
reader must be warned. If not all assumptions of a residuum formula are fulfilled, then it usually give 
a very wrong result, even if the residuum formula itself makes sense. We shall illustrate this by the 
following obvious example, but it is easy to give more subtle examples showing the same phenomena 
in a latent way. 


Example 4.2.3 The rational function 


2 

z 2 + 1 


does not fulfil the conditions of Corollary 4.2.1, because 


the difference between the degrees of the denominator and the numerator is only 1. Clearly, z = i is 
the only singularity - a simple pole - in the upper half plane, where by Theorem 3.7.3 


res 


z 

z 2 + 1 




1 

2 


Then a false application of (135) gives 


r+oo 


x 2 + 1 


da; = 27ri • res 


z 2 +l’ 


;t = 7T i 


which is wrong for several reasons. The indefinite integral is — In (a: 2 + l) —> +oo for x —> +oo or 

x —> — oo, so the improper integral does not exist. And if it did, its value ought to be real and 0 by 
the symmetry of the integrand, and not the complex number iri. Therefore, the reader should always 
be extremely careful to check all the assumptions of the applied theorem before using the residuum 
formula of this theorem. 0 
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We have also the following result of improper integrals, when the integrand is a product of an analytic 
function and a complex exponential. 


Theorem 4.2.2 Let f(z) be analytic in the open domain f l, which contains the closed upper half 
plane > 0 with only exception of a finite number of singularities zi,... ,z n , ll lying in the open 
upper half plane, so ?sZj > 0 for j = 1,... ,n. Assume that there exist positive constants R, a, c > 0, 
such that 

(137) \f{z)\ < for > 0 and \z\ > R. 

\z\ a 

For every real positive number m > 0 the improper integral f_°° /( x) e lmx dx on the real line is 
convergent of the value 

/ +oo 

f(x) e imx dx = 2iri ^ res(f(z)e imz -, Zj ). 

■°° Q Zj >0 


We have emphasized the important assumption that the constant m is positive, because otherwise 
(138) is not true. We shall, however, also deal with negative constants in Corollary 4.2.2 below. 



> Apply now 
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Figure 18: The path of integration in the proof of Theorem 4.2.2. 


Proof. We shall under the given assumptions prove that the limit 
lim lim / 2 /( x)e imx dx 

r l —>+oo r 2 —>+oo J_ ri 

exists and is unique, where the two limits are taken independently of each other. 

Choose r i, r 2 > R and the simple closed curve Ci, 2 (the boundary of a square) of Figure 18. 

According to (137) all singularities in the upper half plane of the integrand f(z) e lmz lie inside Ci, 2 , 
hence by Cauchy’s residuum theorem , 


P n 

(139) j f{z)e imz dz = 2iri^2ies(f{z)e imz -,z j ). 

” C*1,2 j — l 


Clearly, when n, r 2 > R, then the right hand side of (139) is independent of the choices of r\ and r 2 . 


On the other hand, cf. Figure 18, 

/ i"2 rri+r 2 

f{x)e imx dx+ / /(r 2 + it) e im ^ r2+it \dt 

-ri J 0 

/ r 2 pr 1-+T2 

f(x+i{r 1 +r 2 })e im(x+l{ri+r2 }'>dx- / f(-n+it) e im( -~ ri+it) idt. 

-n J 0 


Ci, ; 


(140) 


The positive constant m > 0 is fixed, so we get the following estimates for ri, r 2 > R, 


rri+r2 


. f(r 2 +it ) e^^idt 


rri+r2 


< 


< 


\f(r 2 +it)\ e~ mt dt < — 


C r l+T2 


dt 


mrn 


'2 Jo 

for r 2 —> Too , because a > 0, 
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and similarly for the upper horizontal line integral, 


[ 2 f(x+i{r 1 +r 2 })e lm(x+i{ri+r2}) dx < f ' \f(x+i{n+r 2 })\ e" m(ri4r2) da 
J—r^ J ri 


< e 


—m(r i-+r 2 ) 


l—n (r 1 + r 2 ) 


adx = c- (r 1 +r 2 ) 1 “° • 


Since exponentials dominate polynomials, this tends towards zero if either r\ —> +oo or r 2 —* +oo. 


Finally, 

rri+T2 


rr i-H"2 /■i'i-H'2 

/ /(-n+tt)e iro (- rrHt )*dt < / |/(-n+*i)| e 

Jo Jo 


; dt 


/-ri-p-2 

< — / e _m * di <-> 0 for n + 00 . 

Ti Jo rnrf 

When (140) is inserted into (139), we get by taking the two independent limits r\ —» +oo and r 2 —> +oo 
that the improper integral /( x) e mx dx exists with the value (138). □ 


As mentioned above we shall also consider a negative constant m. This will, however, require that 
the integrand is bounded in the lower half plane, with the exception of in the neighbourhoods of the 
finitely many singularities. 


Corollary 4.2.2 Residuum formula for the Fourier transform. Let f(z) be analytic in the open 
domain S! = C \ { zi ,..., z n }, where none of the singularities Zj lies on the real axis. Assume that 
there are positive constants R,a, c > 0, such that 

(141) \f(z)\ < for \z\ > R. 

Then 

r +oo ( 27r *E^>o res (/( z ) ei2y ;^) fory> 0, 

(142) / f{x) e ixy dx = l 

{ - 2 ^Eq zj <o re s(f(z) eizy -,Zj) fory< 0. 

Proof. If y = m > 0, then (142) follows immediately from Theorem 4.2.2. 

If y = m < 0, then we must modify the proof of Theorem 4.2.2 by reflecting Ci ,2 with respect to 
the s-axis and then change the orientation of the curve, such that the line segment on the rr-axis is 
traversed from r 2 towards — n. In the conclusion we must reverse this direction, which causes the 
change of sign in (142), second line. □ 


Example 4.2.4 A typical transition function f(z) in the Theory of Electric Circuits is given by 


JK ’ 1 + 2 iriRCz’ 

where R denotes the resistance, and C the capacity. The corresponding response function is given by 
the Fourier transform 

/ +oo -t r+oo -i 

2 ttRC 
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The only singularity of the integrand is the simple pole at Zi = 

1 k 

are constants k, r > -—, such that \f(z)\ < — for Id > r. 

2nRC \z\ 

There is no singularity in the lower half plane, hence 

/ +oo 

/( x) e 2nixt dx = 0 for t < 0. 

-OO 

If t > 0, then 

/ 


2 ttRC 


, and it is obvious that there 


h{t) = 


r+0 ° 27T7 

e 2mxt dx = 1 


f(x)i 


2tt iRC 


■ res 


_ a 7 2ttRC 
\ Z 2tt RC 


\ 1 ( 

' t \ 

1 = RC eXp [ 

, RC) 


We mention in particular the case, when f(z) is a rational function. The proof is trivial. 
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Corollary 4.2.3 Let P(z) and Q(z ) be polynomials. If 

1) the denominator Q(z ) has no zero on the real axis, 

2) the denominator Q(z) is of higher degree than the numerator P(z), 

3) the constant m is a real positive number, 

P{x) 


then the improper integral f_ 

r+ °° P(x) 


(143) 


Q(x) 


Q(x) 

dx = 2m 


res 


Sszj >0 


dx along the real axis is convergent with its value given by 

P ( Z ) pimz. 

Q{z) ’ j 


We get another useful corollary, when we assume that the analytic function f(z) has real values on 
the real axis. 


Corollary 4.2.4 Let f(z) be analytic in an open domain f l containing the closed upper half plane 
iAz > 0 with the exception of only a finite number of points z \,..., z n , none of them lying on the real 
axis. 

Assume that there are positive constants R, a, c> 0, such that (137) holds, i.e. 

|/(z)| < for ^z> 0 and \z\ > R. 

Finally, assume that f(x) £ M is real for every iff. 

Under the assumptions above the two improper integrals 

/ +oo r+oo 

f(x) cos (mx) dx and / f(x) sin(ma;) dx 
-oo J —oo 

exist for every positive constant m > 0, and their values are given by 


f+O o 

f _ I 

(144) / f(x) cos (mx) dx = < 

2tt* res(f(z)e lmz -,Zj) \ 

J — OO 

{ 0 J 

and 


f+OO | 

f _ 1 

(145) / /( x) sin (mx) dx = 3 < 

2m Y res(f{z)e lmz ]Zj) \ 

J —oo | 

{ a*f>o J 


Proof. The proof is trivial. Just split (138) into its real and imaginary parts, using that f(x) is real 
on the real axis. □ 

Remark 4.2.1 Warning! the value of e.g. f_°° f(x) cos(mx) da; is not given by 
(146) “27 t i ^ res(/(z) cos(mz) d z; Zj) 

Qzj >0 

which would be natural to expect. By using (146) one makes implicitly the error that one tacitly 
applies Euler’s formula cos mz = - e Lrnz + - e~ lmz in the integral. In the latter term we have the 
constant — m < 0, violating one of the important assumptions of (138). 0 
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z z 

Example 4.2.5 The simple poles of the function -e lz are ±i 4 R, and —^- 

z 2 + 1 z- + 1 

and m = 1 > 0, so the assumptions of Corollary 4.2.4 are satisfied. Thus, 


0 for z —> oo, 



x 

x 2 + 1 


e lx da; = 2-7ri • res 



27Ti • 


i e 1 " 1 
i + i 


7T i 

5 

e 


and we get 

/•-(-OO 


X cosx 
X 2 + 1 


/•-Too 


da; = 0 


and 


x sin x 7T 

2 , 1 da; = 
x 2 + 1 e 


where both improper integrals are convergent. 0 


4.3 Cauchy’s principal value 

We considered in Section 4.2 improper integrals of the type JJ 1 " 00 /( x) dx, where f(z) is analytic in a 
domain f 1 containing {z C C | $sz > 0} \ {z i, ..., z n }, where none of the singularities zi,....z n lie on 
the real axis. 

We shall in this section modify the concept of integral in such a way that we may allow simple poles on 
the path of integration. We shall in the follows show that we also in this case may obtain meaningful 
residuum formula. 

Poles of higher order, or essential singularities will not be allowed on the path of integration. 

We first introduce 


Definition 4.3.1 Let f be analytic in an open domain f l containing R\ {xo}, where the real number 
Xo £ R is a simple pole of f. If the symmetric limit 

e z'+oo'i 

+ / \ f( x ) dx 

J Xq+£ ) 

exists, we say that the improper integral of f(x) from — oo to +00 has its principal value given by 
(H7), and we write 

rx o—s /•+cxd'v 

f + / \ /( x ) dx. 

— OO J XQ+£ ) 


(148) pv 


/•-Too 


/(x) dx := 


lim 

;—> 0 + 




The notation “pv f^°° ■ ■ ■ dx” indicates that there is “something wrong” with the improper integral, 
though it is not worse than that we obtain convergence, if we remove a small symmetric interval 
around the simple pole xo and then let this symmetric interval shrink towards xq- 

Definition 4.3.1 can clearly be extended to the case, where f(z) has a finite number of simple poles 
on the real axis. The details are left to the reader. 

The line of integration does not have always to be the real axis. Any nice curve with a finite number 
of simple poles may be treated in a similar way. The notation will of course be the same as the above. 
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Also, a function may apparently have a simple pole at Xo £ M, where a closer examination would show 

sin z 

that the singularity is removable. One simple example is the function -, where z = 0 is a pole 

of at most order 1, and we know already that it is in fact removable. We shall treat the improper 
integral of this important function in Example 4.3.1. For the time being we only mention that if the 
principal value exists, and the apparent simple poles actually are removable singularities, then “pv” 
can be removed, 

/ +OO r + OO 

f{x) dx = fix) dx. 

-oo J —OO 

In the Real Calculus Definition 4.3.1 can be formulated more generally. We have here restricted 
ourselves to analytic functions in a neighbourhood of R. with exception of a finite number of real 
simple poles, because we want to find the values by using the Calculus of Residua. 

In order to obtain these residuum formulae we first prove the following simple lemma. 


Lemma 4.3.1 Let f{z) be analytic in a deleted open disc B (z 0 , R) \{.zo}, and assume that the centre 
Zq is a simple pole of f. Denote by 

C{e) : z = z 0 + ce l& , ©£[0,7r], 

a family of semicircles of centre Zo and redius e £ ]0, f?[. Then 

(149) lim [ f{z) dz = ni res if; zo). 

E ^°+ Jc{e) 


Proof. We put 


/(*) = 


5(A), 


z ^ zo 

where giz) is analytic in the whole disc B [zo, R ), and where a = a_i in the Laurent series expansion, 
hence a = res(/; zq). Then for e £ 




f fi z ) dz = 

Jc{e) 


'C(s) z ~ z 0 


dz 


>C(e) 


giz) dz, 


where we have the estimate 


giz) dz 


< sup<^ \giz)\ 


' G(e) 


supplied with the computation 


z £ B ( zq, - R ) • 7rer —> 0 for e —> 0+, 


'c(s) z - z 0 


dz = 


1 


e e 


T^ie e lt> d0 = lira = m res(/; zq) , 


i 0 


and the lemma is proved. □ 

From Lemma 4.3.1 we easily derive the following theorem, which also should be compared with 
Theorem 4.2.1. 
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Theorem 4.3.1 Let f : LI —> C be analytic in an open 
{z G C | S2 > 0} \ {zi,.. 

Assume that the (finitely many) singularities on the real 
constants R > 0, c > 0 and a > 1, such that 

(150) |/(2)| < for |z| > R and Zsz > 0. 

\ z \ 

Then Cauchy’s principal value vp /( x) dx is well-defined, and its value is given by 

/ +00 

f(x) dx = 2ni ^ res(f; zj) + ni ^ res(f; zj). 

‘°° SJzj >0 9z 3 =0 


domain fl which contains a set of the form 
axis are all simple poles, and that there are 


Roughly speaking, the integration line through a simple pole cuts its residuum into two equal halves, 
giving one half to the upper half plane, and the other half to the lower half plane, 


7rz res(/; zj) = 2m l - res(/; zj) 
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Figure 19: The path of integration in the proof of Theorem 4.3.1. 


Proof. The proof is almost the same as the proof of Theorem 4.2.1. The only modification is that 
we avoid all the simple poles on the real axis by small semicircles, cf. Figure 19. Using Lemma 4.3.1 
and noticing that the semicircles are traversed in the negative of the complex plane we get 

/ +oo 

f(x)dx-ni ^2 res (/;U') 

3zj =0 

and (151) follows by a rearrangement. □ 

Remark 4.3.1 In the proof above we could of course avoid the simple poles on the real axis by small 
semicircles in the lower half plane. In this case we would get 

/ +oo 

f(x) dx + 7T i ^2 res(/; Zj) = 27ri E res (/; ~?) ■ 

■°° 3zj=0 %Zj> 0 

Then note that £sZj > 0 and not just ^sZj > 0 in the latter sum, and we obtain again (151) by using 
a rearrangement. 0 


= 2tU ^2 res(/; Zj ), 

^szj >0 


Obviously, Corollary 4.2.1, Theorem 4.2.2, Corollary 4.2.2, Corollary 4.2.3 and Corollary 4.2.4 can all 
be extended to the principle value, if there are only a finite number of simple poles on the real axis. 
In each case we add half the residuum to the solution formula. It is left to the reader to formulate 
and prove these simple extensions. 


Example 4.3.1 Important! We shall here compute the important improper integral 


f*+oo 


sin a; 


dx, 


which occurs in some applications in the technical sciences. A straight on attack on the corresponding 


analytic function 


sin z 


is doomed to failure, because 


sin z 


does not satisfy an estimate of the type 


(150). Furthermore, we only integrate along the positive real axis, so the problem apparently does 
not fit into the present theory. 
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The trick is the following: Apply Euler’s formulas to get 


sin x . 1 

-dx = — 

x 2 * 


dec = 



dx. 


This implies that if the right hand side has a well-defined limit for e —> 0+ and R —> +oo, indepen¬ 
dently of each other, then the left hand side is also well-defined, and 


+ °° sin a; 1 /' + °° e“ , 

-dec = — vp / — dx. 

x 2i x 


We shall only prove that vp f_ 


+oo e- 


dx exists, i.e. we shall check the singularities of the integrand. 


The only (finite) singularity of the analytic function — is the simple pole at zq = 0. It lies on the 

z 

path of integration, so it contributes to Cauchy’s principal value with the amount 


ni res —; 0 = ni lim e lz = ni. 

1 x ) 2—>o 


Furthermore, 


= - e llz , where m = 1 > 0 and 


= t pr for 2 7 ^ 0 , so a = 1 > 0 , and it follows 

r N- 


from Theorem 4.2.2 that the limit R —> +oo will not cause any trouble either. Hence we conclude 
that the improper integral J 0 +oc -—- dx is convergent, and its value is 


f + OO 


sin x 1 

-dx = — pv 

x li 


0 

+oo 

—oo 


X 

e ix 1 ( e 

— dx = —7n res 
x 2 1 


Z ’° = 2 ' 


0 


4.4 The Mellin transform 

We shall in this section consider improper integrals of the form 

r + °° dr 

(152) / f{x) x a —. 

Jo x 

When these integrals are considered as functions of a, we get the Mellin transform 9Jt{/}(a) of the 
function /. This is closely related to the two-sided Laplace transform. In fact, if the integral (152) is 
absolutely convergent, then we get by the change of variable, x = e _t , 

r + °° dr r + °° 

(153) M{f}{a) := / /(x)x Q — = / /(e ‘) e at dt, 

Jo x J -oo 

which is recognized as the two-sided Laplace transform of the function g{t) = f(e~ t ) at the point a. 

We shall not go into the applications of the Mellin transform, or the two-sided Laplace transform. 
The purpose of this section is only to use Complex Functions Theory to compute some integrals of 
the form (152). 

In order not to make the theory too complex we have restricted ourselves to real ael. It is possible, 
though far from trivial, also to allow a to be a complex number, but this would require a discussion of 
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the so-called many-valued, functions, which we have postponed to Ventus: Complex Functions Theory 
a-3, and even with this discussion, the generalization becomes far from trivial and only of interest of 
the few. 

In order to combine (153) with Complex Functions Theory we also must require that f(z) is analytic 
in an open domain fh More precisely, f(z) is analytic in all of C with the possible exception of only a 
finite number of singularities, none of which is lying on the positive real axis R + . In the literature one 
mostly assumes that the singularities are all poles, but except for 2 = 0 , which must at most be a pole, 
the proof below shows that all the other singularities are allowed to be even essential singularities. 

Then we shall fix the meaning of the factor x a . If a G Z, then the definition of x a is straightforward, 
and we just apply the methods of Section 4.2. We therefore assume in the following that a G R. \ Z. 

The next problem is to define the analytic power function z a for «gR\Z. first put 

(154) Log 0 2 := In \z\ + i Arg 0 z, Arg 0 z G ]0, 2tt[, 

for 2 G fli := C \ {R+ U {0}}. Note that Log 0 is not the principal logarithm, Log, and Arg 0 is not 
the principal argument, Arg. They of course agree in the open upper half plane, but they are different 
in the open lower half plane. 

Then apply (154) to define z a by the “obvious” formula 

(155) z a := exp(aLog 0 2 ) for z G fli = C \ {R + U {0}} . 

Then z a is a composition of analytic functions in fii, so it is also analytic in fii, and a routine check 
shows that 

- 7 - z a = exp(aLog 0 2 ) • - = az a ~ 1 for z G fli- 

dz z 

Furthermore, since a is real, it follows from (155) that 

(156) |,£ a | = | exp(a{lnr + z©})| = r‘ G 1, 

a formula, which is not true, if a is complex, cf. Ventus: Complex Functions Theory a-3. 

After these preparations we can formulate 


Theorem 4.4.1 Let f be analytic in 12 = C\ {z\, ■ ■ ■ , z n }, where Zj (£ R+, j = 1,..., n. Assume that 
there exist constants a, (3 G R, where a < (3, and C , Rq, Vq G R+, such that 

(157) \z a f(z)\ < C for\z\<r 0 , z G Cl, 

(158) 1 2 ^/( 2 ) | < C for 1 2 1 > R 0 , 2 G 0. 

Then the improper integral / 0 + °° f{x)x a — is convergent for every a G }a, /3[\Z, and its value is given 
by the residuum formula 

(159) f f(x)x a — = - n LXp ( ma ^ res(f{z)z a ~ 1 -, z 3 ) , for a £]a, (3[\Z. 

Jo x sin 7 ra 
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Figure 20: The path of integration in the proof of Theorem 4.4.1. 


Remark 4.4.1 One usually assumes that also a, (3 £ R+. The proof below shows that this is not 
necessary. Formula (159) does not make sense for a £ Z, because then the denominator sin 7 ra = 0. 
As mentioned earlier we may use the theory of Section 4.2 instead, or a limit process on (159). 0 


Proof. We shall as usual find a convenient path of integration Cr,t,v in fi' = f2 \ {R + U {0}}. We 

7r 

on 


choose R > R 0 and r < r 0 and v £ 


°’2L 


, from which we define the integration path Cu,r,v 
Figure 20 composed of two circular arcs of centre 0 and two line segments on the lines through 0 
forming the angles ±v with the positive rc-axis. 


Assume that Cr^,v has been chosen, such that all singularities ^ 0 of f(z)z a 1 lie inside Cr^,v Then 
we get by Cauchy’s residuum formula 

(160) 2ni res(f(z)z a ~ 1 ; Zj) = </ f(z)z a ~ 1 dz = f + f + f + f f(z)z a ~ 1 dz, 

ZjZjLO ^Cr,t,v J — C v J — Cr J c v 

where C denotes a circular arc oriented in the positive sense of the plane seen from 0, and c denotes 
a line segment oriented in the direction away from 0 . 


First note that 


I, 


f{z)z 


dz 



e (°-l)(ln t+iv) e iv ^ 


f(t e iv ) e {a ~ 1)iv ■ e*' 


■t“ _1 dt ■ 


R 


f(x)z a ~ L dx 


Jr Jr 

because the integrand is continuous in the closed bounded interval [r, R ], so 
the limit under the sign of integration. 


we 


for v —► 0 +, 

are allowed to take 


Since the segment c_„ has the parametric description 
z(t) = t ■ e^ 2n ~ v \ t £ [r, R] and 27 t — v £ ]0, 2n[, 
we get 

Logg z (t) = In t + i{2 tt — v) on c_„, 
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hence analogously, 

[ f(z)z a ~ 1 dz = - f f(z)z a ~ 1 dz = + [ R f(te- iv )t a - l e {a - 1){ - 2vi - iv) e i ^- v) dt 

J —c_ v J C— v J r 

pR pR 

= / f(te- iv )t a ~ l e ai ^- v) dt^ / f{x)x a ~ 1 e 2nio dx for v -> 0 + . 

Jr Jr 

Thus, for fixed r and i?, 


lim 

v— >-0+ 


+ j | f{z)z a ~ l dz = £ f(x)x a ~ 1 {1 - e 2 ™} dtr 


= -2 i e l 


gi7ra g —iira pR 


r R 


2 i 


f(x) a ~ i dx = -2ie™ a sm(na) / f(x)x a ~ L dx. 



Brain power 


TO 
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Next we turn to the estimates of the integrals along the two circular arcs Cr and C r . Here we use 
(156), (157) and (158) to get 


[ f(z)z a ~ 1 dz = [ z p f{z)z a - p ~ l dz 

Jc R Jc r 


< 2t tR ■ C ■ R°- /3 " 1 


= 27tC • IT (/3 " q) -> 0 for R —> +oo, 


and 


[ f(z)z a - 1 dz = [ 

JC r JC r 


z a f(z)z a - a - 1 dz 


< 2nr ■ C ■ 


= 2t tC ■ r a ~ a ^ 0 for r —> 0+, 


where both estimates are independent of v £ 


°’ 2 L 


When we apply the three limit processes v —> 0+, r —> 0+ and R —► +oo, and (160), we conclude that 
the improper integral / 0 + °° f(x)x a ^ 1 dx is well-defined and that 


r+oo 

2m y res (f(z)z a ~ 1 \Zj) = — 2i e I7ra sin7ra / f(x)x a ~ 1 dx, 

Zj^O Jo 


and (159) follows by a rearrangement. □ 


Example 4.4.1 The function f(z) = 


z 2 + 1 


, z £ C\{—*,*}, satisfies the estimates 


\z°f{z)\ < 2 for \z\ < r 0 , |z 2 /(~)| < 2 for \z\ > R 0 , 

for some constants 0 < tq < 1 < i?o, and none of the simple poles ivi lies on R + . It follows from 
Theorem 4.4.1 for a G ]a, (3[ = ] 0, 2[ and a ^ 1 that 


r+00 ^ 


dx 


x 2 + 1 x 


7 r e 


sm 7ra 


res 


~a— 1 


z 2 + r 


; i + res 


v a— 1 


£ 2 + 1’ 


7re 


sm 7r a 


{s exp ((“- 1 ) +)-S exp (<““ 1 )i T)} 


^g-iira 

i 

f 1 I 

f. 7T\ 


' o < 

t exp 

, * a tt ) 

sin 7ra 

2 

i* 

< 2/ 


f 3tt\ 

1 a — • 1 

v 2 


1 


sin 7 ra 2 


jexp(-ja|) + exp(ia|)j 


7TCOS 


(+) 


sm 7ra 


2 sin 


in H) 


For a = 1 we get straightforward 


f*+oo 


x 1 dx 


/•+00 


X 2 + 1 X 


dx 7T 

—--= — = iim 

X 2 + 1 2 a — ► ! 


2 sin ( a — 


K)' 


132 


Download free eBooks at bookboon.com 



























Calculus of Residua 


Applications of the Calculus of Residua 


so we have in general, 


f + °° *° _1 H 7r 
Jo ^ 2sin (a\) 


foraG]0,2[. 0 


4.5 Residuum formulae for sums of series 

It is also possible for a large class of convergent series to compute their sum by using a residuum 
formula. We first prove 


Theorem 4.5.1 Let f : Lt —> C be analytic in Lt = C \ {z i,..., Zk}, where Zj Z for j = 1,..., k. 
Assume that there exist constants R, c > 0 and a > 1, such that 

(161) \f(z)\ < for \z\ > R. 

Then the series Y^n =-oo /( n ) convergent, and is sum is 

+oo k 

(162) Y f ( n ) = ~ 7T Y res{cot{TTz) • /(z); zj). 

n =—oo j—1 


Proof. Since a > 1, and Zj ^ Z for j = 1 it follows that Y^n^-oo /( n ) ^ ias a convergent 

majoring series, e.g. 1——7 f° r some constant ci, so it is itself convergent. 

\n\ +1 

Then introduce the auxiliary function g{z) by 
g(z) := 7r • cot(7 rz)f(z). 


Since Zj £ Z, j = 1,..., k, we see that cot(7r Zj) is well-defined for all j = 1,..., k. Therefore, the 
singularities of g{z) are {zi,..., Zk} U Z, where the Zj are the singularities of /, and n £ Z are simple 
poles stemming from cot(7rz). 


We compute 

res (g;Zj) = 7rres(cot(7rz)/(z); Zj ), 

C0S(7TZ) 


(163) 


res (g, n ) = n ■ 


fz sin (^) 


f(n) = f (n), 


2=1 


n £ Z. 


Choose for every N £ N the path of integration Cn as shown on Figure 21. If iV > 42, then the set 
(zi, ..., Zk} lies inside Cn, so we get by Cauchy's residuum theorem, 



g{z) dz = 27 rf 


+N k 

Y fW + Y 7r • res(cot(7rz)/(z); Zj) 
i=i 


n——N 
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-< - 

(N+%o) 






-(N+%.) 



N+%«)v 






- >- 

-»- 







Figure 21: The path of integration Cm in the proof of Theorem 4.5.1. 


and therefore by a rearrangement, 

+N k 

(164) y f(n) = —7r ^ res(/(z) cot(7rz); Zj) + 

n=—N j —1 

The left hand side of (164) converges for TV —> +oo towards the sum of the series. 


1 

2m 


- <b g{z)dz. 
1 J Cn 


Then notice that it follows from a result in Ventus: Complex Functions Theory a-1 that 


| COt(7T2)| 2 = 


C0S(7T3) 


sin(7T2i) 


cosh 2 (7ry) — sin 2 (7ra;) 
COSh 2 (7T7/) — COS 2 (7T3) 


Since x = ±{TV+i}on the vertical segments of Cm, we get on these, 
, 2 / 


| COt(7T3)j“ < 


2 cosh (Try) 


cosh (7ry) 


= 1, 


so cot (773) are bounded on the vertical segments. 

We have y = ± {TV + on the horizontal segments of Cm, so 


| COt(7T3)| 2 < 


cosh 2 (- 7 T {TV + t}) 


= 1 + 


< 2 , 


cosh 2 (7T {TV + i}) — 1 cosh 2 (7T {TV + i}) — 1 

because cosh 2 (7r {TV + t}) > 2 for TV £ No- 

It follows from the assumption (161) that z f(z) —> 0 for 3 —> oo. We can therefore to every e > 0 

£ 

find N e £ N, such that z f(z)\ < — for all n, for which \n\ > N e . Then we get the following estimate 
for TV > TV e , 


7Z— g{z)dz 


27 vi 


Cn 


1 

27r 


z f(z) 

J Cn 


7r cot ( 773 ) 


d3 


< — { - -7T-2 
27t 1 8 


TV- 


_ 4(2TV + 1) 

Cn) ~ 8 (TV + |) ’ 
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and it follows that 


lim - ® q(z) dz = 0. 

N—>+oc 2m J Cn 

The theorem follows by taking the limit N — > +oo in (164). □ 

P( z ) 

Theorem 4.5.1 is mostly applied when f(z ) = ——- is a quotient of two polynomials with a zero of 

Q{z) 

at least order 2 at oo, i.e. degree Q— degree P > 2. The remaining necessary requirement is that 
Q(n) yf 0 for all n £ Z. 

Note also that if Zj ^ Z is a simple pole of f(z), then 
(165) res(cot(7rz)/(2;); Zj) = cot(zjn) ■ res(/; Zj ). 

If in particular 2 Zj C Z is an odd number , then becomes a removable singularity of 
g(z) = cot (7 Tz)f(z), so its residuum is 0. This is in agreement with the fact that the right hand side 
of (165) is 0, because then cot(^-7r) = 0. However, be aware that (165) is not true for poles of order 
> 2 . 
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Example 4.5.1 A simple example is given by 


/(*) = 


1 


{*-§r 


for z£C\|^. 


Then \ f(z)\ < — 7 ^ for \z\ > 1 and some constant c > 0. Since z 0 = h d Z is the only pole, we can 

M 2 

apply Theorem 4.5.1. Here, z 0 = \ is a pole of order 2 for f(z), i.e. it is a pole of at most order 2 for 
the auxiliary function 

COt(7T3) 

{--I} 2 ' 

(The order is of course 1.) If we choose q = 2 in Theorem 3.7.1, we get 

( cot(7rz) l\ 1 d ( f n 2 , 1 d . 

"vbATV = (^I)T 2} = 

= lim (— {l + cot 2 (7T^)} 7r) = — 7 r. 

Then the sum is given by Theorem 4.5.1, 


+ OO 


E =- 


7r • res 


COt(7T2:) 1 


n=— 00 




r! - = 7T . 


Using that 

+ OO 


1 +OO -1 +OO 1 +OO A 

E r _m2 = E r _m2 +E r ,,2 = 2 E ( 2n + 1)2 = ^ 

n=—oo \^ 2 1 n =1 2 1 n—0 2 J n=0 ' 

we have found the sum of the following important series from the Theory of Fourier series 


+oo 

E 


71=0 

Then also 

+oo 1 


(2 n + l) 2 8 ' 


+oo 


E n 2 E 


1 


1 


+oo 




1 


+oo 

E 


71=1 


n^o (2n + 1)2 22 (2n + 1)2 {2 2 } 2 ^ 0 (2n + l) 2 

+oo J ! 


1 1 1 
- < ! 1 +4 + ^2+^3 


E 

71=0 


(2n + l) 2 1-i 


6 ’ 


which is also well-known from the Theory of Fourier series. 0 


1 c 

Example 4.5.2 Let f(z) = -=■, where a > 0 is a constant. Then 1 /( 2 )I < 7 —^ for I 2 I > 2 a for 

z z + a z ' \z\ z 

some constant c > 0, and f(z) has only the two simple poles z = ±ia ^ Z, so we can apply (165). 

First compute the residua, 


res 


1 


z 2 + a 2 ’*/ 2 ia 


1 


and res 


1 


1 2 ’ 2 ia' 
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Then it follows from Theorem 4.5.1 that 

+ OO 


E 


i 


n=—oo 


n 2 + a 2 


= —7r i — cot(*7ra) H-— cot(—Z7ra) 


2 ia 


—2 ia 


7r cos(i7ro) 7r cosh(7ra) n , . . 

— ^“77-7 = • 7-T-T7-77 = ~ ' COth(TTa), 

za sm(z7ra) za z • smn(7ra) a 


so 


+oo 

E 


1 


J_ 1 

- 9«2 + 9 2^ 


n 2 + a 2 2a 2 ' 2 

n=0 n=—oo 


1 1 7T , . . 

= TTo + 77— COth(TTo). 


n 2 + a 2 2a 2 2a 

We get in particular for a = 1 that 


+oo 


1 1 7T 

y „ . , = - + — coth 7r. 


n—0 


n 2 + 1 


0 


Theorem 4.5.2 Let f : Ll —> C be analytic in C\{,Si,... ,Zk}, where Zj Z for j = 1,..., k. Assume 
that there exist constants R, c > 0 and a > 1, such that (161) holds, i.e. 

\f(z)\ < -E f°r\z\>R. 

M° 

Then the series Y(,n=-oo(~ l) n /( n ) convergent, and its sum is given by the residuum formula 


(166) 


+oo 

E (-!)"/(«) 


n=—o o 



Proof. (Sketch) The proof is trivial modifications of the proof of Theorem 4.5.1, where we replace 

cot(7T2:) by In (163) we here get 

Sin(7T2) 


res 


4*1;») 

sin(7rz) / 


(-1 )"/(«), 


so the left hand side of (164) is replaced by X)E_Ar( — l)"/( n )- Finally, the estimate of p—^—yy 011 
Cn is much easier than the estimate of | cot(7rz)| above. □ 


Example 4.5.3 Let f(z) = —^, where a > 0 is a constant. We have already shown in Exam- 

z z + a 1 

pie 4.5.2 that the assumptions of Theorem 4.5.1, hence also of Theorem 4.5.2, are fulfilled. Hence, 


y M) ra - J 1 1 

n 2 + a 2 I 2 ai sin(i7ra) 

no v x ' 


-1 1 
2 ai sin(—77ra) 


Using that 

(-!)"" _ (-1)" 
(— n) 2 + a 2 n 2 + a 2 


for n £ Z, 


7T 1 7T 1 

ia sin(*7ra) a sinh(7ra) 
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we get 

y? (-1)" 1 y? (-1)" _ J_ n_ 1 

n 2 + a 2 2 a 2 2 n 2 + a 2 2 a 2 2 a sinh(7ra) 

n =0 n =—oo 

Then in particular for a = 1, 

V i—^ = 1 + —-—. 0 

^ n 2 + 1 2 2 sinh 7r 

n=0 


Finally, it should be mentioned that even if we are now able to find the exact sum of a lot more 
series than by Real Calculus alone, there are still many series that cannot be treated in this way. The 
simplest example is perhaps 



n= 1 


1 . 202 , 


and more generally, Y^n=i 


1 

n 2p+1 


p £ N, for which no exact formula is known. 



138 


Download free eBooks at bookboon.com 




Click on the ad to read more 













Calculus of Residua 


Index 


Index 


accumulation point, 9, 42 
annulus, 73 

basic power series, 24 

Bessel function, 112 

biharmonic equation, 61 

binomial coefficient, 25 

boundary value problem, 60, 64, 66, 69 

calculus of residua, 84 

Casorati-WeierstraB’s theorem, 87 

Cauchy Integral formula, 7 

Cauchy Integral Theorem, 7 

Cauchy multiplication, 18,28, 111, 115 
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difference equation, 31 
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Dirichlet problem, 7 
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disc of convergence, 12 
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essential singularity, 87, 88, 93, 95 
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Euler’s formula, 123 

existence and uniqueness theorem, 29 
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exponential, 33, 36 
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Fourier series, 44, 66, 71, 79, 136 
Fourier series expansion, 66 
Fourier transform, 121 

geometric series, 15, 22, 37, 76 
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harmonic function, 55 
harmonic functions, 7 
heat equation, 69 
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improper integral, 113 
improper integrals, 7 
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inspection, 32 
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Laplace transform, 7 
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Laplace equation, 55, 69 
Laplace transform, 45, 71 
Laurent series, 7, 40, 45, 71, 72 
Laurent’s theorem, 74, 84, 90 
level curve, 68 
limes inferior, 10 
limes superior, 9, 10, 12 
linear differential equation, 80 
linear differential equations, 29 

many-valued function, 129 
maximum principle, 48 

maximum principle for harmonic functions, 59 

mean value theorem for analytic functions, 60 

Mean Value Theorem for Harmonic Functions, 65 

mean value theorem for harmonic functions, 60 

Mellin transform, 7, 128 

method of inspection, 31 

method of power series, 31 

minimum principle, 49 

multi-valued functions, 7 

multiplicity, 41 

non-isolated essential singularity, 93 
order, 41 

Parseval’s equation, 47, 48 

Parseval’s formula, 79 

Phragmen-Lindeldf’s theorem, 49 

Picard’s theorem, 87, 89, 92 

Poisson’s Integral Formula, 64 

Poisson’s integral formula, 64 

polar coordinates, 45 

pole, 85, 93, 97 

pole at °°, 91 

pole of order 0, 97 

polynomial, 8, 92 

potential, 67 

potential theory, 55 

power series, 9, 11 

principal value, 25, 124 

radius of convergence, 10, 12 

rational function, 75, 98, 112, 117 

real Taylor series, 24 

recursion formula, 31, 35, 80 

removable singularity, 84, 96, 98, 108, 109, 125, 
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residue calculus, 7, 45 
residuum, 93 
residuum at °°, 102 

residuum formula for the Fourier transform, 121 

residuum formulae for sums of series, 133 

response function, 121 

Riemann surface, 40 

Riemann surfaces, 7 

Schwarz’s lemma, 52 

sequential sequence, 11 

simple pole, 86, 98-100, 109, 116-118, 124-126, 
133,135,136 
simple zero, 92 

singular point for dierential equation, 80 
solution procedures for solving linear differential 
equations, 31 
stability, 7 

Stirling’s formula, 33 
Stokes’s theorem, 67 

Taylor series, 7, 20, 23, 76 
temperature field, 69 
trigonometric integrals, 7, 111 
two-sided Laplace transform, 128 

uniform convergence, 16,19 

weak Phragmen-Lindeldf’s theorem, 51 
Weierstra’s double series theorem, 27 
work, 67 

3-transform, 7, 45, 71, 91 
zero, 42 
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